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ABSTRACT 

This thesis deals with the study of on-line measurement 
and control problems associated with the solution 
polymerization of methylmethacrylate in a batch reactor. An 
experimental reactor system with facilities for on-line 
measurements of polymer properties such as density, 
viscosity, torque on the stirrer and molecular weight 
distribution has been built as part of this work. A. process 
control computer was used for on-line data acquisition and 
implementation of control policies on the reactor. An 
automatic sampling system that collects the samples of the 
reaction mixture periodically from the reactor and injects 
them into a size exclusion chromatograph (SEC) has also been 
developed. Computer programs were developed for automatic 
SEC data acquisition and data processing of this information 
to calculate number and weight average molecular weights. 
The analysis of the meaSurement data and demonstration of 
the correlation between density and conversion showed that 
the on-line densitometer can be used to give an indirect 
indication of conversion during the course of 
polymerization. The SEC was used as an on-line instrument to 
measure molecular weight averages. 

A mathematical model of the reactor was developed and 
parameters of the model were estimated using the 
experimental data. Verification of the mathematical model 
was shown by the excellent agreement between model and 


experimental data. Open loop optimal control policies were 
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developed by the application of state-space theory to the 
model of the reactor. These policies were implemented on the 
reactor system. The results show that experimental and model 
conversion results were in good agreement whereas 
discrepancies of the order of 25-30 percent were observed 
for number and weight average molecular weight. Optimal 
feedback control policies were also derived and implemented 
on a simulated batch reactor. The feedback controller 
tracked the nominal conversion and molecular weight profiles 
in the presence of disturbances in initiator concentration. 
This study demonstrated that with the development of 
proper on-line instrumentation, advanced computer control of 
polymer properties such as conversion and molecular weight 


distribution is possible. 
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Arrhenius factors for the terms k, to ke, 


Arrhenius factors for dissociation,propagation 
and termination reaction, respectively. 


3x3 state matrix in linear model, 


constants in viscosity correlation. 


Sx t-cOontrol matrix, in linear model. 
Ort tok / ke) 


ratio of monomer transfer constant to 
propagation constant. 


ratio of solvent transfer constant to 
propagation constant. 


activation energy for the terms k, to ke, 
cal/mol. 


activation energy for dissociation, propagation, 
and termination respectively, cal/mol. 


error vector in parameter estimation. 
initiator effiency. 

3x3 state matrix in discrete model. 
3x1 control matrix in discrete model. 


Hamiltonian 
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3x1 disturbance matrix in discrete model. 
initiator concentration, mol/L. 
initial initiator concentration, mol/L. 


nominal value of the initiator concentration, 
mol/L. 


kpY2£ka/k, = A,exp(-E,/RT) 

kmV/2£ky/k; = Azexp(-E2/RT) 

ksV¥2tka/k;y = Azexp(-E,/RT) 
2kp?k,/V2tkak,':*® = Agexp(-E,/RT) 

2kp7ku7v 2tkyky = v= Agexp(-E,/R?) 
(24+p) kp? /keb= Asexpt-Ee/RT) 

dissociation rate constant, 1/s. 

initiation rate constant, 1/s. 

transfer to monomer rate constant, L/mol-s. 
propagation rate constant, L/mol-s. 
transfer to solvent rate constant, L/mol-s. 
kink ejmb/mol-s. 


termination by combination rate constant, 
L/mol-s 
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Kea termination by disproportionation rate constant 


L/mol-s. 

M monomer concentration, mol/L. 

MW molecular weight of the monomer. 

M, number average molecular weight. 

Me desired number average molecular weight. 

M, true molecular weight of the polymer. 

Mw weight average molecular weight. 

PD polydispersity. 

EP, Ps co-state variables. 

Peery concentration of dead polymer of chain length 
x and y respectively, mol/L. 

Q weighting matrix for states. 

R gas constant, cal/mol-°K. 

R weighting matrix for control. 

S solvent concentration, mol/L. 

t time, Ss. 

et final time, s. 

Ay temperature, °K. 

Vv elution volume, mL. 
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volume of reaction mixture, mL. 


W 3x1 disturbance matrix in linear model. 

wW 3x3 weighting matrix in parameter estimation. 

X conversion 

Pe desired conversion. 

u control vector. 

x state vector 

> = measured values of output vector. 

Yp predicted values of output vector. 

Xo, At, Az zeroth, first and second moment of the live 
radical distribution respectively. 

Mo, 1, M2 zeroth, first and second moment of dead polymer 
distribution respectively. 

p Reece 

n viscosity of the polymer solution. 

PmrPs,Py density of monomer, solvent and polymer 
respectively, g/mL. 

®, volume fraction of solvent. 
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1. INTRODUCTION AND OBJECTIVES 

Currently some 100 million metric tons per year of 
Synthetic polymers are produced in the world in a wide 
variety of polymerization reactors. The total world sales of 
raw polymer is approaching $100 billion per year and sales 
of the final product after processing, molding and 
compounding is many times this figure. Hence the field of 
polymerization reaction engineering is of significant 
economic importance. 

Until a decade or two ago, "polymers were mainly 
manufactured in batch reactors from faithfully prepared 
recipes scaled up from the chemist's beaker" (Ray(1983)). 
However with the growth of demand and increased price 
competition more efficient polymerization methods have been 
required. Manufacturers of high volume polymers are moving 
to fewer product lines, more uniform product and the use of 
continuous reactors. Similarly producers of low volume, high 
quality and high value polymers are realizing that their 
competitive edge comes from a deeper understanding of the 
relationship between polymerization conditions and _ the 
product quality. 

Polymerization reaction engineers have been developing 
improved and efficient processes that can significantly 
reduce capital expenditure as well as operating costs. An 
example is the low pressure fluidized bed process of Union 
Carbide to produce low density linear polyethylene which can 


reduce the plant capital expenditure by 50% and energy 
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consumption by 25% of the conventional high pressure 
polymerization process(Ray(1983)). 

Significant economic advantages can also be obtained by 
operating the polymerization reactors in an optimal manner 
such that consistent product quality is maintained. Since 
any existing process can be operated in an efficient and 
optimal manner thereby reducing the operating costs, more 
attention has been paid to the optimum operation of the 
polymerization processes recently. As a result many 
operations in the polymer production have been automated and 
many of polymerization process variables such as temperature 
and pressure are being automatically controlled with the use 
of computers. The state of art of control of polymerization 
reactors has been reviewed by many authors (Amrehen(1977), 
Hoogendoorn and Shaw(1980), MacGregor et al.(1983)). 

The benefits of automation and computer control of the 
polymerization reactors are many. Some of them are: 

e Saving of manpower. 

ry Better control of quality of product. 

> Shortening batch time in batch reactors. 

Fs Minimum production of off-grade product in the product 
changeover in a continuous reactor. 

= Optimum operation of plant through optimization 
programs. 

= On-line measurement of process variables. 

¥ Logging for process development and trouble shooting. 


re Process alarms and automatic shutdown and hence better 
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safety. 

Inspite of these benefits and much more, the progress 
in automatic control of the polymerization processes has 
been slow compared to other chemically simpler reaction 
systems such aS ammonia synthesis, ethylene production, etc. 

The slow progress in automatic control of 
polymerization processes has been attributed to the lack of 
development of on-line measuring devices. Hoogendoorn and 
Shaw(1980) pointed out the viscous nature of the 
polymerization mixture as the main contributing factor for 
the difficulty of on-line measurements. Polymer quality has 
been traditionally controlled by measuring the properties of 
the polymer product off-line and rejecting the off-grade 
product completely or blending it with good product such 
that the blended product has desirable properties. 

In the past, on-line measurements of polymerization 
Systems were mainly limited to temperature and pressure. 
With the advent of process control, especially due _ to 
digital computers, the impetus to determine product quality 
on-line has increased markedly. In response to this, there 
have been a number of attempts to develop new sensors to 
measure various polymerization process variables on-line and 
to modify the existing instruments so that they can be used 
for on-line measurements. 

The availability of the process control computers and 
the recent development of on-line measurement techniques 


have considerably increased the scope of on-line control of 
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product quality in polymerization reactors. Most of the 
earlier studies, particularly in academia, were limited to 
the theoretical prediction of batch reactor temperature 
policies, semi-batch initiator or monomer feed rate policies 
for the optimization of product quality in polymerization 
reactors. However, the success of these optimization efforts 
depends very much upon having valid dynamic models of the 
physical and chemical processes occurring in the complex 
polymerization systems. Hence considerable research has also 
been directed towards development of mathematical models to 
describe the polymerization processes under various 
conditions. 

Optimization of product quality through feedback 
control can also be implemented on polymerization reactors. 
However the success of the optimal policies through feedback 
control are limited by the availability of on-line polymer 
characterization instruments. In recent years there has been 
considerable research in the area of on-line measurement and 
control of polymerization reactors. In fact the research is 
advancing soerapidly thatit Ys anticipated that the 1980-s 
will see the implementation of higher level optimization and 
computer control schemes throughout the polymerization 
industry(MacGregor et al(1983)). 

From the above it is clear that there is great 
incentive to study the problems associated with real time 
computer applications of optimal control policies to 


polymerization reactors. The present work was undertaken in 
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response to these needs i.e. to develop techniques to 

measure, on-line, important process variables such as 

conversion and molecular weight distribution and to 
implement computer control on a polymerization reactor. 

Specifically the objectives of this study are summarized as 

follows: 

1. Development of analytical techniques to measure on-line 
the density, viscosity, and molecular weight 
distribution of the polymer during the polymerization 
process. 

2. Evalulation of these techniques with respect to their 
accuracy in estimating, on-line, important variables of 
the process such aS conversion and molecular weight 
averages (number and weight average molecular weights). 

3. Development of a reliable mathematical model of the 
process and estimation of the kinetic parameters of the 
model from the experimental data. 

4, Application of optimal control theory to develop open 
loop and closed loop control policies to optimize the 
operation of the polymerization reactor. 

5. Implementation of these optimal policies on a laboratory 
scale reactor using a real time process computer. 

The above stated objectives were attempted on aé_ free 

radical solution polymerization of methylmethacrylate in a 

batch reactor. However the on-line measurement and control 


techniques developed in this work can be eaSily extended to 


many other polymerization systems. 
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This thesis has been divided into eight chapters each 
of which closely follows the chronological development of 
the study. Chapter two contains a description of the 
experimental system which includes a five litre reactor, an 
automatic sampling system, a real-time process control 
computer and a number of measuring devices such asa 
densitometer, a viscometer and a size exclusion 
chromatograph. 

In chapter three on-line and off-line measurement 
techniques used to monitor conversion and molecular weight 
distribution during the course of polymerization are 
described. A computer program developed for automatic data 
acquisition and on-line interpretation of size exclusion 
chromatograms (SEC) to determine molecular weight averages 
of the polymer is also described. 

In chapter Four, extensive experimental results 
obtained using the batch reactor system are presented. The 
effect of initiator, temperature and volume fraction of 
solvent on conversion and molecular weight is investigated. 

A detailed mathematical model of a free radical 
solution polymerization of methylmethacrylate in a _ batch 
reactor is derived in chapter five. The model parameters are 
estimated from the experimental results obtained from the 
reactor system. 

In chapter six open loop control policies are developed 
by applying optimal control theory to the mathematical model 


of the reactor derived in chapter five. These policies are 
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then implemented on the laboratory scale reactor and the 
experimental results are compared with the theoretical 
predictions. 

Chapter seven deals with the closed loop control of the 
polymerization reactor. The nonlinear model of the reactor 
is linearized and optimal feedback control policies are 
derived by working with a linearized model of the reactor. 
The closed loop policies are evaluated by simulation. 

Finally in the eighth chapter the overall conclusions 
and the possible extensions of this work are presented. 

This study is intended to provide the groundwork for a 
series of application studies on on-line measurements and 
real time computer control of polymerization reactors. It is 
hoped that the success of these studies and the knowledge 
gained during this research will have direct application to 


the control of industrial polymerization processes. 
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2. EXPERIMENTAL REACTOR SYSTEM 


2.1 Introduction 

As stated in chapter one, most studies on control of 
polymerization reactors have been limited to theoretical 
predictions of optimal policies required to produce a 
polymer product with desired properties in minimum time. It 
is true that very little work has been done on actual 
implementation of these control policies to a laboratory 
scale experimental reactor system. Chan and Huang (1981) 
first reported the implementation of optimal initiator 
policies on a 300 ml glass reactor for the case of solution 
polymerization of styrene. However their experimental system 
was a Simple one and did not include facilities for on-line 
measurements and interfacing with a process computer. There 
is a great research incentive to study the problems 
associated with the real-time computer implementation of 
optimalrcontrol policies -ototta® laboratory “oscaléstreactor 
system. Therefore a novel reactor was designed and built. 
This system consists of a five litre glass reactor, an 
automatic sampling system, a process control computer and a 
number of measuring devices such as a densitometer, a 
viscometer and a size exclusion chromatograph. In this 


chapter we describe in detail the experimental reactor 


system. 
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2.2 Reactor System 

A schematic representation of the experimental set-up 
is shown in Figure 2.1. It is designed in such a way that 
different measurements from the system can be used for a 
variety of studies in polymerization reactions. 

A detailed list of all the instruments and measurement 
devices used in the experimental system along with the 
technical specifications is given in Appendix A. 

adhe, reactor is a 5b) jacketed cylindrical glass vessel 
of approximately 152 mm inside diameter and 280 mm height. 
It is calibrated and marked to show the volume of reaction 
mixture inside. The top of the reactor consists of a flared 
glass flange which is used to Support the reactor by firmly 
gripping the glass flange by an outer stainless steel flange 
uSing a rubber gasket. The whole reactor assembly is 
Supportd through the outer stainless steel flange to a 
framework consisting of four inclined steel legs and a 
circular plate at the top that is welded to the stainless 
steel flange. The legs are firmly bolted to the floor of the 
laboratory, so that vibrations due to the stirreremotor dre 
kept to the minimum. 

A stirrer-motor assembly is mounted on the reactor to 
mix the reactants thoroughly. The stirrer has two turns of 
helical blades and is supported by a bearing at the centre 
of the lid. It is turned by a totally enclosed DC motor 
which is controlled by a controller that provides filtered 


DC power to the motor, changing the voltage as necessary to 
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keep the speed constant. Because of the action of the 
controller the speed of the stirrer is maintained constant 
in spite of the increase in the load on the stirrer due to 
increasing viscosity of the reacting mixture during the 
course of polymerization. The torque on the stirrer is 
available for on-line data acquisition as a 0-10V output. It 
is found that the helical blades attached to the stirrer 
provided good mixing of the reaction mixture when operated 
at 750 RPM. Also the bearing and the motor support is good 
enough to provide stirring with minimum vibration. 

The reaction mixture is continuously heated by 
Circulating hot water through the jacket using a thermally 
protected centrifugal pump. The hot water flows through a 
rotameter and 1/2" diameter copper tubes. The constant 
temperature hot water tank has a capacity of 25 litres and 
is provided with agitators to keep the temperature uniform 
inside. 

The cooling of the reaction mixture which is required 
during the exothermic reaction is achieved by the cold water 
flowing through 1/4" stainless steel cooling coils provided 
inside the reactor. The cooling water, whose flow rate is 
controlled by a Foxboro feed back controller, flows through 
a control valve and a Foxboro D/P cell that is used to 
measure the flowrate of the cooling water. The inlet and 
outlet temperatures of the cooling water are measured by 
thermocouples. Through the above mentioned facility of 


heating and cooling it is possible to maintain the desired 
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temperature within the reactor within +0.5°C. 

A feedback controller is designed to control the 
temperature inside «the reactor... For. this. purpose the 
temperature of the reaction mixture is measured by a 
copper-constantan thermocouple and the resulting millivolt 
signal is converted to a 10-50 ma signal through a Foxboro 
EMF converter. This 10-50 ma signal is fed to a PI Foxboro 
controller which produces the appropriate control signal 
(10-50 ma) to an electro-pneumatic converter. The output of 
the electro-pneumatic converter (3-15 psi) is fed to a 
diaphram control valve to regulate the flowrate of the 
cooling water. The Foxboro controller can accept remote 
setpoint, a facility that can be used for computer control 
of the reactor. The millivolt signal resulting from the 
thermocouple inside the reactor is also sent to the computer 
for on-line data acquisition of temperature. 

A viscometer (Norcross Corporation Model M8B) is 
mounted on the lid of the reactor to measure viscosity 
during the course of polymerization. This viscometer 
Operates son ica dfalhingerpistonweprineiple.® Apypiston is 
periodically raised by a lifting mechanism drawing a sample 
of the liquid to be measured down through the clearance 
between the piston and the cylinder into the space which is 
formed below the piston rod assembly. The piston is then 
allowed to fall by gravity, expelling the sample through the 
Same path as it entered. The time required for the piston 


rod assembly to drop to the bottom of its travel isa 
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measure of viscosity of the liquid and this measure is 
available as an electrical signal in the range of 1-5 volts 
for on-line data aquisition. The measuring element is 
rigidly mounted with the measuring tube in a vertical 
position as indicated by an air bubble level indicator 
mounted on the base of the measuring element. The 
calibration and measurement of viscosity will be dealt in 
detail in chapter three. 

During polymerization the reaction mixture is 
continuously circulated through a densitometer by a 
reciprocating-revolving pump (Figure 2.2). This pump is 
valveless and can handle highly viscous fluids such as 
polymer solutions. According to manufacturer's specification 
it has an accuracy of one percent. In our experiments we 
experienced plugging of the pump at high viscosities. These 
problems were temporarily solved by removing the pump and 
cleaning the pump head with a solvent. The inlet of the pump 
was connected to the stopstock outlet at the bottom of the 
reactori» through ©1/8" flexible.» Teflon: «tubing. TheTefilon 
tubing was found to perform well whereas Tygon tubing 
Started to break at the higher temperature (about 75°C) of 
the reaction mixture. To ensure that density is measured at 
a constant temperature, the reaction mixture is first passed 
through stainless steel coils immersed in a thermostatic 
bath before entering the densitometer. The stainless steel 
coil is of sufficient length (50cm) so that the reaction 


mixture would reach the temperature of the thermostatic bath 
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(60°C), irrespective of the actual temperature of the 
reactor. 

An on-line digital densitometer (Anton Parr, Austria 
Model DMA 450 marketed by Mettler Company, USA) is used to 
measure the density of the reaction mixture. This meter 
measures density by measuring the change in the natural 
frequency of a hollow oscillator when filled with different 
liquids or gases. It has a precision of 0.0001g/ml when the 
range of density is between -0.50:and 1.5 g/mL. Water is 
Circulated around the density cell from a thermostatic bath 
so that density cell and its contents are kept at a constant 
desired temperature. The densitometer is provided with 
stainless steel adaptors so that continuous measurements can 
be made when reaction mixture flows through the density 
cell. The connections between the Teflon tubing and the 
adaptor are made by inserting the tube firmly into the 
sockets of the adaptors. In the event of the plugging of a 
valve or any part of the tube, the connections break open at 
the adaptors thereby avoiding pressure buildup and 
subsequent damage to the density cell. The measured density 
value is available as a BCD output for data acquisition by 
the computer. The measurement and calibration of the 
densitometer will be treated in detail in chapter three. 

From the densitometer the reaction mixture flows 
through a three-way solenoid valve and back to the reactor. 
This valve is used to collect samples of reaction mixture 


automatically for the SEC sampling system. Initially used 
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Buna rubber and Viton valve seats plugged very often whereas 
the Teflon valve seats were found to perform well. 

A size exclusion chromatograph (SEC Waters Associates 
Model 244) is employed to measure molecular weight 
distribution of the polymer in the reaction mixture. It 
separates the polymer molecules according to their size in 
solution by using microporous solids packed inside the 
columns of SEC. The concentrations of the separated polymers 
are detected uSing two detectors (UV absorption, or 
Refractive Index detector). Even though both of the 
detectors are available in the SEC employed in this work, it 
was found that only RI detector responds to the _ polymer 
polymethylmethacrylate. The RI detector has three possible 
outputs i.e. 10mV, 100mV, and 10mA one of which can _ be 
chosen by a selector switch. The RI responses are amplified 
filtered and transmitted? to the computer HOF data 
acquisition. A computer program has been developed for 
automatic data acquisition and computation of molecular 
weight averages from the SEC measurements and will be 


discussed later in chapter three. 


2.3 Automatic Sampling System 

As mentioned earlier an automatic sampling system has 
been developed to collect samples at regular intervals from 
the reactor and prepare it for injection into the SEC. The 
schematic diagram of the system is shown in Figure 2.2. A 50 


mL conical flask is used to collect the sample from the 
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Table 2.1 Sequence Timer Channel Allocation 


ee ee ee 


Channel Output Channel 
Number Number Output 
J NE es ne e ee 

0~7 1 8 
Bor 2 4 
h2=35 3 4 
16 4 1 
17 5 1 
18 6 1 
19 i) 1 
20 8 1 
a1 9 1 
Ze 10 1 
23 11 1 
24 12 1 
25 ie 1 
26 14 1 
27 15 1 
28 16 1 


reactor. The sampling system also includes 1) two three-way 
solenoid valves 2) two metering pumps 3) a magnetic driven 
Stirrer and 4) an automatic sample injection valve. A 
programmable microprocessor based sequence timer (Potter and 
Brumfield Company) that has twenty-nine channels and sixteen 
Output relays is used for switching the pumps, solenoid 
valves and the magnetic stirrer on or off. The sample 
preparation is achieved by appropriate timing of the above 
devices. The essential steps of sampling system in the order 
of execution are: 
= Step 1: collection of a certain quantity of polymer 
mixture in the sample vessel by switching on the 


threeway solenoid valve (S1) for a predetermined period 
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of time. 
= Step 2: dilution of polymer sample with a predetermined 
amount of solvent by operating the solvent pump (P1) 
Over a period of time. 
= Step 3: mixing of the sample with the solvent and 
circulation of the diluted polymer through the SEC 
sample loop by switching on the magnetic stirrer and 
pump(P2). 
~ step 4: automatic injection of sample into SEC by 
Switching on the automatic sample valve. 
a Step 5: emptying contents of sample vessel to waste by 
Switching on the solenoid valve (S2). 
= Step 6: flushing sample vessel with more solvent by 
operating pumps P1 and P2. 
The above mentioned procedure accomplishes automatic 
preparation and injection of the polymer sample into SEC. 
Table 2.1 gives the assignment of channels to the 
output relays. Each channel can be programmed to switch on 
and off once during a cycle. Since eight channels(0-7) are 
‘assigned to the output relay number one, it can be switched 
on and off eight times in a cycle. Similarly from Table 2.1 
it can be observed that output relays numbered two and three 
can be switched on and off four times in a single cycle 
whereas all other output relays can be switched on and off 
only once in a single cycle. 
Even though there are 16 output relays available in the 


sequence timer, only six of them are needed for automatic 
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Table 2.2 Automatic Sampling System Program 


Cycle Time = 15 min 
eat ee ae eS ee 
Start Stop Total Output 
Time Time Time Number Remarks 
min min min 
Pe 
0.0 ON 0.1 4 Polymer sample collected 


using solenoid valve(S$1) 


Oat 0.4 Ot3 1 Dilution of polymer 
sample by pump(P2) 


O73 ees) 12 = Mixing sample using 
magnetic stirrer. 


O23 ‘5 12 A Pumping sample through 
_SEC loop and to waste(P2) 


0.5 0.8 O23 6 Sample injection by 
automatic sample valve 


OX6 Ono O23 2 Discharging the polymer 
to waste(S2) 


SS. Ny 0.3 1 Pumping more solvent 
for cleaning(P1) 


ye 2 bas o.3 2 Discharging the solvent 
to waste(S2) 


Sampling. Table 2.2 gives the program used in this work with 
the start time and stop time of the various devices used in 
the automatic sampling system. Since it takes 10-12 minutes 
of time to analyze a polymer sample using SEC, a cycle time 
of 15 minutes is chosen for this program. Every fifteen 
minutes the program given in the Table 2.2 is executed by 


the sequence timer, starting a SEC analysis of the polymer 


sample. 
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2.4 Computer System 

A Hewlet-Packard (HP/1000) minicomputer is used for 
data acquisition and control of the experimental reactor. 
This computer is maintained by Data Acquisition, Control, 
and Simulation Centre (DACS Centre) of the Department of 
Chemical Engineering. 

The HP/1000 minicomputer employs a PDP 11/03 computer 
as an I/O subsystem to which various measurements from _ the 
reactor are interfaced. The PDP 11/03 computer performs the 
Following functions; i) A/D conversion of analog signals and 
transmission of the digital signal to the main computer; ii) 
processing of thermocouple signals after providing cold 
junction compensation; iii) reading status of digital 
Switches; and, iv) routing analog and digital outputs to 
appropriate devices. The system supports a plotter that can 
be used to plot data collected during an experiment. A hard 
copy terminal near the reactor site is used for entering 
commands as well as reporting experimental data. The data 
acquired by the computer is also stored in disc files for 
Subsequent analysis. 

The programs are written in Fortran IV, the language 
Supported by the computer system. Real time data acquisition 
is done through CALL-s to standard ISA (Instrument society 
of America) Subroutines. 

In addition another HP/1000 minicomputer which is 
dedicated to the task of gas chromatography has been used to 


do the GC analysis of the polymerization mixture. Details of 
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the GC analysis will be given in chapter three. 


2.5 Experimental Procedure 

A free radical-initiated solution polymerization of 
methylmethacrylate in a batch reactor has been selected to 
study the problems associated with on-line implementation of 
optimal control policies. Toluene is chosen as a_ solvent 
beeause. of its higher boiling=point (120°C) and its low rate 
constant for chain transfer reaction. Because the 
boiling-point is high, the toluene does not evaporate 
Significantly at the temperature (60-80°C) of 
polymerization. The molecular weight distribution of the 
polymer formed during the course of polymerization is not 
affected by the presence of toluene due to the low value of 
its chain transfer rate constant (Brandrup and 
Immergut(1975)). 

Commercial monomer methylmethacrylate (MMA) (supplied 
by Aldrich Chemical Company) contains 65 ppm hydroquinone 
monomethyl ether as inhibitor which is removed by vacuum 
distillation at 32-35°C. Both monomer and solvent are dried 
with anhydrous calcium chloride, since water also inhibits 
the reaction. Both monomer and solvent are purged with 
nitrogen to remove dissolved oxygen which inhibits 
polymerization. The initiator, benzoyl peroxide, supplied by 
Aldrich Chem. is dissolved in chloroform and recrystallized 
in methanol. The reactor is purged with nitrogen during 


polymerization to prevent oxygen from being dissolved into 
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the reaction mixture. 

Measured quantities of monomer and solvent are charged 
into the reactor and the reactor contents are heated to. the 
required operating temperature, while the contents are mixed 
thoroughly by the stirrer. After the contents of the reactor 
have reached the reaction temperature, the initiator which 
has already been dissolved in toluene, is added to the 
reactor. Immediately afterwards the reactor is placed under 
computer monitoring for the purpose of on-line measurement 
of variables of the polymerization process. Off-line 
measurements are also carried out on the samples that are 
collected periodically from the reactor. These samples were 
also stored in a freezer so that they could be used for 
later off-line analysis if necessary. 

A computer program has been developed for the purpose 
of on-line monitoring of various variables of the process. 
Every five seconds the program samples’ the following 
measurements during the course of polymerization:reactor 
temperature, inlet and outlet temperatures of hot water 
flowing through the jacket, inlet and outlet temperatures of 
the cold water flowing through the cooling coils, cold water 
flowrate, density and viscosity of the reaction mixture and 
finally the torque on the stirrer. The above measurements 
are corrupted by the noise of the measuring devices and also 
by the transmission noise and they should be filtered to 
Smooth the data. A moving average filter that smoothes the 


current data by taking an average with certain number 
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(typically 10) of past data is used in this program. The 
filtered data are reported at specified intervals of time 
(one or two minutes) at a hard copy terminal installed near 
the reactor. Also, these data are stored in disc files in 
the computer for later analysis. 

In addition the program reads a digit switch every 
second to check whether a polymer sample has been injected 
into the SEC by the automatic sampling system. If a sample 
has been injected, as inferred from the ON status of the 
digit switch, the program automatically starts the SEC data 
collection and reports these data at specified time 
intervals on the hard copy terminal. After the polymer peak 
has eluted the SEC data collection is automatically 
terminated by the program. The details of the SEC program 


will be given later in chapter three. 


2.6 Conclusions 

The experimental facilities that are required for 
on-line measurement and computer control studies on a batch 
polymerization reactor have been described in this chapter. 
It can be appreciated that the experimental reactor system 
requires extensive instrumentation to measure on-line 
properties like density, viscosity, torque, temperature, 
molecular weight of the polymer etc. A _ process control 
computer was employed for on-line data acquisition and 
implementation of optimal control policies son= the reactor 


system. A description of the software that was developed for 
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data acquisition of all on-line measurements has also been 


given. 
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3. MEASUREMENT STUDIES 


3.1 Introduction 

As pointed out by Hoogendoorn et al. (1980) the 
measuring instruments are often the weakest link in 
polymerization reactor control circuits. Without some form 
of on-line polymer characterization, the monitoring and 
feedback control of polymer quality is not possible. Hence 
there is great incentive to develop on-line sensors to 
measure the polymerization process variables such as 
conversion and molecular weight distribution. There has been 
considerable advances in the development of on-line sensors 
recently(MacGregor et. at(1983)). In this chapter we 
describe some on-line techniques used in this work to 
measure density, viscosity, torque and molecular weight 


distribution in a solution polymerization process. 


3.2 Previous Work 

In the past, on-line measurements in a polymerization 
reactor were mainly restricted to temperature and pressure. 
Even temperature measurements can sometimes present problems 
in polymerization reactors due to fouling of the measuring 
elements(MacGregor et al.(1983)). Traditionally conversion 
has been measured by gravimetric methods which are time 
consuming and necessarily measured off-line. The molecular 
weight distribution of the polymer is usually measured 


off-line using a size exclusion chromatograph. There have 
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been only a few attempts recently to develop on-line sensors 
to measure the above mentioned process variables. 

Jo(1975) used an on-line refractometer supplied by 
Electron Machine Corporation (ECM) to monitor conversion in 
solution polymerization of vinylacetate. He modified the 
refractometer to increase turbulence and prevent polymer 
coating on the prism and reported trouble free operation of 
the instrument as well as accurate reproducible results for 
conversion. Schmidt and Ray(1981) have also used an on-line 
refractometer from ECM to monitor conversion in solution 
polymerization of methylmethacrylate and reported accurate 
results for conversion. 

Abbey(1981) used an on-line densitometer (Anton Paar, 
Austria) for monitoring conversion in the cases of solution 
and emulsion polymerizations of methylmethacrylate and 
reported the following difficulties: 

4 pump failure from either monomer attack or polymer scale 
= monomer phase separation in density cell 

3 lag time for rapid polymerization. 

Schrock and Ray(1981) also used a DMA series densitometer in 
emulsion polymerization of methylmethacrylate and observed 
very good agreement between conversion from gravimetry and 
densimetry. They reported that at least a resolution of 0.5 
per cent conversion is possible in the presence of sampling 
errors. Recently Pollock(1983) used the same type of 


densitometer to monitor conversion in emulsion 


polymerization of vinylacetate. 
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Though gas chromatograph is a sensitive and reliable 
instrument, it has not yet been extensively used in 
polymerization studies (Collins et al. (1973)). Guyot et al. 
(1981) used a gas chromatography with an automatic sampling 
System to monitor conversion in emulsion polymerization of 
vinylacetate. In some cases the gases contained in the head 
Space of the reactor are analyzed and the composition and 
conversion of monomer in the reacting liquid are inferred 
through vapor-liquid equilibrium relationships. Recently 
Garcia-Rubio (1982) investigated the styrene acrylonitrile 
copolymerization kinetics using gas chromatography. Although 
some difficulties were experienced with the monomer 
extraction process at high conversion, in general very good 
descriptions of the copolymerization reaction were obtained 
in terms of the residual monomer composition. 

Even though size exclusion chromatograph (SEC) is 
widely used to measure molecular weight distribution (MWD), 
only a few applications for on-line measurements are 
reported. Gregges et al. (1971) have modified a SEC for 
automatic sample injection and on-line data acquisition by a 
computer. Recently Meira et al. (1981) employed an SEC 
on-line for measuring MWD in the solution polymerization of 
methylmethacrylate and reported good reproducibility of SEC 
results. However they did not include an automatic sampling 
system which can collect samples from the reactor and inject 
into SEC after diluting with a carrier solvent. Roof et al. 


(1980) have field tested an on-line SEC in the case of a 
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solution polymerization process over a period of three years 
and reported that their experience has proven that on-line 
polymer analysis by SEC is practical. They also developed an 
automatic sampling system to collect samples from the 
reactor and inject into the SEC. 

Jo and Bankoff (1976) measured viscosity on-line in 
solution polymerization of vinylacetate and used the 
information to estimate the weight average molecular weight 
of the polymer product. 

From the above it is clear that on-line sensors are 
still under development in polymerization processes. In the 
next section we will deal with on-line as well as off-line 
methods to measure conversion, molecular weight and 


viscosity in solution polymerization of methylmethacrylate. 


3.3 Conversion Measurements 

There are a number of methods available to meaSure the 
extent of conversion of monomer to high polymer during the 
course of polymerization. Probably the most direct method of 
determining conversion of monomer to polymer is to stop the 
polymerization, and isolate and weigh the resulting polymer. 
The chemical determination of the remaining functional 
groups at some stage of the polymerization is often 
practiced, particularly for step polymerization. This method 
not only leads to the extent of reaction but, through the 
use of Carothers equation relating extent of reaction and 


degree of polymerization, to the number average molecular. 
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weight in the case of linear polycondensation. Double bonds 
can also be determined chemically, so that this class of 
methods is applicable to chain polymerization as well. 
Examples of use of this method are sreported» -in» ithe 
literature (Patton(1962), Sorenson(1968), Kolthoff£(1947), 
Lewis(1945)). 

There is a linear relationship between the degree of 
conversion in a polymerization system and its volume, 
provided only that the monomer and polymer mix without 
change in volume. For following the volume change during 
polymerization, the reaction is carried out in a 
Gilatometer, a vessel equipped with a capillary tube in 
which the liquid level can be measured precisely using a 
cathetometer. The decrease in volume on polymerization is 
large being 21% for methylmethacrylate and 27% for 
vinylacetate. In practice a few hundredths per cent 
polymerization can be detected (Billmeyer, JIPSeGte7 15 
Example of the use of dilatometry to study the course of 
polymerization reaction has been reported by many authors 
(Paoletti(1964), Bell (1961), Niezette(1971)). 

The determination of extent of reaction by means of 
refractive index is based upon the refractive index of the 
liquids, their chemical composition and their molecular 
Structure. Since the structure of the polymer differs from 
the structure of monomer because of rearrangement of 
chemical bonds, the refractive index of the polymerization 


mixture ig different from that of pure monomer. Since 
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refractive index is also dependent on temperature, it is 
necessary to keep the temperature of the refractometer 
constant by precise temperature control. Generally, the 
dependence of refractive index on conversion is determined 
by calibration using polymerization mixtures whose 
conversions are known. Conversion can also be calculated by 
using the Lorentz-Lorentz equation (Collins et al. (1973)). 
Although the refractive index method requires removal of a 
Single drop of sample of polymerization reaction mixture per 
determination, it is a Simple and rapid means of following 
the polymerization. Naylor(1953) and Pulley(1968) used 
sealed cylindrical tubes as reaction vessels and measured 
refractive index in situ to follow the polymerization 
process. 

On-line refractometers are available and their use _ to 
measure conversion in polymerization processes has been 
mentioned earlier in this chapter. 

In addition to the above mentioned methods, the extent 
of reaction during polymerization can also be determined by 
densimetry, gas chromatography and infrared spectroscopy. 
Though the number of methods available to measure conversion 
in polymerization reactions are many, only a few of them can 
be used for on-line measurement. Three methods, namely 
gravimetry, gas chromatography and Gensimetry that have been 
used for on-line and off-line determination of conversion in 
solution polymerization of methylmethacrylate (studied in 


this work) are described in the following section. 
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3.3.1 Gravimetry 

This is the most direct method of determining 
conversion of monomer to polymer. The Sample of the 
polymerization mixture is withdrawn from the reactor ina 
vial and further polymerization of the collected sample is 
stopped rapidly by the addition of a fast-acting shortstop 
(in this case a few drops of solution of hydroquinone in 
distilled water). These vials are stored in a freezer so 
that gravimetric as well as gas chromatographic analysis can 
be done at any future time. 

In this study the procedure reported by Jo(1975) has 
been used to isolate polymer from the sample. Samples of the 
reaction mixture are placed in preweighed aluminum dishes 
and weighed accurately using a Mettler electronic balance 
which has a resolution of 0.00001 g. The samples are _ then 
dried in a vaccum oven and finally heated to 60-80°C to 
remove the traces of monomer and solvent that could have 
been entrapped inside the polymer. The isolated polymer in 
the dish is cooled and weighed to a constant weight. Since 
the initiator added to the reactants in the beginning of the 
reaction remains with the polymer after vacuum drying, a 
correction is made in the conversion calculation to account 
for the weight of initiator in the isolated polymer. Because 
only a few drops of the solution of inhibitor are added to 


the sample, a correction for the weight of inhibitor is not 


done in our calculations. 
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This method of determining conversion is necessarily 
time-consuming, involving several hours of analysis, and 
moreover it becomes increasingly difficult to remove the 
traces of monomer and solvent from the polymer at high 
conversions where the viscosity of the sample mixture is 
also high. It was suggested to precipitate the polymer in 
the high conversion samples in methanol or heptane and then 
to isolate the precipitated polymer by vacuum drying to 
avoid entrapped monomer or solvent in the polymer. However, 
it was found that both of the gravimetric procedures, namely 
direct vacuum drying and precipitation, gave almost similar 
results suggesting that conversion results obtained from 
gravimetry are prone to error especially at higher 
conversion levels. The residual monomer and solvent trapped 
inside the polymer sample can be removed by heating the 
vacuum oven to higher temperatures. However heating of the 
samples could also lead to further polymerization of the 
residual monomer to polymer and hence better results of 
conversion need not be obtained. 

The presence of the residual monomer and solvent in the 
dried polymer has been confirmed by dissolving the dried 
polymer and injecting into an SEC. An SEC peak corresponding 
to the monomer and solvent elution volume has been observed 


thus confirming the presence of monomer and solvent . 
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3.3.2 Gas Chromatography 
Though gas chromatography is perhaps better suited than 
any other method to follow the course of polymerization, it 
has not yet been extensively used. As a rule the presence of 
polymer in the reaction mixture does not hinder the use of 
GC. Since the methods for determining volatile components in 
polymer systems have been developed in great detail, GC 
methods can be used directly for determining the reaction 
kinetics from the changes in the concentration of the 
monomers consumed. Berezkin et al.(1977) compared the widely 
used dilatometry with GC methods to determine kinetic 
parameters in polymerization reactions and concluded that GC 
is the preferred method in most instances. 
The following is a list of advantages of using GC 
method in the study of polymerization reactions (Berezkin et 
al .(1977) 0 
= The sample volume required for GC analysis is extremely 
small in the order of microlitres. 

es Since GC analysis directly gives the mass fraction of 
the volatile components in the polymerization mixture, 
the results do not depend upon the temperature of the 
reaction. We need not worry about temperature effects as 
is the case with some other methods. 

ss By using an internal standard in the polymerization 
reaction mixture, one need not worry about the exact 
volume of the sample injected. Yet accurate results can 


be obtained about the compositions of the components in. 
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the sample. 

GC analysis can also be extended to study any side 
reactions that may take place along with polymerization 
reaction, 

The sampling procedure can be automated so that GC 
methods can be used for on-line monitoring of 
polymerization reactions. 

GC methods can be easily extended to _ the study of 
copolymerization reactions with several monomers’ for 
which other methods may not be applicable or require 
extensive calibration. 


following are some of the problems that may be 


encountered when using GC methods in polymerization studies. 


The polymerization samples with higher viscosities 
corresponding to higher conversions can not be drawn 
into GC syringe easily for injection. However these 
samples could be diluted with a_ solvent which will 
facilitate easy injection into the GC. 

Since the polymer does not evaporate in the injection 
port, the polymer particles accumulate and block the 
injection port in course of time. The injection port can 
be removed and cleaned with a solvent. 

The injection port temperature should be such that there 
is no polymerization of monomer arid (*ealsoc tino 
depolymerization of polymer to monomer each of which can 


affect the results obtained from the GC. 
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In this work a HP gas chromatograph (Model 5710A) with 
a thermal conductivity detector is used. Argon is used as a 
Carrier gas. The following operating conditions are employed 
to analyze the sample: 


Column temperature 10226 


Detector temperature 1SOe C 


Injection port temperature : 150 °C 

Carrier gas flow rate ¢ 15 mL/min 

Sensitivity of the detector : 4 

TC current te 150" ma 

Carrier gas : Argon 

For analyzing the samples of the reaction mixture, the 
hardware and software supplied by Hewlett-Packard Company 
have been used. The hardware consists of a A/D converter and 
a HP/1000 minicomputer which is dedicated to the task of gas 
chromatography for a number of GC-s used in the Department 
of Chemical Engineering of University of Alberta. The 
software is a GC program that can automatically analyze the 
GC-readings, find the peaks and report compositions. The 
analysis is started by injecting the sample into the GC 
injection port and pressing a switch located on the A/D 
converter. The results are reported on a teletype terminal 
installed near the GC. 

The solvent toluene does not take part in the reaction 
and so its weight fraction in the reaction mixture remains 
constant throughout the course of polymerization. Hence the 


solvent serves aS an internal standard for the GC analysis 
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of the monomer methylmethacrylate. During the course of 
polymerization the composition of monomer decreases with 
respect to the solvent and conversion of monomer to polymer 
can be easily calculated from the knowledge of initial 
compositions of monomer and solvent and the compositions 
obtained from the GC analysis. Since toluene serves as an 
internal standard, exact volume of injection of the GC 
sample need not be known. 

For the purpose of comparison, the conversion results 
obtained from GC and gravimetry are presented in Figure 3.1. 
In Figure 3.1 T, ©®,, and I 9 denote temperature, volume 
fraction of solvent and initial initiator concentration, 
respectively. It can be observed that conversion results 
from both of the methods agree very well (within 3%) upto 40 
per cent conversion whereas after about 40 percent 
conversion the gravimetric results give higher values 
compared to the GC results. These higher conversion values 
from gravimetry especially at higher conversion range are 
attributed to two possible sources of error which have been 
already discussed in the section on gravimetry namely 
entrapped monomer and solvent in polymer and possible 
polymerization of monomer to polymer due _ to heating of 
samples. These errors are absent at lower conversion levels 
and hence better agreement of results are observed between 
the two methods. 

In the GC method of analysis possible sources of errors 


are also due to polymerization of monomer to polymer inside 
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Fig 3.1 Comparison of Conversion Results fromGravimetry and GC 
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GC columns and the possible depolymerization of polymer to 
monomer in the injection port. The first source of error 
leads to conversion values from GC higher than true 
conversion. However, the vapor phase thermal polymerization 
of methylmethacrylate is negligible thus eliminating the 
error due to this source. The second source of error, i.e. 
depolymerization of polymer would lead to the GC conversion 
results lower than the actual conversions. But the actual 
injection port temperature of 150°C is much lower than the 
temperature at which depolymerization would occur. According 
to Odian(1981), depolymerization would become pronounced at 
temperatures higher than 200°C. It can be concluded from the 
above arguments that the GC method of analysis gives 
conversion results that are relatively accurate. 

Even though, normally gravimetric methods are used to 
measure conversion in polymerization processes, one should 
be careful to avoid the sources of errors described here. In 
particular a quantitative check using size exclusion 
chromatography would be helpful to decide whether the 
amounts of monomer and solvent trapped inside the polymer 
are considerable. Only then would the validity of 
gravimetric results be established. For all subsequent 
conversion values needed for the study of polymerization 
reactions, the results obtained from the GC analysis have 
been used. 

In retrospect, a better way of comparing conversion 


measurements from GC and gravimetric analysis would have. 
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been to prepare solutions of polymer, monomer and solvent by 
dissolving known weights of PMMA in a mixture of monomer and 
solvent and analysing the resulting mixture by the above 
methods. Since the weight fraction of the polymer in the 
mixture is known, we could have demonstrated the accuacy of 


one of the two or both methods. 


3.3.3 Densimetry 

Since the density of the polymer is higher than the 
density of the monomer, the density of the polymerization 
mixture increases during the course of the reaction, thus 
making density measurement an attractive method to follow 
the course of the polymerization reaction. Dilatometry which 
exploits this phenomenon of volume change (and hence density 
change) during the course of polymerization has been widely 
used earlier mainly to study kinetics of polymerization 
reactions. This section will deal with commercially 
available densitometers that can be used for continuous 
measurement of density and on-line control of conversion in 
polymerization reactions. 

However, it should be mentioned that an on-line 
refractometer could also have been used to monitor 
conversion. In fact, we have recently purchased an on-line 
refractometer from the Electron Machine Corporation to 
evaluate its use in monitoring conversion. 

Several companies supply densitometers which are 


Suitable for automatic, continuous operation with sufficient 
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precision for calculation of polymerization conversion. 
These break down into three classes based on the mode of 
operation: 

= y-ray absorption 

= oscillatory frequency of a sample filled tube 

= mass measurement at fixed volume. 

Only one of these, an oscillator-based system 
distributed by Anton Paar has models with dead volume small 
enough for laboratory scale on-line determination of 
density. 

The measuring principle of the instrument is based on 
the change of the natural frequency of a hollow oscillator 
when filled with different liquids or gases. The mass and 
thus the density of liquid or gas changes this natural 
frequency due to a gross mass change of the oscillator 
caused by introduction of liquids or gases. The oscillator 
consists of a hollow elastic glass tube which is 
electronically excited in an undamped harmonic fashion. With 
the simple relationship between the density of the sample 
and sthe tnatural sfrequency-yof sthe filledqaoserllatorpoaitatis 
possible to use the method for the determination of the 
density of the samples which are injected into or flowing 
through the oscillator. 

For calculation of density we may consider an 
equivalent system represented by a hollow body of mass, M, 
which is suspended on a spring with an elasticity constant, 


C, its volume V filled with a sample of density p. 
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The natural frequency f£ of this system will be given by 


1 c 


ih= —> (33:18 
27 M+ pV 


therefore, the period T is given as 


T= 12a 2 (MF pV) /C Gane) 


Taking the square of this expression and introducing the 


constants 


A=421?2V/C B=4n?M/C C303) 


we obtain the equation for the density,p, of the sample from 


the period T as 


p=(T? =°B)/A (3.4) 
Since the constants A and B contain the volume, = spring 
constant and mass, they may be regarded as apparatus 
constants which may be determined from two calibration 
measurements of samples of known density such as air and 
water. 

This suggested procedure by the supplier, Anton Parr, 
was used to calibrate the densitometer. The instrument is 
connected to an ultrathermostat and brought to the desired 


temperature by an adjustment on the thermostat temperature 
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control. The display selector located at the top of the 
instrument is set to "T" to show the period of oscillation. 
Then approximately 0.7 ml of distilled water is injected 
through the lower opening of the oscillator cell using a 
plastic tipped hypodermic syringe. The sample must be 
homogeneous and free of even the smallest bubbles. The 
syringe is left in this position at the lower opening of the 
oscillator tube until the display shows a constant value 
indicating temperature equibrium. 

Next the sample in the cell is withdrawn and the 
distilled water remaining inside the cell is evaporated by 
pumping air through the cell. After the cell has become dry, 
the air inside the cell is allowed to reach temperature 
equilibrium and the period of oscillation due to air is 
noted. From the values of periods of oscillations of 
distilled water and air, and from knowledge of densities of 
distilled water and air, the values of the constants A and B 
can be calculated using Equation(3.3). This procedure is 
repeated for the calibration at different temperatures. 

Table 3.1 gives the results of the calibration for the 
DMA 450 series densitometer employed in this work. Once the 
instrument has been calibrated, normally it does not require 
recalibration until the density cell is replaced. This, in 
fact, has been confirmed by the author. When the instrument 
was recalibrated after several months of use, no significant 


changes in the calibration constants were observed. 
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After calibration of the densitometer, density of the 
polymerization mixture can be measured easily. However, it 
1s important that these density measurements should be 
converted to conversion measurements by developing 
appropriate correlations or equations. For this purpose two 
methods can be used. 

The first one involves the development of an empirical 
correlation between density and conversion and use of this 
correlation to estimate conversion from density measurement. 
Figure 3.2 shows the variation of density with conversion 
for various polymerization conditions. For purpose of 
development of correlations, the density measurement is done 
off-line at a temperatute of 40°C which is significantly 
lower than the reaction temperature(60-80°C) in order to 
avoid further polymerization of the sample in the density 
cell. At 40°C polymerization of the mixture essentially 
stops and hence the calculated conversion values correspond 
to conversions at the times of collection of samples from 
the reactor. The conversion values used for this calibration 
are obtained from the GC meaSurements as the GC method is 
considered to be more accurate compared to gravimetry. 

The density change with conversion during the 
polymerization process is shown in Figure 3.2 for various 
conditions. On closer examination it can be observed that 
two different straight lines can be fitted to the data: one 
for the initiator concentration 0.05 mol/l, the other for 


the concentration .10 mol/l. In fact the two solid lines 
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shown in the Figure 3.2 are the least square fit of the data 
for the two initiator levels. However, it is noted that if 
the volume fraction of the solvent is changed, a different 
straight line has to be fitted for the experimental points. 
It is obvious that this procedure of fitting a straight line 
for each experimental condition is not convenient. Hence a 
second method of relating density and conversion has been 
developed as described below. 

Abbey (1981) assumed ideal mixing of polymer, monomer 
and solvent or water phase and developed mathematical 
equations relating the density of the polymerization mixture 
and conversion in the cases of solution and also emulsion 
polymerization of methylmethacrylate. Schmidt and Ray (1981) 
also used a mathematical relationship to get conversion in 
the case of emulsion polymerization of methylimethacrylate. A 
mathematical equation relating the density to the conversion 
will be developed assuming ideal mixing of polymer, monomer, 
solvent and also the initiator. Since the initiator also 
affects the density of the polymerization mixture, as shown 
in Figure 3.2, the effect of initiator on the density must 
be taken into account. 

Let W be the total weight of the initial polymerization 


mixture.The total weight is given by the Equation(3.5). 


W=WmtW,+tW; (325) 
where W,, W; and W;, are the initial weights of monomer, 


solvent and initiator, respectively. The total weight of the 
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polymerization as well as the weights of solvent and 
initiator remain constant throughout the course of 
polymerization. Since monomer is converted to polymer, the 
weight of the monomer decreases while the weight of the 
polymer increases by an equal amount. The density of the 
polymerization mixture increases due to higher density of 
the polymer compared to the monomer. Let xX be the weight 
fraction of the initial monomer converted to polymer. At a 


conversion X, the following mass balances will hold good. 


the weight of polymer = W,, X. 


the weight of monomer = Wy (1-X) 


volume of polymer = Wmx/Py 
volume of monomer = Wa €145ib/ Ox 
volume of solvent =iWe/p: 
volume of initiator = We/p; 


The total volume V of the polymerization mixture can be 
obtained by adding the volumes of monomer, polymer, solvent 
and initiator assuming ideal mixing of the above components. 


Total volume V is given by: 


WmS Wy Gu &) W, Wi 
V= i (3.6) 


Py Pm Ds pj 


WHETEI Py Pmplpeumand «9;  ,ere the densities of polymer, 
monomer, solvent and initiator respectively.The density of 


the mixture will be given by 
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Poin = W/V (3.7) 
Substituting for W and V from Equations(3.5) and (3.6) 
respectively in #quation(3.7), an expression for monomer 


conversion X is obtained as: 


Wn(1-1/Pm)+W,(1-1/p,)+W; (1-1/0; ) 
Poe 1 ee RRA Le aa (3.8) 
GanreWn/ Dy = Wil Due) 

In the RHS of Equation(3.8) all the quantities except 
density of the mixture are known at the beginning of the 
reaction. Thus by measuring density, conversion X can be 
calculated. However it should be noted that density also 
changes with temperature. To eliminate the effect of 
temperature, the density of the polymerization mixture is 
measured at a constant temperature. 

Equation (3.8) has been used to calculate conversion 
from the measured values of density at a temperature of 40°C 
for various conditions of polymerization. These results are 
compared to the conversion results obtained from GC 
measurements in Figure 3.3. It can be easily observed that 
the conversion results from both of these methods agree very 
well, confirming our earlier deduction that the gravimetric 
results are prone to error at higher conversion level. This 
equation can also be used in Situations where the density 
measurements are available under changing temperatures of 
the reaction mixture if the densities of monomer and solvent 


are known as functions of temperature. Hence its use should 
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Fig. 3.3 Comparison of Conversion Results from Densimetry and GC 
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be superior to that of the empirical relationship shown in 
Baqurencs .21 

The densitometer is capable of measuring density to a 
resolution of 0.0001 g/mL. The resolution with respect to 
conversion in the case of solution polymerization of 
methylmethacrylate can be calculated by knowing the initial 
density and the final density at 100 per cent conversion. In 
the case of bulk polymerization of methylmethacrylate a 
resolution of 0.04 per cent conversion can be expected 
whereaS in the case of solution polymerization with 60 per 
cent solvent (toluene) a resolution of 0.20 per cent 
conversion can only be expected. Schrock and Ray (1981) have 
reported that in the presence of meaSurement errors, they 
eouldargetsta leconversion dresdolbutiron cof i005 per scent «in 
emulsion polymerization of methylmethacrylate. 

The above discussion indicates that the densitometer 
has sufficient resolution to measure conversion in 
polymerization reactions and can also be used for on-line 


monitoring of conversion. 


3.4 Size Exclusion Chromatography (SEC) 

SEC has become one of the most important instruments in 
polymer engineering because of its ability to measure 
molecular weight distribution (MWD). SEC separates polymer 
molecules with molecular weights from 1000 to as high as 10 
million according to their size in solution. The separation 


is accomplished by passing a solution of molecular mixture 
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of interest through a column containing particles of a 
microporous solid which has solvent in the interstices and 
in its pores. Since any pore accessible to larger molecule 
is also accessible to a smaller molecule, larger molecules 
are excluded from more pores than are the smaller ones and 
come out of the column first. This separation is illustrated 
schematically in Figure 3.4. 

Moore (1964) of Dow Chemical Company first disclosed 
the use of cross-linked polystyrene gels for separating 
synthetic polymers soluble in organic solvents and with that 
event conventional size exclusion chromatography (SEC) was 
born and came rapidly into extensive use to ee MWD 
information for synthetic polymers. Modern high-performance 
SEC has resulted from the development of small, more rigid 
porous particles which can withstand relatively high 
pressures (2000-3000 psi) and provide high performance in 
terms of resolving polymer species. A number of very good 
references are available that deal with SEC separation 
techniques in polymerization studies (Yau et al.(1979), 
Boni(1976), Hamielec (1976)). 

A large number of equipment manufacturers market 
systems and components for performing SEC separation. The 
information concerning the vendors and the equipment that 
can be used for high-performance liquid chromatography are 
listed by McNair and Chandler(1976)| and Snyder and 


Kirkland(1979). 
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An SEC Supplied by Waters Associates (Model 244) was 
used in this work to measure the molecular weight 
distribution of the polymer polymethylmethacrylate. The 
Waters Associates model consists of the following 
components: 
= a high pressure solvent delivery pump 
= a UV absorbance detector 
= a refractive index (RI) detector 


x a sample injection system 


u-Styragel columns 
= a syphon at the outlet of the columns. 
Tetrahydrofuran is used as the carrier solvent and also as 
the solvent to dissolve the polymer. This SEC has been 
modified by replacing the sample injection system with an 
automatic sample valve and by incorporating an automatic 
sampling system which collects the samples of the reaction 
mixture from the reactor and injects into the SEC after 
diluting with the carrier solvent. The automatic sampling 
system which facilates the on-line measurement of MWD has 
already been described in detail in chapter two. 

Even though the SEC is provided with a UV as well as a 
RI detector, only the RI detector was used in this work 
because of its response (being a universal detector) to the 
polymer methylmethacrylate. The RI detectors are universal 
in the sense that they respond to a wide variety of 
chemicals. The UV detector does not respond to the polymer 


PMMA. The RI detector responses (0-10mv) are amplified, 
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filtered and sent to the HP/1000 computer for on-line data 


acquisition. 


3.4.1 SEC Program 

A computer program has been developed in order to 
collect the SEC data automatically and process them to 
calculate the number average and weight average molecular 
weights. The SEC program accomplishes the following 
functions: 
baselnvt ration tof ‘SECs datavacquiksi tion 
2. Identification of start of SEC peak 
3. Identification of end of SEC peak 
4, Termination of SEC data collection 
5. Processing of SEC data for molecular weight averages 
Initiation of SEC data acquisition iS accomplished by 
reading a digit switch associated with the switching of the 
automatic sample valve. This digit switch which is located 
in the programmable sequence timer turns on at the instant 
of injection of polymer sample into the SEC. The program 
reads the on-status of the switch and starts the process of 
SEC data: acquisition. Since the actual SEC peak which 
contains the essential data elutes a few minutes after the 
injection, a dead time can be specified in the program 
during which SEC data will not be collected. However time 
count is started from the moment of injection. Few minutes 
after the injection the digital switch is turned off so that 


at the time of next injection, it can again be turned on 
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starting another SEC analysis. 


Logic to Identify the Start of the Peak 

The SEC data that are necessary for the calculations of 
molecular weight averages lie between the start point and 
the end point of the SEC peak. Hence it is important to 
devise proper logic to correctly identify the start point of 
the SEC peak. For this purpose initially 20 SEC readings are 
collected at the appropriate sample interval which in this 
case is three seconds. These readings are stored in an array 
and starting from the first reading, each reading in the 
array is analysed to check five consesecutive readings 
showing an increasing trend. The first reading which 
Satisfies the above criterion is considered to be the start 
point of the SEC chromatogram. This criterion works well 
upon proper differentiation between baseline drift and SEC 
peak. For this purpose each successive SEC reading after the 
Start point is required to be higher than the previous value 
by a certain value (10mV), which is usually higher than the 
baseline sdriftisvalue(s5mV). 81 £ call «the ‘pointsomfaibeto 
Satisfy the criterion for the start point of the SEC peak, 
it implies that the SEC peak has not yet started eluting. 
Since the SEC data outside SEC peak are not needed, the 
program rejects the initial five data and proceeds to 
collect five more data points so that there are again twenty 
points to repeat the same logic to find the start point of 


the peak. After othe) start point of ~SEC peak ‘has v«been 
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identified a logical variable LSTAT which is initialized to 
-FALSE. 1S set to .TRUE. and this variable helps the program 
to recognize the occurrence of the peak and the future 


readings are not analyzed for the start point of the peak. 


Logic to Identify the End Of The Peak 

After the start point has been identified, the next 
important point is the end point of the chromatogram. For 
this purpose, the program first identifies the peak point 
corresponding to the maximum value of the chromatogram 
readings. For identifying the peak point of the 
chromatogram, the SEC data are arranged in an array and 
Since the peak point can occur only after the start point, 
each reading in the array after the start point is analyzed 
to check whether it satisfies the criterion for the peak 
point. A peak is considered to have been identified if five 
consecutive readings show a decreasing trend. The first 
reading after the start point which satisfies the above 
criterion is the peak point of the chromatogram. If a peak 
point is not found, the program stores a point from where 
the search for the peak point should start next time and 
proceeds to collect five more SEC data. The same criterion 
for the peak point is repeated but this time starting with 
the stored value of the search point, thus reducing the 
computer search time. After a peak point has been identified 
a logical variable LPEAK which is initialized to .FALSE. is 


set to .TRUE. and the program proceeds to identify the end 
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point of SEC chromatogram. 

To identify the end point of the SEC chromatogram the 
SEC data are placed in an array and each reading after the 
peak point is analyzed to check whether it satisfies the 
logic for the end point. A_ point is considered to have 
Satisfied the criterion for end point if five consecutive 
readings have SEC readings that lie within a bound. The 
first reading after the peak point that satisfies the above 
criterion is considered to be the end point of the 
chromatogram. If all the readings fail to satisfy the 
criterion, a search point is stored and the program collects 
five more SEC data to repeat the criterion. After the end 
point has been identified a logical variable LSTOP which is 
initialized to .FALSE. is set to .TRUE. and the program 
stops data collection. It is possible that sometimes a 
solvent peak or initiator or monomer peak may start eluting 
before the end of polymer peak settles to a constant value 
in which case the SEC chromatogram ends in a valley. The 
identification of end point of SEC chromatogram 
corresponding to an end in a valley is also handled by the 
program. However by selecting SEC columns which can achieve 
good resolution of polymer peak, an end in valley can _ be 
eaSily avoided. 

After identification of the end point of the SEC peak, 
the program stops the data collection and reinitializes the 
above mentioned logical variables to .FALSE. so that another 


search for the start, peak and end points can be made in the 
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next polymer sample injection. Then the program proceeds to 
calculate the molecular weight averages. 

The SEC data that had been collected can also be 
plotted using an HP plotter. Figure 3.5 gives a typical 


chromatogram obtained in one of the SEC analyses. 


3.4.2 Molecular Weight Calculations 

The SEC data are processed to calculate number average 
and weight average molecular weights taking into account the 
following: 
= Base line drift 
= Dispersion 


- Variation in solvent flowrate 


Base Line Drift 

It is noted that the number average molecular weight is 
affected considerably if the base line drift is not 
accounted for in the molecular weight calculations. Base 
line drift is corrected by noting the initial base line 
value and the final base line value corresponding to the end 
of the SEC peak and assuming a linear drift from the initial 


to the final base line value. 


Normalization 
Since the chromatogram above the base line is only 


needed in the calculation of molecular weight averages, the 
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base line values are subtracted from the corresponding SEC 
readings. The area of the resulting chromatogram is 
calculated and a normalized chromatogram is obtained by 
dividing each reading by this area. The normalized 
chromatogram has the area of unity and simplifies further 
calculations to get molecular weight averages. In _ the 
normalization procedure it is assumed that the detector 


response factor is independent of molecular weight. 


Dispersion Correction 

Axial dispersion of solute molecules of a= single 
species results in the elution of a single species occurring 
over a range of elution times, giving a distribution of 
elution times Por each species of the solute. The 
chromatogram of the sample as measured by the detector is 
the superposition of these distributions. Interpretation of 
a chromatogram must therefore account for this 
Superposition. 

If there were no column dispersion, i.e. at infinite 
resolution, the weight concentration and the molecular 
weight at any SEC elution volume, y, are W(y) and M,(y) 
respectively, where M,(y) is the true calibration 
relationship. If we let a Gaussian function G(v-y) describe 
the fraction of species at y which gets spread over to v, we 
can consider the molecular weight mixture at any elution 
volume v to be made up by species coming from all different 


y elution volumes. The molecular weight of these species are 
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M; and their weight concentrations detected at v are 
W(iy)G(v-y). The total detector response observed at elution 


volume, v, iS given by 


co 


F (v) = J Wly)G(v-y)dy (3.9) 
—c 

The above is the well known Tung's equation for axial 
dispersion and has the form of the convolution integral as 
shown by Hamielec and Ray (1969). The solution or 
deconvolution of the above equation is needed to correct for 
the effect of dispersion. A number of methods have been 
proposed (Balke et al.(1969), Provder and Rosen(1971), Yau 
et al.(1977) ) to correct for the dispersion effect in SEC 
and a method proposed by Yau et al.(1978) is used in this 
study. The details of the mathematical derivation of this 
method is presented in Appendix B. 

Narrow PMMA molecular weight standards were used to 
obtain o? values of SEC chromatograms. However, it should be 
noted that the observed variances of SEC chromatograms are 
not only due to the column dispersion but also due to the 
polydispersities of the narrow standards. The PMMA standards 
obtained from the manufacturer were characterized only for 
the peak molecular weights, whereas they were not accurately 
characterized for their polydispersity values. Hence, it was 
not possible to accurately take account of the variance 
contribution of the polydispersities of the narrow standards 


in calculating number and weight average molecular weights. 
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In this work, dispersion correction was done based on 
the observed variances of SEC chromatograms of the narrow 
standards only. If the variance contribution of the known 
polydispersity values were taken into account, it should 
change the results of number average and weight average 
molecular weights slightly. A computer simulation experiment 
was done uSing a raw SEC chromatogram to study the effect of 
variance valueS on number and weight average molecular 
weights. For a variance value of 0.40, the number and weight 
average molecular weights, and polydispersity values 
obtained were 8.23x10*, 1.35x10°5 and 1.65, respectively. 
When the variance value was changed from 0.4 to 0.0, the M,, 
Mw and polydispersity values obtained were ASSExX HOST 
1.40x10*° and 1.75, respectively. Thus the changes in M,, Mw 
and polydispersity values were 4.8, 3.7 and 11. percent, 
respectively. We observe that the variance affects the 
polydispersity values considerably. However, in our 
calculation of number and weight average molecular weights, 
the errors due to not taking account of the polydispersity 
of the narrow samples will be less than the error values 
just reported, because the o* values due to column 
dispersion will be only smaller than the o? values of SEC 
chromatograms of the narrow standards, but will be never 
zero. 

Correction for Solvent Flowrate Variations 
Since the molecular weight calculations are sensitive 


to the variations in the flowrate of the carrier solvent, 
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the average flowrate for each run is measured and used in 
the calculations. A syphon that would discharge its contents 
when it is full, is fitted at the outlet of the SEC columns. 
At the time of discharge of the SEC solvent from the syphon 
a digit switch is closed. The average flowrate is easily 
calculated from the volume of the syphon and the readings of 
the computer times at which the syphon discharged the 


solvent. 


3.4.3 Molecuar Weight Calibration 

The traditional technique for relating elution volume 
to molecular weight has been to determine the elution 
volumes of sample with relatively narrow MWD that have been 
characterized as to molecular weight by independent means. 
The most commonly used standards are anionically polymerized 
polystyrenes, prepared by Pressure Chemicals 
Inc.(Pittsburgh, Pa.). Since the relation between the size 
and molecular weight strongly depends on _ the polymer 
structure and it is the size that determines elution volume, 
a calibration curve obtained with a particular polymer will 
not be accurate for a polymer with a different structure. 
Since narrow MWD sample are only available for a limited 
number of polymer types, means for generating a universal 
calibration curve that would be applicable to _ others, 
independent of their structure, has been the subject of 
considerable study (Boni(1976)). Hamielec(1976) and Yau et 


al.(1979) discussed in detail the methods that can be used 
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to calibrate the SEC. 

Since marrow MWD standards are available (Polymer 
Laboratories FAC. ; Stow, Ohio), the molecular weight 
calibration curve is established directly by injecting 
different standards into the SEC and finding peak elution 
volumes for a particular set of operating conditions. In 
this work the SEC iS operated under the following 


conditions: 


Columns su-Styragel Columns 10°, 10*, 10° °A 
Solvent : Tetrahydrofuran(THF) 
Flowrate Th 7 a 1 


RI detector: 8x 

Six narrow MWD standards (molecular weights ranging 
from 45,000 to 680,000) were used in this study to obtain 
the molecular weight calibration curve for the SEC. In Table 
3.2 peak elution volumes (along with duplicate values) and 
the variance of the SEC peaks of the narrow standards are 
reported. A semilog plot of molecular weights and elution 
volumes shown in Figure 3.6 suggests that a semi-logarithmic 
relationship of form given in Equation (3.10) can be fitted 


to the data, 


lit Me )=ao. + a.V (3.10) 


where a, andaj, are constants and v is the elution volume. 
In fact the solid line shown in the Figure 3.5 is the least 


squares fit of all the data. 
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Fig. 3.6 SEC Calibration 
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Table 3.2 Molecular Weight Calibration 


Molecular Retention Variance 
Weight Volume (m1) o? 
45200 23256 043316 

Zoo 0.4084 
72000 22.82 0.3390 
22280 0.3753 
96000 2220 O..4247 
Zola 0.3987 
128000 21.56 023527 
Ze 0.3663 
280000 20. 16 0.3485 
Pl Vea (is: areca eee s: 
640000 18.52 0.3782 
iowa! 0.3898 


It should be mentioned here that PMMA molecular weight 
Standards are available in the range 45,000 to 680,000. 
Hence it is important to establish whether the set of SEC 
columns could be used for molecular weights below 45000 and 
whether the semi-logarithmic relationship could be 
extrapolated for molecular weights less than 45000. For this 
purpose, we injected polystyrene molecular weight standards 
ranging in molecular weights from 2000 to 200,000. The 
molecular weight-elution volume relationship for the 
polystyrene standards are shown in Figure 3.6. From this 
relationship, it can be concluded that the SEC columns were 
able to resolve molecular weights as low as 2000. According 


to the manufacturer's (Waters Associates) specification, the 
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Set OL columns 50, 10* and 10° 8A) could be used for the 
molecular weights from 500 to 8,000,000. 

From the elution volume and molecular weight plot shown 
in‘ Figure 3.6 for polystyrene standard, we can fit a 
Straight line for the molecular weights ranging from 2000 to 
200,000. Even though we do not have narrow PMMA standards 
for molecular weights lower than 45000, we can reasonably 
assume the semilogarithmic relationship for molecular 
weights up to 2000 by analogy with polystyrene standards. Of 
course, we can expect that this assumption can lead to some 
error in the calculation of number average and _ weight 
average molecular weights particularly at molecular weights 
less than 45000. 

Two typical SEC chromatograms, one corresponding to a 
high molecular weight sample (M,=8.83x10‘, Mw=1.48x10°) and 
the other corresponding to a low molecular weight sample, 
(M,=3.04x10*, Mw=5.74x10*) are shown in Figure 3.5. We can 
observe that for the low molecular weight sample, the 
elution volumes range from 20 mL to 29 mL. Since the elution 
volume for the 45000 molecular weight sample of PMMA is only 
23.6 mL, we note that extrapolation of the molecular 
weight-elution volume relationship is necessary for elution 
volumes greater than 23.5 mL. We could have used a universal 
calibration curve based on polystyrene standards to avoid 
extrapolation. 

Usually the molecular weight-elution volume 


relationship is expressed in the form: 
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M, = D,exp(-D2v) tse) 


From Equations (3.10) and (3.11) the constants D, and D2 are 


obtained as, 


D,=exp(ao), D2=-a2 (32872) 


The D, and and Dz values obtained from Equation(3.12) 
along with the reported value of average variance in Table 
3.2 are used to calculate the number and weight average 
molecular weights as shown in Appendix B. 

It should be pointed out that the SEC calibration 
should be repeated often to compensate for column 
deterioration. While small fluctuations in temperature or 
flowrate may not upset the calibration, changing columns, 
solvent, or the nature of the polymer samples will require 
recalibration of the SEC experiment. After several months of 
use, we observed only a slight change in the calibration 


curve of the SEC used in our study. 


3.5 SEC Error Analysis 

A number of experiments have been conducted to check 
the reproducibility of SEC results. A sample of the 
polymerization reaction mixture diluted with the solvent THF 
was injected seven times into the SEC and the number and 
weight average molecular weights obtained are reported in 
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From Table 3.3 it can be observed that the averages of 
the number and weight average molecular weights are 42000. 
and 76400., respectively. The percent relative standard 
deviation for the number and weight average molecular 
weights are calculated as 5.3 and 4.1, respectively. These 
results are comparable to the reproducibility studies 
reported by Balke(1973). 

The accuracy of molecular weight average determined 
from a SEC experiment is dependent on factors which are 
discussed by Yau et al.(1979). It is important that the 
samples are diluted enough to prevent the overloading of the 
columns and Significant band broadening due to viscous 
streaming. A rough guide for sample dilution is that an 
injected sample should have a viscosity no greater than 
twice the viscosity of the mobile phase. For high molecular 
weight polymers concentrations less than 0.10 weight percent 


are often required to eliminate undesirable effects. 


3.6 Viscosity Measurement 

Even though number average and weight average molecular 
weights can be measured using the SEC, there is a delay of 
10-15 minutes before these meaSurements become available to 
take any control action during the course of polymerization. 
Hence it is desired to have more rapid on-line measurements 
of the molecular weight averages. In this context the 
viscosity of the polymerization mixture has been measured to 


get an estimate of the molecular weight average of the 
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Table 3.3 Reproducibility Results of SEC 


Injection Number Average Weight Average 


Number Mol. Weight Mol. Weight 
1 3.9584 7.50EB4 
2 3.97EB4 7.31B4 
3 4.33E4 7.70E4 
4 4.4554 8.04E4 
5 4.13E4 8.15E4 
6 4.01E4 7.48E4 
a 4.544 7.31E4 
Average 4.20E4 7.64E4 
Std. deviation 2e23s SeLSEs 


polymer. 

It is known that the viscosity of linear polymer 
solutions are related to the molecular weight of the polymer 
(Flory(1953)). Since PMMA is a linear polymer an expression 
relating the molecular weight and the viscosity of the 
polymerization mixture could be developed. In addition to 
molecular weight, the viscosity of the polymerization 
mixture is also affected by monomer conversion, temperature 
of the reaction and concentration of the solvent. Influence 
of all of the above variables on the viscosity of the 
reaction mixture should be taken into account when 
developing mathematical models to estimate the molecular 
weight average from measurement of viscosity. 

Jo(1975) measured the viscosity on-line for the case of 
solution polymerization of vinylacetate in a CSTR and 


estimated the weight average molecular weight from the 
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following correlation: 


In(n)= b, +b2XMw" C3013) 


In the above equation 7, X, Mw are the viscosity, monomer 
conversion and weight average molecular weight respectively. 
The constants b, and bz and n had been determined from his 
experimental data. 

Alhough a number of on-line viscometers are available 
in the market (Huskins(1982)), only a few of them can be 
considered suitable for use in the laboratory-sized reactor. 
Most of them required large sized process streams or large 
sized reactors, the emphasis of application being on 
industrial-scale processes. Finally, an industrial-type 
viscometer supplied by Norcross Corporation which could also 
be used for laboratory application was chosen for this 
study. 

This viscometer operates on a falling piston principle 
which is better explained by referring to Figure 3.7. The 
measuring element consists of a cylindrical tube which is 
mounted vertically on the reactor lid. The lower part is 
immersed in the reaction mixture and a piston-rod assembly 
is located inside the measuring tube. During the filling 
phase the piston-rod assembly is lifted by an air-actuated 
mechanism drawing a sample of the polymer mixture down 
through the clearance between the measuring tube and the 
piston. During the measuring phase, the piston rod assembly 


is allowed to fall under gravity, expelling the same sample 
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through the holes in the measuring tube. The time required 
for the piston-rod assembly to drop to the bottom of the 
tube iS a measure of viscosity and this measure is available 
imerne range “of 1-5V for on-~line-> data acqusition. The 
element is continuously self-cleaning owing to the reverse 
flushing that takes place during the normal up-and-down 
motion of the piston. In this instrument viscosity 
measurement is taken every minute which is _ considerably 
faster than the 10-15 minutes analysis time for SEC 
measurement. 

The viscosity range that can be measured uSing a single 
DisStoniseehimeted-—tom a» asingles decade ««(ise«s ep-10cp, 
10cp-100cp etc.). Since during the course of polymerization 
the viscosity can increase from icp to as high as 100cp, it 
is necesSary to change the piston when the viscosity values 
exceed the maximum range. However the piston can be changed 
easily in two to three minutes time which is negligible 
compared to the total reaction time of 1-5 hours. 

The viscometer is calibrated for each piston separately 
uSing viscosity standards available in the market(Can-Am 
Instruments Ltd, Mississauga, Canada). The calibration curve 
obtained for two pistons are shown in Figure 3.8. The 
computer readings are converted to viscosity values using 
this calibration curve. In the next chapter the viscosity 
measurements along with the analysis to estimate the 


molecular weight of the polymer will be presented. 
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Fig. 3.8 Viscometer Calibration 
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3.7 Torque Measurement 

Since the viscosity of the polymerization mixture 
increases during the course of polymerization, the torque 
exerted on the constant RPM stirrer by the reaction mixture 
also increases. By measuring the torque exerted on the 
Stirrer, it is possible to estimate the molecular weight 
average of the polymer produced during the course of 
polymerization. In fact, Jo(1975) monitored on-line the 
torque on the stirrer in order to measure the viscosity of 
the reaction in a solution polymerization process. In turn 
the viscosity measurement was used to estimate the weight 
average molecular weight of the polymer as described in the 
previous section. 

The stirrer-motor assembly used in this experimental 
study has been already described in chapter two. The torque 
exerted on the stirrer is available in the range of 0.10V 
output for on-line data acquisition. Since torque 
measurement is continuous, molecular weight averages could 
be estimated instantaneously if a reliable correlation 
relating torque and molecular weight of the polymer is 
developed. 

Though the torque measurement could have served as an 
additional measurement for estimation of the molecular 
weight of the polymer produced during the reaction, in this 
study the torque measurement could not serve that purpose. 
The torque measurements were found to be noisy due to 


unpredictable friction at the bearing site of the stirrer. 
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Since it was not possible to predict the torque exerted on 
the stirrer only due to the polymerization mixture, these 


measurements were not analysed to develop correlations. 


3.8 Conclusions 

Various on-line as well as off-line techniques to 
measure monomer conversion, molecular weight averages, 
viscosity and torque on the stirrer for a free radical 
solution polymerization process have been described in this 
chapter. Comparison of conversion measurements from 
gravimetry, gas chromatography and densimetry confirms that 
gas chromatography and densimetry results are more accurate 
than gravimetric results. Further it has been concluded that 
gravimetric conversion results are prone to error especially 
at higher levels of conversion and hence these results 
should be accepted only after careful consideration. 

A description of computer software that was developed 
for on-line acquisition of the SEC data and analysis to 
calculate the number and weight average molecular weights of 
the polymer is also given. This software along with the 
automatic sampling system form an important contribution to 


the automatic analysis of polymerization reaction mixtures. 
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4, EXPERIMENTAL RESULTS 


4.1 Introduction 

In chapters two and three the experimental reactor 
system and the techniques of measuring various process 
variables in the solution polymerization of MMA in a batch 
reactor were described in detail. It is important to 
generate experimental data on the polymerization processes 
for the following purposes. Various on-line measurements 
such as density, viscosity and SEC can be analyzed to 
develop correlations to estimate important process variables 
such aS conversion, mumber and weight average molecular 
weights. These experimental results can also be used to 
check the validity of the mathematical model of the reactor 
as will be shown in chapter five. Also by designing 
appropriate experiments we can understand the effect of 
various inputs such as temperature and initiator 
concentration on the polymerization process. 

A total of 15 experiments were conducted with different 
combinations of initiator concentrations and temperatures. 
The experimental procedures described in section 2.5 were 
followed to measure process variables such as _ temperature, 
density, viscosity, torque, number average and weight 
average molecular weights. All of these variables were 
recorded on the HP/1000 computer through an on-line data 
acquisition program. The experimental data obtained are 


Summarized and reported in Appendix G. For experiments (1-7) 
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the number and weight average molecular weights are not 
reported because of the experimental difficulties with the 
SEC in getting reliable values. Some of the SEC columns 
Plugged and had to be replaced with new columns. Later it 
was found that some of the earlier measurements done using 
the old columns could not be reproduced well and hence they 
were discarded. Although torque measurements were recorded, 
they were not reported because of the unreliability of these 
data, as discussed earlier in section 3.7. In the following 
sections the analysis of the experimental results will be 


presented. 


4.2 Reproducibility of the Experimental Data 

The results from a polymerization experiment are 
dependent on many different factors: purity of monomer and 
solvent, the condition of the initiator, the presence of any 
inhibitor”, inthe réaction mixture} ~the’ “purity °of° “the 
nitrogen blanket, the actual temperature of the reaction 
mixture, etc. A good reproducibility of the experimental 
results gives confidence in the accuracy of the results in 
spite of the difficulties in maintaining identical 
experimental conditions. Hence it is important to repeat 
some of the experiments keeping identical conditions to 
check if the results are reproduced at all. This lends 
credibility “to the experimental set-up as well as the 


measurement techniques and the instruments. 
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Figure 4.1 shows the reproducibility results for 
conversion, viscosity, mumber average and weight average 
molecular weights for a case where the temperature, solvent 
ratio and initial initiator concentration were 70°C, 0.6 and 
0.05 mol/L respectively. From this figure it can be observed 
that we are able to duplicate conversion results very well, 
the maximum difference being only 5 per cent. There are 
differences to the” order of 10-15 pen. cent ~in® the 
reproducibility of number average and weight average 
molecular weights and even larger differences of the order 
of 20-25 per cent can be observed in the reproducibility of 
viscosity values. The same trends are observed in Figure 4.2 
in which reproducibility results are presented for the 
following experimental conditions: temperature=75°C, solvent 
ratio=0.6 and initial initiator concentration=0.10.mol/@. 

From the above results i1t can be .concluded “Sthat 
conversion, number average and weight average molecular 
weights can be duplicated within reasonable accuracy when 
identical experimental conditions are maintained. However, 
viscosity results could not be reproduced well enough and 
the reason for this can be explained as follows. The 
viscosity of the polymerization mixture depends on monomer 
conversion, molecular weight of the polymer and also the 
temperature of the polymerization mixture. Since the errors 
in any of the variables affect the viscosity, the poor 
reproducibility is most probably due to _ the cumulative 


effect of small errors in all of these variables. 
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4.3 Effect of Solvent Concentration 

A number of experiments were conducted to study the 
effect of solvent concentration on the rate of 
polymerization in the batch reactor system. These 
experiments are helpful to decide whether the gel effect, in 
which the rate of polymerization increases with increasing 
conversion is significant, and thus whether it should be 
taken into account in developing the mathematical model of 
the reactor. In Figure 4.3 the conversion-time histories 
obtained for solvent volume fractions of 0.4, 0.5 and 0.6 
(temperature=75°C I,9=0.05 mol/L) are presented. No marked 
differences can be observed in the conversion-time histories 
for different solvent fractions. The same result was also 
observed for experiments conducted at 65°C for different 
solvent fractions. The shape GF the conversion-time 
histories in Figure 4.3 shows a decreasing slope of the 
curves throughout the course of polymerization which means 
the rate of polymerization decreases during course of the 
reaction. If a pronounced gel effect were present, the rate 
of polymerization would increase with increase in 
conversion. From these experimental results, we can neglect 
the gel effect phenomenon in the modeling of the 
polymerization reactor. The gel effect can be neglected in 
this work because of high solvent fractions (0.6) employed 
and also because of the limited extent of conversion at 
which the viscosity of the polymerization mixture is not 


very high. On the other hand it is known that in bulk 


i 


eos q2ata, eS Gedwubnes Siew © ner tags: ° 
tet tt Be) 4 90 (exeaant i dieitoe 
wedeys TESST Avda oad | 

tie tap ety Ieditsdy shboat ba ine ona 

a Lewsnigt ; Adee ambinn soti ere Bt . 
h Cusde  -edsere avis bee Gabsphingte ee ofa: | 

io febor te iam vat y eriqotewe® tb drwospa 
vein oni a-aoheasiwne. ods Suh grupis at 1s | 

14 to enehote+) erutey hi oleae ii bee : 

ye» ihe / nap erieg ou6 atte Bcd 5*2y % 7 pr 


oa 11289 wn att’. deal sno ia 
pemerers, wm? fsa 28 Be soub«o? “eifiimtssaxe 508 
wit s-nolepewnoe!” aie | Ip eget, pat “ethottoe37) | 
to siete sidaileitie avout eb iy ng | 
angot doldy roi texbpamytog te araiss a swedguesda 
fe sans coer ban abapn se cotdqed aga fog 30 SSRD 
aoe" | ot? tania aw soetts twp Besoviono7q 6° 5 tM <node 
rr! senmaaete’ go ee ehnwson -blinow nolsssivemgfog. a 
froiger, “wes yu aghoass titoomtegne head? ao7t,, coterevao> ay x 
wit. do gai Calne ays nl consaeneng spetie fee | ‘io a 
ni paszelgen ac Soo anette top ont Vaodoae7 noiseaisemytog | 
bevel pie £8 1) anecsoan? aoawins styitit te envased “tow eind 
1 folsveunso” to diberae bed imi eas Yo savin cals bas 
ton ei stein 062 1881 tomy hog pers +9 gitsovty | 9d fotde 


divd a say nwaek ak we) Bead teddo' odd nO 2Abie Yrev 


oe 


83 


UOTPEZLIaWIATOY JO ayeY UO OTWeY YUAAIOS Jo JOaTTY EP By 
SOYNUIP, UT OUIT, 
0-001 0°08 0°09 0'0F 0:02 00 


0'09T oO Ort 0O0eT 


SO’ 

CO: 

SO’ 
T/1OW 


0 Ss 


I 


O0;00T 0:06 O08 O0LZ 009 00S O0F ONDE ONDE ODOT O00 


UOISI3AUO) 


ba a us 


WELtisgriov 


OFAGHE 


A 
a 


EIS FETCH OF : 


exypTo Of 


+ 
? 
- 


ag Oe, fies a ot cond 


84 


polymerization of MMA, the polymerization reaction shows gel 
effect at conversion as low as 10 per cent (Odian(1981)). 
The presence of solvent in the polymerization mixture 
reduces the viscosity of the medium and postpones’ the 


occurrence of gel effect to high conversion values. 


4.4 Effect of Initiator Concentration 

A series of experiments were conducted to study the 
effect of initial initiator concentration on the course of 
polymerization. Two levels of initial initiator 
concentration (0.05 and 0.10 mol/L) and three levels of 
temperatures (65°C, 70°C and 75°C) were employed to generate 
the necessary experimental data. Figures 4.4 to 4.6 show the 
effect of increasing the initial initiator concentration on 
the conversion-time histories. It has been established that 
the polymerization rate is directly proportional to the 
Square root of initiator concentration (Odian(1981)). Hence 
increasing initiator concentration should increase the rate 
of polymerization. This can be clearly observed in Figures 
4.4-4.6. 

Also Figures 4.4 to 4.6 show the effect of initiator 
concentration on the development of number average and 
weight average molecular weights, respectively. Sacks et 
al. (1973) proved that in conventional free radical 
polymerizations taking place at constant temperature, the 
instantaneous value of number average molecular weight of 


the polymer produced would decrease during the course of 
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Fig. 4.4 Effect of Initiator Concentration on Polymerization at 65°C 


85 


7 


Soe 


0 oes P o9oe ne Pepa 4 oos : 
! | eos ote 


Pi: 


86 


Mew oe 
rm © 
x 

| 

S 

Gi 

s satus 

° 

oO ie) 

So 

~ 

ie a eel 
0.0 50.0 100.0 150.0 200.0 250.0 

~ 6° Time in Minutes 
=) 
xe tl 

° 

Sa Oo 

° oO O a) 
se Oo 0 O 

et) 

°o O 

+ 

o 

ch A  eaiaettse Une Sen 

2 res SS See ee ie FS eee Oe 
0.0 50.0 100.0 150.0 200.0 250.0 
Time in Minutes 
@ 
° 


mS 
T tO ¥el, 
4 © Mol/l 
S O= 7012.68.05 
O=70. .6 .10 
° 
° 


100.0 we. <*s0me | Yoom.e | asat o $6 2oomie 25 0 250.0 
Time in Minutes 


<3 
Fo) 
So 
° 
v 
so 
nN 
° 
° 
°o 
; 00.0 


ae yt 


Fig. 4.5 Effect of Initiator Concentration on Polymerization at 70'S 
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polymerization, when transfer to monomer and solvent do not 
control molecular weight. The same trend would be observed 
by the weight average molecular weight. Also, the kinetic 
chain length, and hence degree of polymerization, is shown 
to be inversely proportional to the initiator concentration 
(Odian(1981)). Hence with increase in ind trator 
concentration, the number average as well as weight average 
molecular weights should decrease. We observe in Figures 4.4 
to 4.6 that with increase in initiator concentration, the 
number and weight average molecular weights decrease. 
Similar experimental results for number and weight average 
molecular weights have been reported by Chen and Huang(1981) 
for the solution polymerization of styrene in a batch 
reactor. 

A master plot is drawn by plotting conversion versus 
(/Io)t for different temperatures as shown in Figures 4.4 to 
4.6. Since the data, "peimes Siror “<difterene initiator 
concentrations coincide, we can infer that the rate 
polymerization is proportional to the square root of the 


initiator concentration. 


4.5 Effect of Temperature 

A number of experiments have been conducted in order to 
study the effect of temperature on different variables of 
the polymerization process. The effect of temperature on the 
rate of polymerization is shown in Figures 4.7 and 4.8. It 


can be observed from these figures that the rate of 
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polymerization increases with increase in temperature. 

The effect of temperature on the rate of polymerization 
is determined by the value of activation energy for the rate 
of polymerization. The activation energy for the rate of 
polymerization can be derived as (EptE,/2-E,/2) where Ep, 
Eu, and E,~ are activation energies for propagation, 
dissociation and termination reactions, respectively. The 
Ep, Eg and BE, values for MMA polymerization are 9600, 33700, 
1000 cal/mol, respectively. From these values, the 
activation energy for rate of polymerization is calculated 
as «25,500 cal/mols Hence, it (is. «clear. that. the. rate. of 
polymerization will increase with temperature. This can be 
observed in Figures 4.7 and 4.8. 

The effect of temperature on number average and weight 
average molecular weights is also shown in Figures 4.7 and 
4.8. Odian(1$81) has shown that except for photochemical 
polymerizations where the degree of polymerization increases 
with temperature, in all other cases the degree of 
polymerization and hence the molecular weight averages 
should decrease with increase in temperature. This decrease 
in number and weight average molecular weight can be clearly 
observed in Figures 4.7 and 4.8. 

These experimental results establish that initiator 
concentration and temperature significantly affect the rate 
of polymerization, number and weight average molecular 
weights. These two variables can be employed to control 


conversion and molecular weight averages. Open loop control 
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policies using initiator concentration and temperature as 


control variables are developed in chapter six. 


4.6 Analysis of Viscosity Data: 

In chapter three the method of meaSuring viscosity 
on-line uSing an industrial viscometer was described. 
Viscosity values were measured at various temperatures and 
initiator concentrations during the course of 
polymerization. Figures 4.9 to 4.11 present the experimental 
values of the viscosity of the polymerization reaction 
mixture. From these figures it can be clearly observed that 
at the same conversion and temperature, the viscosity values 
are less for the experiments which employed higher initiator 
concentration. It has been explained (see section 4.4) that 
at higher initiator concentrations, lower molecular weight 
polymers are produced. Hence the effect of molecular weight 
on the viscosity is clearly indicated in these figures. 

Viscosity is a function of many variables. Though 
mathematical equations have been developed relating 
molecular weight for dilute solutions of the polymer, for 
concentrated solutions such relations are lacking 
(Rodriguez(1982)). 

Jo(1975) developed a semiempirical model to estimate 
the weight average molecular weight of the polymer from the 
measurement of viscosity during the course of 
polymerization. Since all the experiments were conducted at 


a constant temperature, they did not take into account the 
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variation of viscosity with temperature in their model. 
Recently Hamer (1983) developed two types of 
correlations to estimate viscosity of the polymerization 
mixture in solution polymerization of vinylacetate. 
The first and the simpler version was based on the 


correlation suggested by Spencer and Williams(1947). 


n = A exp[ F(X,;E) ( + CY(.001Mw + D)] (4,1) 
RT 
where 
P(x GE ye ak oR Rye 4B? 
=X /4E, Mes eB 


where xX, is the weight fraction of polymer in solution. A, 
Bepe, i DilandusEigare constants. Equations of this form, 
although difficult to support theoretically were found to be 
usefuolesitor) bfiteine early experimental results with 
concentrated polymer solutions. (Flory(1940), Ferry et 
ad (1951008 

The second viscosity correlation was based on the free 
volume concepts and chain entanglement theories discussed in 
detail by Berry and Box(1968), Graessly(1974) and 
Ferry(1980). We will not deal with the details of the 
mathematical derivation here aS our aim in this study is to 
develop empirical relationships that can be used to _ obtain 


molecular weight values from measurements of viscosity. 
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Since in this study the objective is to estimate the 
molecular weight from the measurement of viscosity, the 
empirical correlation(Equation 3.13) used by Jo(1975) was 
modified to take account of the variation of viscosity with 


temperature. The following correlation was used. 


a2 ay 
) at (a3 ae 
db ah 


)XYM, (452) 


where 7, T, X and M, are viscosity, temperature, conversion 
and number average molecular weight respectively. The 
ParameterS a;,, a2, a3, and ay, must be estimated from the 
experimental data. 

Viscosity values obtained from six experiments were 
used to eStimate the four parameters appearing in the 
empirical equation. 

Figure 4.12 shows the comparison between the 
experimentally measured viscosity and the viscosity values 
estimated from the Equation 4.2. From these results it can 
be observed that the model predicts the viscosity reasonably 
well. The correlation is reasonably good at least in the low 
viscosity region. Equation (4.2) can be rearranged to get M, 


from the measurement of viscosity as: 
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Fig. 4.12 Comparison of Measured and Predicted values of Viscosity 
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t inta)i=-a,-a;/T)}? 
mM, = { -—————————} (4.3) 
Cher tag TK) 4 
The experimental M, values and the M, values predicted from 
Equation (4.3) are plotted in Figure 4.13. From this figure 
it can be observed that M, is predicted well at high 
viscosity solutions, whereas considerable discrepancy exists 
at molecular weights predicted for low viscosity solutions. 

From these results it can be deduced that molecular 
weight of the polymer can be determined on-line from 
measurement of viscosity, conversion and temperature. We 
have observed that the M, is estimated with better accuracy 
at higher viscosity values whereas in dilute solutions 
having lower viscosity the estimation of M, is poor and 
unreliable. 

Our analysis of viscosity data is not exhaustive and 
complete by any means. We tried various forms of 
semi-empirical functions to estimate number average 
molecular weight from viscosity measurement. However, none 
of them could estimate the number average molecular weight 
with better accuracy over the whole viscosity range. 
Perhaps, application of free volume theory to the viscosity 
data could have provided better estimates of molecular 
weight. 

It should be pointed out here that the control action 


is not completely dependent on the molecular weight 
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estimated from the viscosity measurement. These estimates 
could be updated when the molecular weight averages become 
available from SEC measurement. This combined viscosity and 
SEC measurements of molecular weight averages could provide 
better control of the polymerization reactor. As an 
alternative recently introduced, Bondagel columns by Water 
Associates are able to perform SEC analysis of polymer in 
about 5 minutes. With these columns, on-line molecular 
weight determination using SEC has become possible in very 


short analysis time. 


4,7 Torque Measurements 

Figure 4.14 shows the torque meaSurement obtained using 
the experimental reactor system. The torque exerted on the 
Stirrer by the polymerization mixture increases during the 
course of polymerization. The torque meaSurements were done 
with an intention that they would provide additional 
meaSurements to estimate viscosity and therefore indirectly 
the average molecular weights. However these measurements 
turned out to be unreliable because of the friction at the 
bearing site of the stirrer. Hence these measurements were 


not analyzed to estimate the molecular weights. 


4.8 Conclusions 
In this chapter extensive experimental data obtained 
using the reactor have been presented. The analysis of the 


data revealed that initiator concentration and temperature 
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Gfie the reaction considerably affect the rate of 
polymerization and the number average and weight average 
molecular weights of the polymer formed during 
polymerization. Hence these two variables, temperature and 

initiator concentration, can be employed as manipulated 
variables to control monomer conversion and molecular weight 
averages in the polymerization reaction. The analysis of 
viscosity data showed that at a higher viscosity of the 
polymerization mixture, the M, values are predicted well 
whereas considerable discrepancy exist in M, values 
predicted in dilute solutions of the polymer, where 
viscosity is also low. The experiments on reproducibility of 
the results showed that conversion and molecular weight 
averages are reproduced well when identical experimental 
conditions are maintained. This lends credibility to the 


instruments and measurement techniques used in this study. 
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5. MODELING AND PARAMETER ESTIMATION 


5.1 Introduction 

A reliable mathematical model of the polymerization 
process is required in order to calculate the optimal 
control policies that will improve the performance of the 
batch reactor. Free radical polymerization of vinyl monomers 
is a well studied process and a description of the 
elementary chemical reactions occurring in the process is 
readily available in most standard texts on polymer 
chemistry (Odian(1981)). The more fundamental properties, 
often of importance to the chemical engineer, such as 
conversion and molecular weight distribution are the result 
of elementary reactions occurring in the polymerization 
process and their relative rates. By developing a 
mathematical model of the process from these elementary 
reactions, it is possible to find the effect of control 
variables such as temperature and initiator concentration on 
monomer conversion and molecular weight distribution(MWD) 
during the course of polymerization. 

In the past, a number of investigators have reported a 
variety of models to describe the free radical 
polymerization of methylmethacrylate (MMA) in’ alebatch 
reactor. Matheson et al.(1949) measured the average life 
time of PMMA radicals as a function of temperature in 
photosensitized polymerization of liquid monomer. By 


combining their experimental results with those of Schulz 
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and Blaschke(1942), and Schulz and Harborth(1947), they 
obtained the rate constants for propagation, termination and 
transfer to monomer for conversions up to 10 percent. They 
attributed the accelerated polymerization rate occurring in 
the latter stages of the reaction to a decrease in the 
termination rate constant. 

Nandi(1957), Ferington and Tobolsky(1957) have reported 
experimental results that were obtained with a wide variety 
of initiators under different isothermal conditions for bulk 
polymerization of MMA. Using their experimental results they 
calculated the rate constant for chain transfer to monomer 
and the ratio kp/VYk, where kp and k, are the rate constants 
for propagation and termination reactions, respectively. 
Their results were in good agreement with the results of 
Matheson et al.(1949). 

Hayden and Melville(1960) measured the increase in 
temperature in an adiabatic polymerization reactor as a 
function of time and reported that after about 10 percent 
conversion the polymerization rate and the life time of the 
polymer radical chains increased. This increase was 
attributed to the increasing viscosity of the reaction 
mixture resulting in a decrease of the collision rate of the 
growing radicals. They found that beyond 40 percent 
conversion, the activation energy for the propagation 
reaction also increased while the corresponding velocity 
coefficient decreased due to the monomer addition step 


becoming diffusion controlled. Paul et al.(1973) also 
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observed a Similar decrease in the propagation rate 
constant. 

Balke(1972) investigated the bulk free radical 
polymerization of MMA to high conversions. Using size 
exclusion chromatography (SEC) and gravimetric analysis he 
followed changes in the MWD and conversion during 
polymerization and developed a mathematical model for this 
process. 

Mahabadi and Meyerhoff(1979) proposed a new model to 
estimate the rate constant for primary radical termination 
using the radical life time, the rate of polymerization and 
the rate of initiation. By applying this model to high 
conversion polymerization experimental data, they showed the 
dependence of conversion on termination and initiation rate 
CONnSLanes. 

Schmidt and Ray(1981) developed a mathematical model to 
describe the free radical solution polymerization of 
methylmethacrylate in a CSTR. They also developed 
expressions to account for the gel effect based on the free 
volume correlation of Ross and Laurence(1976). 

A mathematical model to describe the free radical 
solution polymerization of MMA in a batch reactor is 


developed in the following sections. 
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5.2 Description of the Reaction Mechanism 
A general description of the chemical reactions taking 
place during the free radical solution polymerization of MMA 


initiated by benzoylperoxide is as follows: 


Initiator Decomposition: 

I kg 2£Ro 
Initiation: 

Ro + M kj R, 


Propagation: 


Rx + M kp Ryo 
a 


Transfer to Monomer: 
R, + M Km Py ot 4 
—____> 
Transfer to Solvents 


Ry + 9 Key Pee on4 
_ 
Termination by Disproportionation: 
Ree ty Kea IS Ne 
SS 


Termination by Combination: 


Regan tag Fy Ree Prey 


where I denotes the initiator molecules, Rp. the initiator 
radicals, Rx and R, denote the live radicals of chain 
lengths x and y respectively. M is the monomer, and P, and 
P, represent the dead polymer of chain length x and y 
FESDECCIVELY wha pmikitn kD; km key hae, end k,care the mate 
constants for initiator dissociation, thictiation, 


propagation, transfer to monomer, transfer to solvent, 
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termination by disproportionation and termination by 
combination respectively. 

Using the above kinetic scheme, general population 
balance equations that describe the time variation of the 
concentration of different species in the batch reactor can 
be written. To simplify the mathematical formulation of the 
problem the following assumptions are made. 

1. All the reactions are irreversible and elementary. 

2. Reaction rate constants are independent of chain length. 

3. Reactor contents are perfectly mixed and there are no 
temperature gradients in the reactor. 

These assumptions are common and are well documented in 
the literature on the modeling of free radical 
polymerization reactions. Although chain transfer to monomer 
and solvent are usually neglected in most other modeling 
Studies reported in the literature, they are included in 
this model development. 

The volume of the reaction mixture decreases during the 
course of polymerization due to the increase in the density 
of the reaction mixture as the result of monomer conversion 
to polymer. This volume change is accounted for in the 
modeling of the batch reactor. By assuming that the reaction 
volume varies linearly with the monomer conversion, the 
volume of the reaction mixture, V, at a conversion, xX, can 


be expressed as: 
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Ve Vo (1+eX) CS) 
where Vo is the volume of the reaction mixture at conversion 
X=0 and e¢« is the fractional volume change at X=100%. The 


monomer conversion X is defined as: 


X=(MoVo - MV)/ MoVo (542) 


From Equations (5.1) and (5.2) it can be easily shown that 
the change in the reaction volume is related to the monomer 


conversion by 


| av € dx =e aM 
—_  — = —S = GS!) 
V dt eX+1 dt Mo + eM dat 
The general equation for the species, i, in the batch 
reactor can be written as: 
1 aN, ie elite) aC ; Cj dv 
R; a —-Se  eee (Gee a5) 
V at V dt dt V dt 
Combining Equations(5.3) and (5.4) we obtain: 
aC ; C; dv 
Bin R; See — (5-7 5)) 
dt V dt 


where N,; represents the number of moles of the species, i, 


in. the reactor, )C; the concentration, and R; denotes the 


‘nos ¢& sundae rol joet. 
LOU =N th, apaediy emetov, I 


Che ferry i ry nolan 
, a "7 ae ee 


eh +6 Ve | > tap hz, i 
:) | enaiNin 


f 


, % 


ie ye 
ae 
AWOTe Va ahaa at bis + (sie) ane u. e}.. a0 


vonoaem. edd’ oF kasueee ar. smu toy a one Pia e 
| La anak ie ee 


i veelvege ong fo sefem Bo tedain ad3 azrnosetqed (M sredw 
pan eetores ,f Ohe gk pe Teeoned edt -3 .yososer odd ai 


ey 
— 


— a4 


110 


rate of appearance (disappearance) of species, i. In the 
following section, the mass balance equations for the 
different species in the reactor are derived directly from 


Equation(5.5). 
5.3 Derivation of the Model Equations 

According to the above kinetic scheme, the mass balance 
equations, for initiator, monomer, live macroradicals and 


dead polymer chains are expressed as: 


ab/dt t= =at 0)=.0CF /¥) “(av seb) (5.6) 


@M/dt = -k;MRo - (kptk,)M ERs - (M/V) (dv/dat) Chin) 


dR,/dt =k MRo-kpMR)-k«Ry ER,+(k,StkqM) ER, 
=0RA/V.) oCdv/dt) (Se) 
oo 
aR /ate= =kpMR.+kpMR. 4 —CK,Mtk,S) Ry KkiRx Z Ry 
= GRy/V) Cav/det) (5.9) 


co x-1 
GP, /dt = (k,StkmM)Rx + keaR, Z Ry + (Kg eR Ree 
y= 


Yaqui 


-(P,/V) (dv/dt) (510) 


In the above equations k,; is the sum of kyy and ke. 


The last two equations represent an infinite number of 
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differential equations since x can vary mathematically from 
2 to ~. Therefore the above description of the process 
results in a model with an infinite number of differential 
equations. To reduce the order of the mathematical model we 
introduce the moments of the live and dead polymer 
distributions. 

It should be pointed out here that the concept of 
instantaneous molecular weight (Hamielec(1976)) could also 
be used to develop the mathematical model of the reactor. 
This method gives the molecular weight distribution in 
addition to the number and weight average molecular weights 
of the dead polymer. 


Define the k-th moment of the live radicals as: 


AXA HO) eet: Ry, (taht) 


oo 
i a aere APL (5.12) 


From Equations (5.6)-(5.10) we can easily obtain equations 
for the leading moments of the live and dead polymer 


distributions. 
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Ako /G0 = 2ZEk yt - kite * — (Ag/V) (avV/dt) Gon? 


da ,/dt =2£ikyIt+kpMAot+(k,StkmM) (Ao-A,) i kiAodA2 
-(X,/V) (daV/dt ) (5. 14) 


Gr2/dt = 2fkygI + 2kpMA, + (kpM+k,,Mtk.S)Ao 


= (ke eM+k JS) As - chido earm -ORe/V) (AV/at) (52 ho) 


icy dt eens PK pM) Ago tek gho & kag Ao /2 


-(uo/V) (avV/dat ) (5.16) 


Gpe/at cenGkieStk aM) eich kee Xd; 


= (it {/Mi i GAV7at) G5 5 157) 


du>/dt = (k,Stk,M)X2 + kiAgAo * kee hy? 

=(p2/V) (dv/dt) Cai.) 

The leading moments uo and uw; denote respectively the 

concentration of the terminated polymer and the total amount 

of monomer associated with the terminated polymer. The 

moments Ao and A; have a similar meaning for the growing 
radical distribution. 

From Equations (5.3), (5.7) and (5.11) we obtain the 


rate of monomer consumption as: 


GM/dt = -kpMAo(Mo + Me) /Mo (5.19) 


In Equation (5.19) it is assumed that the monomer 
consumption and monomer production due to initiation and 


chain transfer to monomer and solvent are negligible. Thus 
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the rate equation for the change of volume can be obtained 


from Equations (5.3) and (5.7) as: 


(1/¥) (dN /dt) =“kproeM / Mo (5.20) 


In order to simplify the kinetic equations a quasi 
Steady state approximation (QSSA) for live radicals is 
usually made. According to this approximation the net rate 
of production of radicals is assumed to be zero. This 
implies that the number of radicals produced by the 
initiation process must be approximately equal to the number 
of radicals consumed through the termination reactions. In 
practice it is found that QSSA is usually valid for all but 
the first few seconds of the free radical polymerization 
reactions. 

Using the QSSA assumption the LHS terms in Equations 
(5.13)-(5.15) can be set to zero. Further neglecting any 
volume change, the resulting algebraic equations are solved 


to obtain the moments of the live polymer distribution. 


Nomaay DEKE TEL Ks (oe 2) 
Mw Si ZEkGt ls] “CRDMEKEM*K S)AG) / Ck, StkmMtki Ao) (5522) 
Ne = Nyt CZ2KpMA, )/ Ck, StkpMtki Ao) (523)) 
Substituting Equations(5.21)-(5.23) into Equations 


(5.16)-(5.18) we finally obtain: 
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du,/dt = Cees or KmM) A, tk Aor 4 C5)..25;) 
du2z/dat = (k, Stk ~M)A> te Kako ks af Kaen @ (5.26) 


The total monomer consumption will be equal to the sum of 
the first moments of the dead and live polymer 


distributions. 


Mo = M = (Or or v4 (S227) 


Differentiating Equation(5.27) with respect to time and 
neglecting the contributions of the term A, compared to u,, 


we obtain: 


ane, /arl =) Moo (Vax / dts) (5.28) 


From Equations(5.28) and (5.19) monomer conversion rate can 


be expressed as: 


ax kpuls2) 
=~f2£k,I ——————_- 5. 29) 
dat Vk 


The physical and mechanical properties of the polymer are 
related to the number average(M,) and weight average(Mw) 
molecular weights and these averages are related to the 


moments of the dead polymer distributions as: 
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Moa MW (it hy) / Cus + No) (5.30) 


Mw = MW (ue x4 rz) ZL (uy a r+) OS a1) 


where MW is the molecular weight of the monomer. Since A, is 
negligible compared to u,, Equations(5.30) and (5.31) can be 


Simplified to the following equations: 


bis 
" 


MW(XM,) /to €5°.32) 


MW Mi /Ho 


Mw 


" 


MW u2/uy, = MW u2/(XMo) (5.33) 


From -Bquations(5.6),° (5.24), €5.26)and (5.29) one can 
determine the time variation of I, X, uo and uz which are 
initiator concentration, monomer conversion, zeroth and 
second moment of dead polymer respectively. 

be can be appreciated that the original 
Equations(5.13)-(5.18) have been considerably simplified by 
the introduction of the moments and the quasi-steady state 
approximations. These four nonlinear differential equations 
now can describe the solution polymerization of MMA in a 
batch reactor. 

In the development of the above mathematical model it 
has been assumed that the reaction rate constants do not 
vary during the course of isothermal polymerization. However 
experimental results reported on free radical polymerization 
of MMA show the presence of an autoacceleration (gel effect) 


in the polymerization rate as the reaction 
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proceeds(Odian(1981)). This dramatic increase of the 
polymerization rate is due to the increase in viscosity of 
the reaction mixture which reduces the mobility of the live 
macroradicals and therefore the rate of termination 
reaction. It should be noted that at very high conversions 
(above 60%) the propagation reactions may also become 
diffusion controlled, reducing the rate of polymerization. 

It iS important to include the decrease of the 
termination reaction due tofyvthe tgelrveffect ine:the 
mathematical model of the polymerization reactor. Many 
mathematical expressions have been developed for the gel 
effect based on empirical and free volume approaches (Friis 
and Hamielec(1976); O'Driscoll et al.(1979); Ross and 
Laurence(1976), Marten and Hamielec(1979)). 

Recently Schmidt and Ray(1981) extendended the free 
volume correlation of Ross and Laurence(1976) to the 
solution polymerization of MMA to predict the gel effect as 
a function of monomer conversion, temperature and _ solvent 
concentration. The gel effect is calculated by assuming that 
glass transition occurs at ae free volume of 09025 
(Beuche(1962)). 

In this study the gel effect was initially taken into 
consideration after the free volume correlation developed by 
Schmidt and Ray(1981). However the experimental results 
presented in Figure 4.3 show that there is no significant 
gel effect present at the conditions of polymerization. High 


solvent/momomer ratios and lower conversions used in this 
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Study assure negligible gel effect during the course of 


polymerization. 


5.4 Parameter Estimation 

A nonlinear state space model for the solution 
polymerization of MMA ina batch reactor was developed in 
section 5.3. These differential equations can be solved to 
obtain the time variation of I, M, wo, and wz during the 
course of polymerization by specifying the numerical values 
of the kinetic parameters of the model. Using the reported 
values of the kinetic parameters by Schmidt and Ray(1981)for 
solution polymerization of MMA, Significant discrepancy 
between the experimental and model predicted results was 
observed. Since there was considerable scatter of the values 
of the kinetic parameters reported in the literature, it was 
decided to reestimate the parameters of the model using our 
own experimental data. 

The parameters that appear in the set of four nonlinear 
differential equations that represent the reactor model are 
RByapcekpfulk: an kee hiek) pak, Bandsieeincthestollowinguwe will 
discuss how the numerical values of various parameters of 
the model were obtained. 

In MMA polymerization it is generally assumed that 
termination reaction occurs predominantly by 
disproportionation mechanism. However polydispersity 
values(Mw/M, 1.7-2.0) obtained in this work suggest that 


both modes of termination reactions should be taken into 
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account. Odian(1981) reported that at 60°C, termination by 
combination account for about 40 percent of the total 
termination reactions in MMA polymerization. 

When Kee wWaS set equal to kiya, the theoretical 
polydispersity values varied from 1.88 upwards, matching the 
experimental polydispersity values better. If the 
termination by combination mechanism was neglected the 
theoretical polydispersity values would vary from 2.0 
upwards. However, it should be noted that the contribution 
of the disproportionation mechanism to the total termination 
mechanism increases with increase in temperature in MMA 
polymerization (Odian (1981)). This effect of temperature on 
the modes of termination reactions was not taken into 
account in our parameter estimation studies. 

ft bis Aknowninthat @notagald: the initiator radicals 
produced in the dissociation reaction are able to initiate a 
polymer chain because of many Side reactions. The initiator 
efficiency f is defined as the fraction of radicals formed 
in the primary step of initiator dissociation which are 
successful in initiating a polymer chain. The values of f 
for most initiators lie in the range of 02.3=0..8 
(Hamielec(1976)). In this work a value of 0.5 was chosen 
which is also the value used by other workers (Schmidt and 
Ray(1981)). 

It is noted that the rate constants for chain transfer 
reactions to monomer and solvent are very low compared to 


propagation rate constant. The numerical values of the chain 
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constants were taken from Brandrup and Immergut(1975) and 
are reported in Table 5.1. The kinetic parameters which were 
estimated in this work using the experimental data are k,, 
kp cand ky 

But the kinetic terms kp and k, appear in the model 
equations in such a way that their values can not be found 
individually, but only the ratio kp?/k, can be estimated 
from the experimental data. Since the individual kinetic 
constants are reported in the literature, we re-adjusted all 
the parameters simultaneously so that the experimental and 
model results match. Also we checked that the values of the 
estimated parameters are well within the range of values 
reported in the literature. 

Each kinetic constant consists of two parameters 


through the Arrhenius relationship: 


k; = A; exp(-E,/RT) (5.34) 
where 
k; = kinetic constants ky, kp, or k¢ 
A; = Arrhenius factor A, , Ap , or A, 
E; = Activation energy E, , Ep , or E, 


According to Equation (5.34) a total of six parameters 
consisting of three Arrhenius factors and three activation 
energies have to estimated. The above problem reduces to a 
standard problem of estimating numerical values of 
parameters ina system of Bir terencial equations. Numerical 


procedures for solving this problem have been reported in 


‘Wredomesag cet Fe sseitno> ally oni 


So 


“ Lee YP : Pa a 
: oes) bm Ve ' = } a5) 
ear et De) Sh a ot. 
7 a “Ot [= a ee a a) 


bi « (ETEL) ) sup eman bre qu tot 2 7 pe 
ovaw dotdy areemezeT of sone oan fee ofdst at 8 c 
O36 lead Hetnomiseqns od goku bee ates t 
va i) ac ‘ _ 
kebom arith” ay Teeqae” i bn yt ening alienka| 
Ktifet ed 20m Hes am hii aibe2 ted) gew" w dawe ay 
fae ea vas . aw le sitar gy gino dud: » videyboy 
clients Fsokielank ati egniz Pr Lomenee 
lis beragrhe-m aw pPussnreril wad ih. neroney es) Ler 1 
hina imei nT a wit jitie me atipemenbum: a ; TT 
chy 36 waaay wits $na2 OP 91> Bw ante jted om avi 
aol ev Sor-eapoe? ade diets bw ilsw oie ages eae “s | 


eh  aruunanes it 


: " % 
- ‘ 


a idsnotanl f aisinon sk 


eae - ae a atid 
ict) Be } eS ‘ ae 2 
(ee, 8) oh’ | RD ax ne 
—- 2 co cd p oiseninen ote . x 


oa 


of vs BA, a seine) * euinads3A * tr, 
*, 7 Be ial ad noluevisoa = 


- j 


areca’ 
atprtomayag ae to twee a $B ah at ia 5 ; = 
neliqvisaa oe bra euovosa aeinonssh sertd 2opr aiano> 
e Od ae wetdesg veda aft Bassa sao o3 sats 
8): aeltoe bee erage ené samba fo’ wi ed nae 
Lae rani cats sb Dea: o mesays 5 @iesems | 


mi bes sogar" “aie sh woken br gation *° 
} io an io a ae Ee Se 
ae : i 7 


120 


Table 5.1 Numerical Values of k, and k, Parameters 
(Brandrup and Immergut(1975)) 


Temp Coax tos Co x1u* 
mS Cy=ke/kp Cm=Km/kp 
60 . 190 O05 15 
70 Belay) = 
80 7918 .0240 


the literature and are given below. 

Rosenbrock(1966) dealt with the problem of estimating 
parameters in differential equations by minimizing the sum 
of squares of errors between the meaSured values and 
computed values of the states. He also gave an analytical 
procedure to calculate the confidence intervals for the 
parameters by linearizing the model equations. 

Donnelly and Quon(1970) used a quasilinearization 
method to identify parameters in systems of ordinary 
differential equations derived for chemical engineering 
problems. Hosten(1979) carried out a comparative study of 
short cut procedures for parameter estimation for an 
experimental catalytic reaction system. He reported that 
Maximum Likelihood treatment is a laborious task requiring 
up to 100 times the computer time required by short cut 
methods. However since’ short cut estimates do not possess 
desirable statistical properties, it is important to follow 
up any such procedure with rigorous treatment such as 


Maximum Likelihood. Recently Wantenbe and Himmelblau( 1983) 
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examined a number of aspects of estimating coeffecients in 
kinetic models. They reported about a class of methods in 
which the derivatives in the differential equations are 
approximated by a proper difference relation to yield a set 
of algebraic equations that can be solved to obtain values 
for the coefficients. 

In this work a finite difference Levenberg-Marquardt 
method (Becks and Arnold(1977), Bard(1974)) available in the 
IMSL Library (ZXSSQ) has been used in combination with a 
routine (DGEAR, IMSL Library, Gear(1971)) to integrate the 
system of ordinary differential equations for estimating the 
parameters. 

The idea behind the Levenberg-Marquardt (L-M) method 
which combines the Gauss method with the steepest descent 
method is the following. The steepest descent method often 
works well in the initial iterations, but the approach to 
the minimum becomes progressively slower. On the other hand, 
the Gauss method works well and converges’ fast in the 
neighbourhood of the minimum but gives troubles or _ even 
diverges when the initial guesses are away from the minimum. 
This L-M method starts as the steepest descent method and 
progressively switches to the Gauss method when the minimum 
of the objective function is reached thus ensuring 
convergence as well as fastness of the algorithm. 

For the parameter estimation of the reactor model, the 
differential equations representing the batch reactor are 


numerically integrated assuming an mnitval set of 
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parameters. At times t; at which measured data are 
available, the computed values of the outputs Yp; are 
compared with the measured values of the outputs Y,;. An 


error vector, e(,;) at the time t; is formed as follows: 


Clerk bi) (5.35) 


ej 


If ais the vector of the six parameters consisting of 
the Arrhenius factors and activation energies, then a scalar 
function J(a) can be formed by summing the weighted squares 


of the errors as: 


(5.36) 
where w is a 3x3 weighting matrix and n is the total number 
of times meaSurements are made. Note that at each time t,; 
three measurements namely conversion, number and weight 
average molecular weight are made. This is a case of a 
multiresponse problem in contrast to single response problem 
in which only a single measurement is made at each time. The 
function J(a) is minimized by a proper choice of the vector 
a of parameter using the L-M method. 

One of the main difficulties involved in the least 
Squares estimation procedure is how to carry out 
successfully the optimization of J(a) for a highly nonlinear 
model. Because of the complex nature of the model solutions, 


local optima can be expected to be encountered, hence the 
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choice of the initial guesses for the parameters is 
important. Since the numerical values of states in the model 
are of different magnitude, a Gear routine (DGEAR in IMSL 
library) which can handle stiff differential equations with 
the provision for automatic step size and error control was 
used. 

The kinetic parameters reported by Schmidt and 
Ray(1981) were used as the initial values, since it is known 
that kinetic parameters that would be estimated would not 
differ significantly from these values. The data obtained 
from six different experiments were used for the purpose of 
parameter estimation. Table 5.2 gives the details of the 
experimental conditions and the number of data points used 
from each experiment to estimate the parameters. A total of 
159 data points were used. 

Since the initial guess of the parameters was good, the 
program converged fairly quickly in five iterations. These 
parameter estimates along with the values of the kinetic 
constants are listed in Table 5.3. 

Figures 5.1 and 5.2 show the "goodness" of fit of the 
parameters in predicting the conversion-time histories for 
various experimental conditions. The experimental points 
along with the model predictions (solid lines) clearly show 
the agreement between predicted values and the experimental 
values are quite good in all the six experiments. 

Figures 5.3 and 5.4 compare the experimental and 


predicted values of number average molecular weight. It can 


father ote ab eetiete! Bo nau tew Sees 
TEV ni TAG earawes 7880 oy 


43009 Toase ots ‘ania gota 
ack 


bivitos ag | besiege 


sa Shwe ve tant dated, otuoy sents 
thetletoo sope ant .esuhso oaeds ‘nost Atsheuiah 
azogauts ots 368 bemu: wre~ c (osiabaaes pera. | 
wit %o.. ehiezeh one ene de ca aiden pnolyantzaa’® 
+s ayniog wnat So Hodye on\(haa amano gan E 
» £eToe & leis ad a domi xages mE! 


a ate 


wna ‘ean. 7 
ett Song ‘qam csatamiting nid a ng 
path .Rno! hemes ave, oft yfsslup Cae 
ident he sonnel ot fhe grate es | 
me ob sid nh ry eas. 
oda Bo WER ae: *amwtBoo" atid woite 2 Bast. 2 ueruy 
WR) 29 isneagiitit sahithititeiias eid gritsatbexq nite 
as) 9 ‘Loatemizeges edt. enoizibaem a 
vont yizdely qian et bites) anoisotbeag hebom edt fziw - profs 7 
Caiideh tits edt bn, aouilev bessbheq: Abewssd inemesipe‘sdt 
| Teytemizegue cta oft Lhe ab bovp stiup sta geukev 
‘bas foinemeseqye ad »équee £2) Bre eve eee it 
Asa ot .siliptee talvwelom spsteve cedaun Ic estibav ‘beyotbesq 


“es... 


124 


J-UOISIAAUOD IO} sj[Nnsey Jepow pue TeyueuItiedxy Jo uost1redulog ['G ‘sty 


SONUIP{ UT SUIT], 


0'082 0'0%e 0'002 O'O9T 0'0eT 0'08 0 OF 00 


0'0 


S}[Nsey Jepowm —— 
SO’) 9 ‘G4=x 
sO. CU COL =U 
rere) * ‘6¢9=G 
1/1OW 2. 


t) s 


I 6 J 


0 Z°0 
UOISIaAUO) 


9°0 


8°0 


O'T 


=! 


Pre) 
' 
yer 
9 : ‘ 
boom « ie 
% H 
4 vs , 
— , 
i> 
44 
| ci 7 


nr 
aa * 


5 
—_ eae a 


Sli{st 


* 


jase ro Po Apee 


Vibciin 


91 Cousbuton 01 F 


a7, 


ye 


a8) 


[I-UOISIA@AUO) IO] S}[Nsay [epoy{ pue TeyUsawItIedxy Jo uostiIedulog 2G SIy 
SONUIP, UT SUIT], 


0'082 0'0%2 0'002 O'O9T 0'OeT 0°08 0O'OV 0'0 


0°0 


S}[Nsey [Jeapo_y —— 
OF 9 SiZ=~x 
OT 9 ‘OL=V 
OT ' “¢9=0 
1/1OW 2. 


1 @ XL 


v0 20 


UOISISAUOZ) 


a0 


8°0 


OT 


por 


fOURSLSICI-Ty 


*~ 
Lb 


rd 
= 
~ 


hie as coushsii2on 0] pebci wens} Si Wo 


qay Bens 


“7 
4 


yw mo yUnped 


$a0'C 


Svo'0 


sooo 


<> 


5800 


wy 


vi ue : : 


7 


126 


1-319 M Te[NoeTOW eFe19ay JoqUINN Jepow pure [eyUaUItiadxg yo uostieduroy gg ‘Sty 


SOYNUTIP UT SUIT], 
0'082 0'0F2 0'002 0'09T 0'02T 0'08 0'0F 0'0 


s}[Nsay Japow —— 


sO’ 9° ‘SL=x 
SO 9 ‘O4=V 
so 9 ‘S9=0O 
1/1OW 0; 


I ¢@ J 


0 


0? 0°0 


O'F 


1USIOM TOW eseseay JequIny 


08 0°9 


omen! 
Ole 


edivest labo ~—— 3 a. = *, 


0.088 AGS 6005 e03f COS: 208 a ee 
roe geiucit atemiT | + 


l-tgis¥ 1chyoeioM ageieyaé odie Isbolt ban isinedtivsqus to Beehagme? €a grt 


y1-}U5TOM JE[Nofop eseIsaAY IJaequINN [epow pue [e}yUaeuIt1edxg so uostiedwiog FG ‘Sty 


SOPNUTP, UT SUIT, 


0'082 0'0rve 0'00¢ O'O9T 0'OeT 0'08 0'OV 0'0 


00 


sy[nsey Tepow 


Or’ G SL=x 
or F ‘0L=¥ 
OF 9 ‘S$9=0o 
1/TOW D; 


0 $s 


I 6 J 


08 09 oF oz 
1431IOM TOW esetaaAy Jaquiny 


oot 
nORE: 


/ i * 


T 
= a ee ap=c a ee a ere SS a> ae oe 
+-€tt. =o 


etigee® ilebok -——— 


0.085 ods  a00S oo Gos 008 ‘om’ ; 0 
aeturiM mi seit | as 


l- tigie® 1slwssioM egatsvé tedme lshoM bas [sissoheged to geetisg 09 BE mt 


5 
4 


ae eee 


128 


J-7USIOM IeTNIe[oW eses1aay YUSIAaW [epoy pue [eyuUauItiadxg so uostiedurog Gc ‘BIy 


SONUI{ Ul SUIT], 


0'082 0'0ve 0002 0'O9T 0'0eT 0:08 0'OF 00 


0°0 


Szy[Nsay [Te pon 
SOF (95. °“Sic= x 
co; 6s OL V 
sO 9 S9=a 
1/1OW Je 


i ae Noa 2 


O-9T oI 08 OF 
JYSIOM TOW edetoay yUSIOM 


0°02 
lx 


oc6s cos 000 O08: UeSI ¢98 (oor; . 60 
estima ai soul ; 


l-JdgioW taluosloM sgsievA tdgio¥ isholt Sac isinemiisqud to Hogitagaics €£ att 


129 


TI-}US19 My TeINoeToW aseraay FIO Tepow pue [eyUaewtiadxg so uostieduro) gc ‘B1g 


SOYNUIP UT UIT], 
o'082 0'0F 0002 0'09T 0'02T 008 0'OF 00 


S}[Nsay Topow —___ 
OT 9 ‘SL=x 
OT 9 ‘OL=V 
ot 9 ‘“$9=0 
1/1OW Je 


0'8 


O'eT 
WUSIOM TOW eserseay pYSIOM 


0-0 


O'V 


O'9T 


0°02 
,O lx 


est 


‘>it 


~~ 
~ 


See g 


porter yastats Wopscayst yet 


¥ 
s 


gst pg combenieos oy Exbexmmenss) VIG HOSS 
978 


390°C 


ite 1p apres 
7809 e 


130 


Table 5.2 Experimental Data used for Parameter Estimation 


Exp wig solvent Io X M, Mw Total 
no 2G Ratio mol/L Data 
1 65 .6 2-05 9 9 9 ZT; 
2 65 oo = 0 9 9 9 27. 
3 70 36 OS 9 9 9 Ze 
4 70 «6 ako 9 9 9 2 
5 TS 0 Qo 9 9 9 vg 
6 TAS ab = 1,0 8 8 8 24 


Total number of data points =159 


1 0 0 
w= 0 1/50000 ) 
0 0 1/80000 
be 


Table 5.3 Estimated Values of the Kinetic Parameters 


Age @ to Uk 10. S Es = 33740 

Ap = 1.201x10’ Ep = 9630 

A, = 2.113x10° ee = 1018s 
Typical values of kz, kp and k, at different 
Temperatures 
4g koe, 1 O:* kp kx oe * 
.C 1/s L/mol-s L/mol-s 
65 .0490 68a, 4.648 
70 O22 840.6 4.751 
75 «2085 1030.0 4.854 


be observed that though the model predicts the number 


average molecular weight well, some of the experimental 
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points do differ from the model results. The accuracy of 
number average molecular weight measurements from _ size 
exclusion chromatography has been discussed earlier in 
Section 3.4. A meaSurement error in the range of 10 percent 
can be expected. 

Figures 5.5 and 5.6 show the comparison of experimental 
and predicted weight average molecular weights for various 
experimental conditions. A similar argument that has been 
advanced with respect to number average molecular weight 
holds good since weight average molecular weight has also 
been measured using SEC. 

It should be mentioned here that the model prediction 
of conversion, number average and weight average molecular 
weights considerably improved after the parameter 
estimation. The values of kinetic constants taken from 
literature did not predict the experimental values well. 
Therefore the newly estimated parameter values are used in 


the derivation of optimal control policies in chapter six. 


5.5 Conclusions 

In this chapter a mathematical model has been derived 
for the free radical solution polymerization of MMA in a 
batch reactor. By using the moments of the polymer chains 
and the quasi steady state approximation (QSSA), the order 
of the model was reduced significantly. 

The model parameters were estimated from our 


experimental data by using a nonlinear parameter estimation 
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routine (ZSSQ IMSL Library) due to Marquardt(1963). It is 
found that when the estimated parameters are used, the model 
predicts experimental values of conversion, number and 


weight average molecular weights reasonably well. 
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6. OPEN LOOP CONTROL STUDIES 


6.1 Introduction 

In the polymerization industry there is considerable 
economic incentive to develop real time computer control 
policies that will increase production and produce a polymer 
with desired final physical and mechanical properties. As 
discussed in chapter five a homopolymer is mainly 
characterized by the molecular weight distribution and its 
mean values e.g. number average (M,) and weight average (Mw) 
molecular weights. These two variables greatly affect the 
final chemical, physical and mechanical properties of a 
polymer product. Therefore if we want to produce a polymer 
with specified final properties, we must be able to exercise 
control of these variables (M, and Mw) during the course of 
polymerization. In the past most of the studies on the 
optimal control of polymerization reactors have been limited 
to theoretical predictions of the optimal control policies 
such aS optimal temperature and optimal initiator addition 
policies to produce a polymer with desired final conversion 
and desired molecular weight averages. Few if any studies 
have reported actual implementation of such optimal policies 
to laboratory scale polymerization reactors. Hence there is 
an obvious need to demonstrate and evaluate computer 
implementation of these optimal control policies to a 


laboratory scale reactor. 
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In this chapter we develop optimal control policies by 
applying state-space based control methods to the model of 
the batch polymerization reactor derived in chapter five. 
These policies are implemented on the experimental batch 
reactor system and the results compared with the theoretical 
predictions. However in open loop control of the reactor no 
feedback corrective action is taken should there be any 


process upsets or disturbances entering the system. 


6.2 Control of Polymerization Reactors 

Compared to the successful application of optimal 
control theory to other chemical processes (distillation, 
evaporator etc.) there have been only a few successful 
attempts at the optimal control of polymerization reactors 

(Kiparissides (1978)). Several valid reasons exist for this 

lack of progress of application of the optimal control 

theory to polymerization reactor systems: 

1. Polymerization reactions are highly complex and at least 
until relatively recently it was difficult to describe 
such systems by low order models. 

2. The formulation of a mathematical objective function to 
optimize the performance of the polymerization reactors 
is not easy. 

3. On-line measurements of important polymerization process 
variables such aS conversion and molecular weight 
distribution until fairly recently has not been 


available. 


vi asixzifoq lerinos Samiyqo g 
ic Ishom eA? oF Bhoddom consis aaa soma 
svit tesqedp ae aiadtlaie 10} 980% WoL sexi se 
fozad leineniiseqaem a3 “90 betnemalgai e216 mens — 
ispéiseiced? ofe 423m Bexaqnop apiipet an tne aaa 
sti sor odes add Be salimiind qool. rade, ad peer 
vin oa e2edt Biwede) nade. ef Monae awitoay ! 
ieteve ed gphiwearne asonad td Sai> 16 @: 


j 
oo 


avesiceA aodteatzemylo4 to 1 > 
ehizge to neigeot igus ~tyteeenpue | it on eT GD 
ellsaeib) eseeenenq [soiesda, * merigo a yioady + 
Ltasessva wal © 46h nmed» dived) meds” (aaa 
@toteses. notiss Lapmpheg (ao Lorsriog Lamigo ld ae 
2id2 wo dalew acomeag bt ine: fetevee i ARES) par 
i2ao>  famidge eae jo noidea! lqge ae eaangord \Me 
sepesmye soazannneiboastsomytog of 93 
reaol 3c Sas xeiques gidpid gis enoiiinsss nots smyiod 
sdiicesb O31 Jiuae ERG gee 41 ghonegee ‘eviaieaeen 
- i) 1 8 6 = «@iebom aehae Ol yd-emeteye’ di 5 | 
os ooltaqu? evitoette Ieoizancitom @ Se fottslumied ent) ra 2 
220;%es" a@liesivamyfeq edz 30 samemegotssq oft. v6 | 
7 ase ton. at! * 7 | 
easco'rg woltagitesrleg tnai:0ogqml, io 2inemewesem snik-n0 we 
sAyieyw seluseiom ban Hoiersvneo as aoe —eeldsiret 
need 66 «6©@mt 6yiveese: «6©6 Nisin? ©«6 bien nobtudiade tb 


-sidaliave 


6.2.1 Previous Work 

In the past several investigators have applied optimal 
control theory to polymerization reactors to determine the 
temperature and initiator addition policies that minimize 
the reaction time and produce a polymer with desired final 
molecular weight distribution. Hoffman et al. (1964) 
developed policies for monomer or initiator addition to a 
batch polymerization reactor to produce a polymer product 
with desired final conversion and number average molecular 
weight and at the same time having minimum polydispersity. 
They concluded that in polymerization reactor systems in 
which gel effect is absent, the narrowest molecular weight 
distribution of the polymer is obtained when the 
instantaneous degree of polymerization is maintained 
constant throughout the course of polymerization. Beste and 
Hall(1966) extended this study of molecular weight 
distribution control to the batch polymerization reactors 
and determined addition policies of initiator and monomer. 

In one of the earliest applications of optimal control 
theory to polymerization reactors Shatkan(1965) determined 
the temperature histories to obtain a polymer with desired 
final conversion in minimum time for the case of an ideal 
batch polymerization reactor. According to this policy, the 
temperature must be raised at a predetermined rate until it 
reaches an upper limit and maintained at this level until 


the desired target set is reached. 
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Ray (1967) calculated the optimal temperature and 
initiator concentration policies fOr a series of 
polymerization reactors by minimizing a performance index 
which 1s a function of deviations of conversion, number 
average degree of polymerization and polydispersity from 
their desired values. He applied a peak-seeking method to 
compute the best temperature and initiator concentration in 
each tank. Hicks et al.(1969) formulated a similar objective 
function and applied Pontryagin's Maximum Principle to 
optimize the performance of a free radical polymerization 
reactor: (i) by controlling the polymerization temperature 
while keeping a constant initiator concentration, and (ii) 
by using both temperature and initiator concentration as 
control variables. The latter attempt proved unsuccessful 
due to computational difficulties. 

Yoshimoto et al.(1968) applied the Maximum Principle to 
a free redical solution polymerization reactor to obtain the 
desired average degree of polymerization and conversion in 
minimum time. Their results are similar to those of 
Shatkan(1965). In a later study Yoshimoto et al.(1971) 
determined optimal temperature policies to obtain specified 
conversion and average molecular weight in minimum time for 
the case of a thermally initiated polymerization reactor. 
They demonstrated that there exists a unique initial 
temperature and thermal history that permits to reach the 
selected final conversion and average molecular weight in 


Minimum time. However their model did not take account of 
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the gel effect caused by the diffusion controlled 
termination reactions. 

Nishimura et al.(1968) also applied the Maximum 
Principle to an ideal radical polymerization reactor to 
obtain the optimal temperature and initiator concentration 
policies. They implemented these policies to a_= solution 
polymerization of styrene in a batch reactor and found good 
agreement between the experimental results and theoretical 
predictions. 

Osakada and Fan(1970) calculated sub-optimal 
temperature and initiator feed rate policies in an attempt 
to obtain a desired molecular weight distribution. The near 
optimal policies were represented by two polynomial 
functions of time the coefficients of which were found by a 
simplex pattern search technique (Nelder and Mead(1969)). 

Sacks et al.(1972) applied the Maximum Principle to 
determine temperature and initiator policies that minimize 
the time to obtain a final desired conversion and number 
average molecular weight. They showed that the optimal 
initiator policy is one of maintaining constant initiator 
concentration throughout the course of polymerization. The 
Same authors (Sacks et al.(1973)) studied the effect of 
temperature variations on the molecular weight distribution 
of the polymer. By using the Maximum Principle they found 
that the temperature policy which minimizes the breadth of 
the molecular weight distribution should produce a polymer 


with constant instantaneous number average molecular weight 
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throughout the course of polymerization. The temperature 
policy which maximizes the breadth of the molecular weight 
distribution is one of requiring step changes in temperature 
between prespecified maximum and minimum values. 

Chen et al.(1978,1980,1981) conducted a series of 
studies on optimal control problems related to free radical 
polymerization of styrene ina batch reactor by considering 
the polymerization temperature and the initiator 
concentration as two control variables. They derived optimal 
temperature and initiator policies that are required to 
obtain a polymer with desired final conversion and number 
average molecular weight in minimum time. Experimental 
verifi¢ation |\"Yof } thein: spolieress ifor ithe icasewof solution 
polymerization of styrene ina batch reactor revealed that 
there was good agreement between experimental and calculated 
conversion values throughout the course of the 
polymerization reaction. However appreciable deviations 
existed between experimental and calculated number average 
and weight average chain lengths. 

Masterson(1977) and Clough et al.(1978) applied the 
Maximum Principle to solve the minimum time problem for 
batch polymerization of styrene. They obtained a_ product 
with desired conversion, desired number average and weight 
average molecular weight by controlling temperature and 
initiator feedrate. The numerical solution of the above 
problem has proven challenging. They used a target set 


concept in the sense that they searched for final values of 
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conversion, number and weight average molecular weights in 
the "target set" which satisfied the initial conditions of 
the state variables. 

Wu et al.(1980) presented a graphical solution to _ the 
minimum time problem of styrene polymerization. Their 
performance index depended only on conversion and number 
average molecular weight at the end of reaction. Applying 
the Maximum Principle, they were able to obtain the the 
optimal variation of temperature with time. Their 
theoretical predictions were in agreement with experimental 
measurements of conversion and number and weight average 
molecular weights. 

Farber(1983) applied optimal control theory to bulk 
polymerization of styrene in a batch reactor to obtain a 
polymer with desired conversion and molecular weight 
averages in minimum time using temperature as the control 
variable. Recently Thomas (1981) and Thomas and 
Kiparissides(1984) applied the optimal control theory to a 
batch polymerization reactor for MMA to calculate the 
near-optimal temperature and initiator policies that are 
required to produce a polymer with a desired final 
conversion and desired number average and weight average 
molecular weights. They numerically solved the resulting 
two-point boundary problem to obtain the optimal control 
policies by introducing a "target set concept" which 


Significantly simplified the numerical solution. 
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6.3 Statement of the Control Problem 

From the previous review of literature on the optimal 
control of batch polymerization reactors it becomes evident 
that most of the studies were directed to the derivation of 
theoretical optimal control policies using the model of the 
reactor and only very few implemented these policies to an 
actual reactor in an effort to verify experimentally the 
theoretical Pesutcrec. Such tamiverification is not only 
important but also necessary to successfully close the wide 
gap existing between the theory and practice in chemical 
process control (Foss and Denn(1976), Athans(1976)). 

In the following sections optimal control policies for 
the solution polymerization of MMA ina batch reactor will 
be developed using the mathematical model of the reactor 
derived in chapter five. 

A variety of optimization problems can be formulated 
for the batch solution polymerization reactor depending upon 
the objective function that has to be minimized and the 
selection of control variables. The reaction temperature and 
the initiator concentration are usually selected as the 
control variables. These two variables affect the rate of 
polymerization and the molecular weight distribution of 
polymer as seen in chapter four. 

The following problems have been formulated and solved 
for the free radical solution polymerization of MMA in a 
batch reactor. 


Problem 1: Find the optimal isothermal temperature and_ the 
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optimal initial initiator concentration to produce a polymer 
with desired conversion and desired number average molecular 


weight in minimum time. 


Problem 2: Find the time optimal initiator addition policy 
for a given polymerization temperature to produce a polymer 
with desired conversion and desired number average molecular 


weight in minimum time. 


Problem 3: Find the optimal temperature policy for a given 
initial initiator concentration to produce a polymer product 
with desired final conversion and desired number average 


molecular weight in minimum time. 


Problem 4: Find the optimal temperature policy for a given 
initial initiator concentration to produce a polymer product 
with desired final conversion and number average molecular 
weight and with minimum polydispersity. 

The first problem reduces to a static optimization 
problem. Problems two and three require a minimum time 
formulation based on the Pontryagin's Maximum Principle. The 
fourth, minimum polydispersity problem, turns out to be a 


free terminal time problem. 
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6.4 Optimal Isothermal Temperature and Optimal Initial 

Initiator Concentration Policy (Problem 1) 

The objective of this static optimization problem is to 
find the optimal isothermal temperature and the optimal 
initial initiator concentration to produce a polymer with 
desired final conversion and desired number average 
molecular weight in minimum time. A solution to this problem 
is obtained by using the Lagrange multiplier technique. 

The mathematical model of the free radical solution 
polymerization of MMA in a batch reactor has been derived in 
chapter five. The model equations were further simplified as 
shown in Appendix C. 


The model equations are 


Gb/dt..=34kyl (6.1) 

aM/dt = -k,M/I Cee2) 

duo /dt = Ckal + k-M/1I +7k;Sy71 Com) 

Gis 7dt =e CkYS + kaM)As tue hone + kee AG? (6.4) 
where 

No = Zeke. (6.5) 


2Eky TP CKDM + uk,M +k. S)iXe 
A, = —_-_O (6.6) 
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k, = kpy2fka/k, = A,exp(-E,/RT) (6.8) 
k 2 = K mV 2Ek a/R, = A2,exp(-E2/RT) (6.9) 
k3 = Ki yaZekecks = A3;exp(-E3/RT) C6510) 
C = 2£ 001 = Ke? Zi ae) C6214) 


Equations(6.1) through (6.11) completely describe the 
solution polymerization of MMA ina batch reactor. For a 
constant temperature Equation(6.1) can be analytically 


solved to give 


I = Io exp(-kyt) CGfk2o 
where Ip5 is the initiator concentration.Substituting 
Equation(6.12) into Equation(6.2) and integrating the 


resulting equation we get 


kK; kgt 


M = Mo expl VYIo 2 { exp(- ——) -1}] (6e%3) 
2 


d 


where Mo is the initial monomer concentration. Equation(6.3) 


can be integrated as: 


t 
fo = [Ckal at + fk >M/I dt + fk3SVI dt (6.14) 


From Equation(6.12) and (6.13), the integrals in the 


Equation(6.14) are evaluated as: 
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[Shee Ree Ciel i-expe@ ket) | Pez 


Sy tol l—-exp(-k.t/2)] 


Kg 


ke ky 
Mo{il -exp{2 


1 d 


VIo Lexp(-kgt/2) -1)}} (6.15) 


The desired conversion X* and desired number average 
molecular weight M,* can be expressed in terms of desired 
monomer concentration (M*) and desired zeroth moment of the 


polymer(uo*) as: 


M* =5Me 61-23%) (6.16) 
Lo™ = MoX* MW/M, * (6.17) 
where MW is the molecular weight of the monomer. 


Equation(6.13) can be solved for the final time as 


2 kaln (a) 
ef Inf 1 + ————— } (6.18) 
kg 2kiV1o 
where 
a = M* vA Mo (6.19) 
Substituting Equation(6.18) into  Equation(6.15) and 


rearranging we get 
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{ k gin ¢a)c},? kK 3 
Ho” - CI,o [1-{1+ —— }] + Sln(a) 
{ 2kiV1Io } k ; 
K 2 
- My (1-a)=0 (6.20) 
ky 


Now the optimization problem can be formulated as follows: 


J= min, tf subjected to Equation(6.20) C6221) 


An objective function J(Io,T) can be formed incorporating 


the constraint Equation(6.20) using a Lagrange multiplier. 


2 kgln(a) 
he hy) = In(1+ ————_—__—_ 
kg 2kPl 6 
a olen K 3 
+y{uo* -erett-G+ My (1-a)+ Sln(a)} 
2k,V1o6 k, C622) 


where y is the Lagrange multiplier. From Equation(6.22) the 


optimal conditions are 
dJ/AIo = 0 (6.23) 
dJ/oT = 0 (6.24) 


From Equations(6.23) and (6.24), Equation(6.25) can be 


derived after considerable algebraic manipulation. 
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-2k ,IoBy-Ck2Sln (a) (E;-E2)+Ck3M,(1-a) (By-E3)+kgVIolna(E,-E,) 
2 hipioantekeele lita) 
~ka(E,-Eqg)V¥1Ioln(a) —————————____ 

2k,1o 


2K abo tk avi ein (a) 


" 


sfeskaVioln(a) bint 0 (6.25) 


2kiVIo 

The optimal values of initial initiator concentration 
ie and isothermal temperature T were obtained by 
Simultaneous solution of Equations (6.25) and (6.20) for 
which the Newton-Raphson method was used. 

Table 6.1 gives the optimal values of I, and T that 
must be used to obtain a polymer product with the desired 
conversion and number average molecular weight in minimum 
time. 

From the results in Table 6.1 it can be observed that 
the combined policy leads to almost total depletion of 
initiator, the initiator conversion being 99 percent in all 
the cases. This situation in which the initiator is fully 
consumed before all of monomer is converted to polymer is 
called dead-end polymerization. Chen and Cheng(1978) also 
obtained similar results in batch bulk polymerization of 
styrene. The isothermal temperature obtained are high in the 
range of 103=1232C whereas the initial initiator 
concentration are very low in the range of 0.006-0.004 
mol/l. It can also be observed that when desired number 


average molecular weight is increased keeping desired final 
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Table 6.1 Simulation Results of Optimal Isothermal 
and Initial Initiator Concentration Policy 


xs M.* he Tr cé Initiator 
mol/l ae min Conversion 
-70 50000 .0066 t2s o.5 ~99 
a7 0 FOO0Ohs 0047 117 19.7 .99 
Be 8 90000 .0036 112 34.0 .99 
.80 50000 .0075 114 25°80 ~99 
- 80 70000 .0054 108 52.0 .99 
-80 90000 .0041 103 90.0 99 


conversion constant, the isothermal temperature as well as 
the initial initiator concentration decreases. However total 
time to produce higher molecular weight polymer increases as 
expected. 

These optimal temperature and initial initiator 
concentration policies reported in Table 6.1 were not 
implemented on the experimental reactor system because of 
difficulties of maintaining high temperatures inside the 
reactor. 

These results should be interpreted with caution. 
Probably the model itself is not valid at these high 
temperatures and hence one may question the accuracy of the 
Quantitative results. However the qualitative predictions 
obtained from the model should be valid. Since the rate of 
reaction is high at high temperatures, considerable amount 
of heat will be released by the exothermic MMA 
polymerization. The reactor should have enough cooling 


capacity to carry out the reaction at high temperatures. 
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6.5 Time Optimal Initiator Concentration Policy(Problem 2) 

The objective in this problem is to find an optimal 
initiator concentration policy for a given temperature of 
the reaction to produce a polymer with desired final 
conversion and final number average molecular weight in 
minimum time. This problem reduces to the standard minimum 
time problem in optimal control theory for which standard 
procedures of solution are available. 

It should be pointed out here that only desired 
conversion and desired number average molecular weight are 
specified whereas the weight average molecular weight is 
allowed to be free from any constraint. In polymerization 
systems in which the gel effect is absent, as is the case 
here, the polydispersity values do not change very much 
suggesting that the weight average molecular weight is an 
almost constant multiple of M,. Hence the number average and 
weight average molecular weights can not be changed 
independently when the gel effect is absent. 

Now let us apply optimal control theory to the model of 
the reactor given by Equations (6.1)-(6.11). 


For the isothermal process the objective is 


tt 
minCte) = “sde (6.26) 
0 


The Hamiltonian for the above problem is given as: 


H = -1 -P,k,M/I + P2{CkgI+k2MVYI+k3Sy1I} (6227) 
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P, and P22 are co-state variables which must satisfy the 


following canonical equations. 
dP,/dt = -dH/dM (6.25) 


aP,/dt -dH/duy5 = 0 (6.29) 


From the Equation(6.29) it is clear that Pz is a constant. 
If the control variable is unconstrained the conditions for 


optimality are given by: 
dH/aI = 0 (6.230) 
H = 0 (6:31) 


From Equations(6.30) and (6.31) I can be obtained as 


T = =1/(€kGP> ) = "constant (6.32) 


Here I is constant because Pz is constant and ky is constant 
for a given temperature. 
Integrating Equation(6.3) after substituting for dt 


from Equation(6.2) we get 


k 3 Mo K2 
Uo = —— Sln(——) + —M—— (Mo-M) + CI kat (6.33) 
ke M ky 


From Equation(6.2) 


In(M/M,) = In(1-X) = -k,yVIt (6.34) 
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From Equation(6.34) I is obtained as 


{1n(1-x) }? 


Sms ees (6.35) 
Chita 


From Equations(6.35), (6.17) and (6.30) tf can be obtained 


as 
(1-v)k,{1n(1-x*) } 
ea Bees (6.36) 
kz Mo REMW/Mati be eSin( isk" jo- vCaMgx* 
where 
Core lke kD Cyr a-kesko paki (6.37) 


Subsiti tut sng tebaquat ton (6636) ovink Equation(6.35)"s1 acan ibe 


obtained as 


kp? {MoX*MW/M,* + C,S1n(1-X*) -CmMoX* }? 
t= — 6735) 
( Lap) Mb2t keokic { Eni KE 

Given X* and M,* the optimal initiator concentration 
which is constant for a given isothermal temperature can be 
calculated from Equation (6.38). The minimum time required 
is obtained from Equation(6.36). Since the initiator 
dissociates during the course of polymerization, fresh 
initiator has to be continuously added to maintain I at 
constant level. Deriving a simple mass balance for the 


initiator, the initiator feedrate F is obtained as 
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F=Ikyg=-1/CP2= constant (6.39) 
The above equation indicates that the optimal initiator 
feedrate is also constant. 

However in conventional polymerization reactions, the 
initiator decomposition rate is so slow that the initiator 
concentration remains constant throughout the duration of 
polymerization. Therefore the initiator policy turns out to 
be single charge of initiator at the beginning of the 
reaction. But in dead end polymerization reactions where the 
initiator dissociates at an appreciable rate, fresh 
met iator has to be added during the course of 
polymerization to keep the initiator concentration constant. 
It has already been mentioned that dead end polymerizations 
occur in systems that employ high temperatures and low 
initiator concentrations. 

Table 6.2 gives the optimal initiator concentration 
values that were obtained to produce a polymer with desired 
X* and M,* in minimum time for the given isothermal 
temperature. 

The results in Table 6.2 shows that at constant X* and 
Mag when the temperature is increased, the optimal 
initiator concentration as well as the total time of the 
reaction decreases. This suggests that to have appreciable 
time savings the isothermal temperature should be chosen as 
the maximum allowed temperature of the reactor and the 
optimal initiator concentration shoud be employed to. get 


polymer with desired X* and M,”. 
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Since the optimal initiator concentrations are almost 
constant for the temperatures employed in our’ reactor 
system, separate experiments were not carried out to verify 
the results shown in Table 6.2. However some of the 
isothermal experimental results reported in chapter four 
correspond to the conditions shown in Table 6.2. Table 6.3 
compares the experimental results with the results shown in 
Table 6.2. The above results clearly show that there is very 
good agreement between the experimental results and the 


theoretical predictions. 


6.6 Minimum Time Temperature Policy (Problem 3) 

In this problem the objective is to find the optimal 
temperature policy at a given indtiad rntepator 
concentration to produce a polymer with desired conversion 
and number average molecular weight in minimum time. Sacks 
et al.(1973) solved the above problem for the bulk free 
radical polymerization systems and deduced that the optimal 
temperature policy should keep the rate of production of 
primary radicals constant. Since the concentration of the 
initiator decreases with time, the temperature should be 
increased to keep the initiator dissociation rate constant. 
However in conventional polymerizations in which the rate of 
dissociation of the initiator does not decrease 
Significantly, the optimal temperature policy turns out to 
be one of constant temperature throughout the polymerization 


process. 
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Table 6.2 Optimal Initiator Concentration Policy 


Temp x* Me” i tt 
ne mol/l min 
65 50 50000 - 1040 161 
70 ~ 50 50000 SAO eas 109 
75 250 50000 -0520 74 


Table 6.3 Comparison of Experimental and Optimization 
Results: 
Optimal Initiator Concentration Policies. 


Time iit I xan x Mee My, 
min 2¢ mol/L 
161 65 5 10 510 -46 5.0E4 4.8E4 
74 75 -05 250 .49 5.084 4,.8E4 


The mathematical derivation and numerical techniques 
necessary to calculate the optimal temperature profile for 
the given problem are discussed in detail in Appendices C 
and D. 

In Appendix C, the Pontryagin's Maximum Principle is 
applied to the mathematical model of the batch reactor in 
order to minimize an objective function representing final 
time and weighted sums of squares of deviations of M and uy 
from their desired values at the final time. The objective 


function is given as: 
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min J= tft w,{M(tf£)-M*}? + wofuoltf£) - uo}? (6.40) 


In the above equation w, and wz are weights. A two point 
boundary value (TPBV) problem involving the state and 
co-state variables is obtained as a result of application of 
optimal control theory. 

There are a number of numerical methods available 
(Rirk (1970) ) to ‘esolves*this,: difficult: TPBV problem. «A 
gradient method (Thomas(1981)) is used to obtain numerically 
the optimal temperature profile. The minimum time problem is 
solved as a series of fixed terminal time problems. For this 
purpose a fixed time and an initial temperature profile is 
assumed at first. The temperature profile 1S updated by 
solving the TPBV problem using the gradient technique. The 
minimum time is found by solving the TPBV problem for 
different final times and selecting the smallest one for 
which convergence is possible. 

This method to find the optimal temperature profile has 
many disadvantages. The initial guess to the optimal 
temperature profile has to be a good one, otherwise serious 
convergence problems could occur. The step changes that are 
made in the steepest descent direction has to be properly 
controlled. Small steps in the steepest descent direction 
result in too much computer’ time before the solution is 
reached whereas large steps could lead to oscillations and 


divergence from the solution. 
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In Appendix D the optimal temperature profile for the 
given problem is developed by using a simplified model of 
the batch reactor, in which transfer reactions to monomer 
and solvent are neglected. However it should be pointed out 
that the transfer reactions are neglected only to derive the 
optimal control policy; they are included in the state 
equations for model description during the simulation runs. 
From earlier simulation studies it has been established that 
neglecting transfer reactions does not affect conversion 
results whereas molecular weight results do show changes of 
approximately 10 percent. 

As shown in Appendix D the simplified model facilitates 
the direct solution of the minimum time problem of the batch 
reactor. Solution of the optimal control problem is obtained 
by iterating upon the initial value of one co-state variable 
P2. The TPBV problem is thus avoided. 

Computer programs have been developed for both of the 
methods to find the optimal temperature policies. Since both 
methods gave almost similar temperature profiles, the method 
obtained from the simplified model was used for all _ the 
problems because of its computational efficiency. 

In Figure 6.1 the optimal temperature profile obtained 
is plotted for different desired number average molecular 
weights. It can be observed that the temperature profiles 
are almost isothermal which is an expected result for 
conventional polymerization systems. It has been explained 


that the optimal temperature policy should be an increasing 
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policy in order to keep the rate of dissociation of 
initiator constant during the course of polymerization. The 
increasing trend of the temperature policy can be clearly 
observed in curve no.3. From Figure 6.1 it can also be 
observed that when the desired M, 18S decreased, the 
temperature of polymerization increases whereas the total 
time of the reaction decreases. However as pointed out by 
Chen and Cheng(1978) these optimal temperature policies do 
not achieve considerable time saving compared to isothermal 
policies. Figure 6.2 shows the conversion profiles that are 
obtained when optimal policies are used. 

Figures 6.3 and 6.4 show the number average and weight 
average molecular weights development, respectively, in the 
polymerization process when the optimal temperature profiles 
are used. It can be observed that for this policy, high 
molecular weight polymers are produced at the beginning and 
progressively lower molecular weight polymers are produced 
such that at the end of reaction, the desired M, value 
polymer is obtained. Obviously this leads to broadening of 
the molecular weight distribution. 

Two of the temperature policies shown in Figure 6.1, 
corresponding to desired M, values of 50000 and 60000, have 
been implemented on the experimental reactor system. The 


experimental results are documented in section 6.8. 
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6.7 Minimum Polydispersity Policy (Problem 4) 

The aim in this problem is to develop the optimal 
temperature policy for a given ia Cc Lad miitrator 
concentration to produce a polymer product with minimum 
polydispersity and with desired final conversion and number 
average molecular weight. 

Hoffman et al.(1964) derived initiator and monomer 
addition policies to produce a polymer product with minimum 
breadth of molecular weight distribution. He reasoned that 
to have the minimum polydispersity, the instantaneous number 
average molecular weight of the polymer produced should be 
constant throughout the course of polymerization. Sacks et 
al.(1973) later proved mathematically Hoffman's (1964) 
reasoning. They also showed that in conventional 
polymerization the optimal temperature policy is one of 
decreasing temperature during the course of polymerization. 
It can be readily observed (Figure 6.3) that at a constant 
temperature, the number average molecular weight produced in 
the polymerization process decreases with time. In order to 
keep the instantaneous number average molecular weight 
constant, the temperature of the reaction must be decreased 
during the couse of polymerization. 

Optimal control theory has been applied to the 
simplified model of the batch reactor in which transfer 
reactions have been neglected. The detailed mathematical 
development along with the numerical procedure is given in 


Appendix E. As mentioned in Section 6.6, through the use of 
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the simplified model, the computational procedure to obtain 
the optimal temperature policy is considerably simplified by 
avoiding the TPBV problem. 

In Figure 6.5 the optimal temperature policies to 
obtain desired final conversion (X*=0.5) and desired number 
average molecular weight are shown. The optimal temperature 
decreases with time during the course of polymerization, a 
result that is expected from the earlier discussion in this 
section. When higher molecular weight polymer is desired, 
the temperature decreases and the total reaction time to 
obtain the desired conversion is considerably increased. 
Figure 6.6 shows the conversion profiles for the optimal 
temperatures used in the reaction. 

In Figures 6.7 and 6.8 the development of number 
average and weight average molecular weights during the 
polymerization process are shown. It can be observed that 
for the minimum polydispersity temperature policy, constant 
instantaneouS number average and weight average molecular 
weights are produced throughout the course of 
polymerization. Figure 6.9 shows polydispersity values that 
are constant throughout the polymerization. It should be 
noted that polydispersity will increases in the case of 
isothermal polymerization. 

Two of the policies were selected and implemented on 
the reactor system in order to compare the experimental 
results with theoretical prediction. These results are 


discussed in the next section. 
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6.8 Experimental Verification of Optimal Control policies 

In the previous sections, theoretical optimal control 
policies were obtained for a variety of problems in the 
batch polymerization reactor using the mathematical model. 
However it iS important to verify experimentally these 
policies to demonstrate and evalulate on-line experimental 
open loop control of the polymerization reactor. 

For this purpose the open loop temperature policies, 
which are stored in a disc file, are read by the computer at 
the time of implementation and sent aS remote set point 
values through a current output station to the local Foxboro 
temperature controller. The controller mainpulates the 
cooling water flowrate through the coils inside the reactor 
and brings the reactor temperature to the desired set point 
value. 

The experimental procedure described in section 2.5 has 
been followed. Conversion was measured using gas 
chromatography. Number average and weight average molecular 
weights were measured using the SEC. 

In total four policies were implemented, two from 
minimum time temperature policies and the other two from the 
minimum polydispersity temperature policies. The 
experimental results along with the theoretical results are 
presented in Figure 6.10 through 6.13. It can be observed 
that in all the cases the local temperature of the reactor 
closely followed the set point temperature provided by the 


computer. 
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The experimental conversion results agree closely with 
the model conversions in three of the four experiments in 
which a maximum difference of only three percent is 
observed. However in Figure 6.13 it can be observed that the 
experimental and model conversions agree well in the first 
one-hour of the reaction and afterwards the experimental 
conversion values are less than model values, the difference 
being about 10 percent at the end of the reaction. This 
discrepancy observed for conversion results in the Figure 
6.13 can be explained as follows. The parameters used in the 
model have been estimated based on the experiments conducted 
in the temperature range 65-75°C whereas in this particular 
experiment the temperature is varied from 63°C to 52°C. 
Since the parameter values are extrapolated, an error in the 
values of model parameters and consequently the model 
prediction of conversion values are expected. It should be 
noted that in the same experiment, in the initial period 
when the temperature was around 63-60°C, the experimental 
and model conversions agreed very well. For the other three 
experiments, the temperatures employed are within the range 
65-75°C used to estimate the parameters and hence very good 
agreement between model and experimental conversions are 
observed. 

In all the four cases, the experimental number average 
and weight average molecular weights are lower than the 
theoretically predicted values from the model. The 


experimental number average molecular weight values are less 
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by about 20 percent whereas the weight average molecular 
weight values are smaller by about 20-30 percent. However, 
though Systematic error is observed, the experimental 
molecular weight averages follow the trend predicted by the 
model in all cases. 

For estimation of the parameters of the model the 
molecular weight averages were measured using only two 
columns at a flowrate of 2 mL/min. At the time of open loop 
control studies some columns were replaced with new columns 
and an extra column was added to obtain better resolution. 
Also the flowrate was increased to 3mL/min so as to achieve 
the SEC analysis in about 10 minutes. All these changes 
could have contributed to a different resolving power of the 
column set compared to the original set used and hence a 
systematic error could have crept into the SEC analysis. 

Also errors in the calculation of M, and Mw may have 
resulted due to the following reasons. We assumed a 
semi-logarithmic relationship of the SEC calibration curve 
for the molecular weights less than 45000. The dispersion 
correction was based on the variance of the SEC 
chromatograms of the Marrow standards. The variance 
contribution due to the polydispersity of the narrow 
standards was not taken into account. 

Recently some investigators have reported experimental 
verification of the model results in polymerization reactor 
systems. Chen and Huang(1981) implemented optimal initiator 


policies on a batch solution polymerization of styrene. 
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Their experimental conversion results agreed very well with 
the theoretical results whereas their experimental M, and Mw 
values were in error with the theoretical predictions by 
about 20 percent. Cutter and Drexler(1982) carried out 
experimental studies to verify their simulation results in 
large scale polymerization reactors. Their experimental 
conversion values differed from model conversions by as much 
as 15 percent. The experimental M, and Mw values were in 
error with model results upto 25 percent. 
Pittman-Berger(1982) reported errors of as much as_~ 150 
percent between experimental and predicted values of number 
average and weight average molecular weights in _ batch 
copolymerization of styrene and MMA. Hamer(1983) carried out 
model verification studies in solution polymerization of 
vinylacetate in a pilotmeiscale —tubukamiereactor. His 
experimental number and weight average molecular weights 
were in error with model results to the extent of 35 
percent. Both of the authors attributed the discrepancy to 
the poor values of kinetic constants used in their models. 
From the above studies it can be concluded that the 
reported errors between the experimental and predicted 
number and weight average molecular weights fall in the 
range of 25-50 percent. The errors observed in this study 
which is 25-30 percent is at the bottom end of the range. 
The error reported in this study could have been improved if 
the same conditions for the SEC analysis could have been 


maintained. 


itiw ilew grey bese ativeot 
Wi ins. i Ladmemitaque riage ie: 
vi anoisalketa fapttencads add Aa r0r nb | 
‘wo 6etwao ASSET Pael xed Bria abst imma 
esluwmys mn itefunte nied qJitBy of eoihuse Ber: 
. Leanomt ogre vied qetiotonan eas somy So 
dsute ae. vd co] a raweere Letom soa boweitib woul 
et even. aso lees ol Gow iM Lo Setmsé tegen sit 3 - 
| «dears 25 o9qM | adloae: tebow 138 
an. AOL Be io. nibaee pan oqert (Saét }xsecon-nemaat i ° 
vedmon To eeyhin sasoibieng fine bagnomireas. resW eds; wa 
it. RE aoe lee amhigaton . epee. Wigiew ona 
geo Qe hosed (EGR i) ageeihy, ail pn onan yi a0 nok 
lo nolissisaayhog Aeidudos ad -apibusa roignad 
iH ..2OeDse2 tediedy dy ateom “gota: Ne iil ae 
pict iaw ve luositom: “pemtevs argon: paceoen nee 
ee to. whese, add, vel ea tuens Bebow. date | ‘WITS take e 
of oneueapebhl wid Aesudii2t9, sredipe sus 20 cdgpliw nena 
aLobeen he nd beat eanetanow alsenti. 6 asutev mee 3” / 
ois seig iebcddmmie) oil may 22, <aiboue svods ons moth: : 
bess theag Bain tesnaal eqns ods inital 2720378 ae 
odd Hk £tet aitigiew  walvaslom sapavers — ee 
you te aid? at havieado BxO3ts od neqeq oe -ag, to spret 
opcas she to Bae mostod oft ae @b Anoo19q. 0E-8R) at dotdw 
12 bevo jeune pisd eved. bluoo ybuite pins! ni. boxroges s9xts edT 
aiedt eved . bheoo: each 942 eft 103 enoltibsos emne edd 


bend stntam 


oa | 


i 
rig 
Li 


- 


E,, — 


176 


6.9 Conclusions 

In this chapter optimal control theory has been applied 
to the mathematical model of the batch polymerization 
reactor to derive open loop control policies. Four different 
control problems were formulated and solved. 

Some of the optimal control policies corresponding to 
the minimum time temperature policies and minimum 
polydispersity temperature policies were implemented on the 
experimental polymerization reactor system. The experimental 
conversion results agreed very well with theoretical 
predictions whereas errors in the range of 20-30 percent 
were observed between experimental and theoretical molecular 
weight averages. These errors are well within the range of 
errors reported by other investigators who carried out model 


verification studies in polymerization processes. 
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7. CLOSED LOOP CONTROL STUDIES 


Zak Imtreduction 

Ene chaptenjesixd openanloopaicontro) policies were 
developed by the application of optimal control theory to 
the model of the batch reactor. These policies were 
implemented on the experimental reactor system with no 
provision for feedback control as corrective action in the 
event of process upsets or disturbances or process-model 
mismatch. A more Suitable mode of control is to carry out 
closed loop Optimal | eEcntcoheal Fé i haetol emonitor @e polymer 
properties on-line and then take necessary corrective action 
should they deviate from the optimal trajectories. In this 
chapter we will develop and evaluate by simulation the 
feedback control policies by the application of linear 


regulator theory to the model of the batch reactor. 


7.2 Previous Work 

MacGregor et al.(1983) have reviewed the application of 
on-line control of different polymer properties in 
polymerization reactor control. Zeeman and Adams(1967) have 
presented an on-line control scheme which included on-line 
measurement of MWD and a digital computer. Roquemore and 
Eddy(1961) developed a discrete deterministic model and 
implemented feedback control schemes to regulate final 
conversion in a series of 9 to 12 continuous SBR emulsion 


polymerization reactors by manipulating temperature and the 
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number of reactors in the reactor chain. Fellows(1969) 
studied theoretically the regulator control of conversion 
and molecular weight in continuous SBR reactors where the 
feedback was subjected to random disturbances. 

Jo and Bankoff£(1976) studied experimentally the on-line 
estimation of states in the solution polymerization of 
vinylacetate in a CSTR using various forms of Kalman filter 
toebeonused “fore the ‘icontrolsmeot the: process? Hyun “and 
Bankoff£(1976) used a semiempirical model of the kinetics of 
a vinyl polymerization reactor to obtain estimates of the 
performance of a linearized Kalman filter in early detection 
of on-line process drifts. 

Ahlberg and Cheyne(1977) implemented a computer control 
scheme on a full scale continuous solution polymerization 
for the production of rubber. Some process model parameters 
were updated by using extended Kalman filter and the control 
action was based on the solution of the Riccati equation. 

Kiparissides(1978) solved a regulator problem of a 
continuous latex reactor. He derived suboptimal stochastic 
control policies by solving a linear quadratic optimal 
control problem and concluded that the control algorithm can 
be used to achieve better performance of the reactor. 
Kiparissides and Ponnuswamy (1981) applied the linear 
regulator theory to control conversion ina train of three 
CSTR-s. By using a locally linear state space model of the 
reactor, they obtained optimal temperature profile that 


brought the monomer conversions to the desired values over 
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the reactor train. 

Cluett et al.(1982) used a globally stable Adaptive 
Predictive Control System(APCS) for setpoint tracking and 
constant reaction rate control of a simulated PVC batch 
reactor. They obtained good, robust control despite 
nonlinear dynamics of the system. Kiparissides and 
Shah(1983) evaluated two types of adaptive controllers and a 
PID controller to a simulated batch suspension PVC reactor. 
They found that despite the time varying and nonlinear 
characteristics of the reactor, excellent control was 
obtained uSing either one of the two adaptive controllers. 
Both adaptive techniques proved to be very robust and 
outperformed a well-tuned PID controller under noisy 
conditions and time-varying parameters. 

Timm et al.(1982) studied feedback control of the 
molecular weight and production rate of polymer using the 
model of the styrene polymerization reactor and reported 
that PID controller worked well in tracking the setpoint 
changes in production rate and molecular weight. 
Schuler(1980) used a Luenberger observer and a Kalman filter 
to estimate non-measurable states in a simulated 
polymerization reactor. He also described how SEC 
measurement can be used to estimate the states that can be 
used for feedback control. 

However there have been only a few experimental 
evaluation studies with respect to on-line measurements and 


control of polymerization reactors. Meira et al.(1979) 
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Studied the control of MWD of the polymer product in free 
radical solution polymerization of MMA in a CSTR. They 
varied the feed flowrate of monomer and initiator ina 
feedforward manner and measured the effects of these 
variables on MWD of the polymer product and reported 
considerable success in controlling the MWD of the polymer 
product. 

In the following sections we will describe the 
development of optimal feedback control policies using the 
mathematical model of the free radical solution 


polymerization of MMA in a batch reactor. 


7.3 Development of Feedback Policies 

From previous experimental studies (chapter four) we 
know that temperature and initiator concentration affect 
conversion and molecular weight averages in the free radical 
solution polymerization in the batch reactor. The nonlinear 
model developed in chapter five can not be directly used to 
develop feedback control policies as the discrete linear 
regulator theory can be applied only to discrete linear 
systems. Hence it is necessary to linearize and then 


discretize the nonlinear model of the reactor. 


7.3.1 Linearization and Discretization 
There are two methods available for the linearization 
of nonlinear models. The first one, called apparent 


linearization, was proposed by Pearson(1962) and later 
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successfully used by Weber and Lapidus(1971) and McGreavy 
and Vago(1975). In this scheme the state equations are 
recast in a form that appears to be linear. The local values 
of the state variables are used to get the coefficients of 
the linearized equations. Kiparissides(1978) used this 
method to linearize the non-linear model derived for’ the 
case of emulsion polymerization of vinylacetate in a CSTR. 
This method is useful in situations where the steady states 
around which linearization is done is not known in advanced 
and is useful for batch processes where the operations 
values always change with respect to time. 

The second method uses the truncated Taylor. series 
approximations of the nonlinear terms. The linearization is 
done around known steady state variables of the system. The 
objective in this feedback control problem is to follow the 
nominal trajectories of state variables as closely as 
possible. The nominal trajectories are obtained by solution 
of open loop control problem as described in chapter six. 
Since the nominal values of state and control variables are 
known in advance, the linearization using Taylor. series 
could be used for this problem. Sage(1968) dealt with the 
procedure of linearizing around the nominal trajectories. In 
the following the Taylor series approximation is used to 
linearize the model of the batch reactor. The simplified 


model equations as shown in Appendix C are: 
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al 
——_ = -kyl Cite) 
dt 
aM 
—— .= -k,MyI 752) 
dt 
Ao 
= CkygI + k2MYI +k3SyI1 (73) 
dt 
Guz = k,SM? k 5M? 
eee et CK EM? (754) 
dt v1 VI 


Let=1,.,. Men #o. sana ws... Teter Jto the «nominal values of 
trajectories of initiator, monomer, zeroth moment and second 
moment of dead polymer concentrations respectively. Let T, 
be the nominal temperature profile. Linearizing the 
Equations (7.1)-(7.4) around the nominal trajectories in 
terms of deviation variables as denoted by A gives the 
following state-Space model in a vector-matrix form as shown 
in Equation(7.5) given in the next page. Equation(7.5) can 


be compactly written in the standard state-space form as: 


= AAY + BAu + WAI (726) 


AY = (4M ’ Auo, Au2)* ieee 
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Av i= eAT (7.38) 

Matrices A and the vectors B and W can be identified 
easily from Equation(7.5). The elements of A, B and W are 
dependent on the nominal values of state of the system and 
hence are time varying. 

For the purpose of digital control of the process, it 
is necessary to discretize the linearized model of the 
reactor. Discretization of the continuous differential 
equations has been treated by many authors (Ogata(1967), 
Franklin and Powell(1980)). 

After discretization the continuous linearized 


Equation(7.6) can be written in the standard form as: 


AY ed, 2 BPAY pot GrAws 4+ THAD, (7.9) 


where the discrete state-space coefficient matrices F, G, 


and H are given by: 


F = exp(At) (72-80) 
G = J exp(Ar)Bdr (7.11) 
H = j exp(Ar) Wdr C7ate) 


and t = sampling time 
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To check the validity of the linearization and 
discretization procedure, the responses of the linearized 
and discretized models were compared to the response of the 
Original nonlinear model when disturbances in the initiator 
concentration from the nominal values are introduced. The 
nominal values obtained from _ the open loop optimal 
temperature policy developed in section 6.4 has been used. 
The discretization is done based on a sampling time of five 
minutes which is the expected SEC analysis time when 
efficient and fast columns are used. 

First an increase in the initial initiator 
concentration by 0.025 mol/L from the nominal value (0.05 
mol/L) is introduced. The conversion and number average 
molecular weight responses to the initiator changes 
predicted by the nonlinear, linear and discrete models along 
with the nominal responses are presented in Figure 7.1. It 
can be observed that conversion and number average molecular 
weight responses almost coincide for the linearized and 
discretized models whereaS a small discrepancy exists 
between the nonlinear model response and the other two model 
responses. These results confirm that the linearized and 
discretized model represent the nonlinear model quite well. 
In Figure 7.2 the conversion and number average molecular 
weight responses are given for the case where the initial 
initiator concentration is decreased by 0.025 mol/L from its 
nominal value (0.05 mol/L). Again the responses confirm that 


the linearization and discretization procedure represents 
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Initiator Concentration Change=-—0.025 Mol/1 
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the nonlinear model quite well. 


7.3.2 Control Law Design 

Our objective is to design a feedback controller for 
the batch polymerization reactor so that the states of the 
reactor track the nominal trajectories of states despite 
process upsets, disturbances, process model mismatch etc. 
The nominal trajectories were previously obtained by solving 
the open loop minimum time or minimum polydispersity problem 
as described in chapter six. 

One can formulate a quadratic objective function for 


the closed loop control as: 


mised) =fAYSZAY, b> LiebcAYeOQAY, tHeR Aus) 
k =0 CTs) 
The following physical interpretation can be given to 
the objective function (Kirk(1970)). It is desired to 
maintain the state vector close to the nominal or desired 
values without excessive expenditure of control effort. The 
weighting function matrices Zz and Q are positive 
semidefinite and R is positive definite. If R is singular, 
the control effort would become infinite. 
The above is a well studied problem in discrete 
linear-regulator theory. The control law that results from 
the minimization of J with respect to u, is: (Sage(1968), 


Franklin and Powell(1980)): 
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Au, = K fi (7614) 
where 

K = UR PuGuce) ! €P, -O) CIS) 

PP. = O-rer 2: .., (1 + GR°"G'*Pk+1)F G7 16) 

P,.¢ <=«Z C7T tT) 


The elements of K which are time varying gains of the 
feedback controller can be computed from Equation (7.15). 
The elements of the 3x3 matrix P, in Equation (7.16) can be 
obtained by solving the discrete-time matrix-Ricati Equation 
(7.16) using the boundary condition given in Equation 
(7.17). It can be observed that the gain matrix K can be 
calculated off-line by using only the nominal trajectories 
of the states. These time varying gains can be stored ina 
computer and then used to implement the feedback control 
after measuring the values of the state variables as shown 


int Figure’ 7.3% 


7.3.3 Simulation Results 

Closed loop control policies were derived using the 
discretized model of the reactor. For all the cases the 
nominal control and state variables corresponding to _ the 
minimum time policy developed in chapter six was chosen. 
Although initiator concentration could have been chosen as 


another control variable, in this study it is treated as a 


eer ea i As Siren 
" to Hanae Re tee ys ; 
Pa We fe 


Cer .d) ) a 


ri) ‘ ; x 
2 (o- ia 
Ss, 1 pa _ my : 
re al) Ane ae oem % Se eh , 
iar.y)C- shrergro! ta + b) ‘ rt 4 +" 9 ot gr 
“a 2 7 mle re. — | ie 
ae - ; a a 7 ‘ 4 Ly = Ge 3 ' 


ade |p antes: gnieaey) Sint ste ‘Hidde ® a0 B: ; 

v} WO ES GLE nowt Bastygno> od” sang sitet 

ad nee bate) wo aap eh (hae ey eds Yo’ 

moet wane byes tax siey, “anks-eseroett alg gatvion a 

iotrenps ni weavig pemINR BOD: yeabaived ; od rine (a 

neo t eave wiape eit rads savaisede ae at ithe 

aehrovsntgs? Lain pla oat eho: ental geen ote 5 

it: hexose ad nae ‘prime oni¥sey re oaok? vasgage s #3 

leoxgaes Avaaltve it ata “Ponwmised et votibeen ned! bos 168 | 

yea seldalaan dae ‘eit do ae ssc 

| peeieen was - 

o77° pba. bovbiss ‘e304 sotodten terse goat: ‘eenreinaunee! a 

oA estes. .olt? Lie 2 Stodose1 » ada ten Labom! bosidersetb | 
ed: of  Biidbacgeestes  geldsixev pegete bne forgnes! Leaimon 
 MOROS cee niz Tataalts at \beqelewet: qoifeq omig mumi nim 
#¢ nsaots ‘ceed eval bigod moldersnasnos roseistiat dpvoislA 


a 
oh 


‘wy os Beteeds oi d¢-yOure eidz as vafdsiicy fovinos 1si20na 


190 


disturbance. Temperature is treated as the control variable. 

In Figure 7.4 a disturbance in the form of increase in 
the initiator concentration has been introduced. It is known 
that in free radical polymerizations, increases in initiator 
concentration reduces the molecular weight of the polymer 
and increase the rate of polymerization. These effects can 
be observed in the responses of the uncontrolled reactor in 
figure 7.4. When the reactor is put under feedback control, 
the control variable temperature decreases from its nominal 
values to counteract the effect of the disturbances. It can 
be clearly observed in Figure 7.4, that the conversion and 
number average molecular weight responses move towards’ the 
nominal responses. 

Figure 7.5 presents the results when a disturbance in 
the form of decrease in initiator concentration is 
introduced. Again it can be observed that the feedback 
controller brings the responses to the nominal desired 
values. In Figure 7.6 the feedback control results are given 
for the case where initiator disturbance is introduced 20 
minutes after the start of reaction. Results similar to the 
previous cases can be observed. 

In all of the feedback control results it can be 
observed that there are offsets between the nominal values 
and the controlled responses. From Equation (7.14) and 
Figure 7.3 it can be seen that the optimal feedback 
controller has the structure of proportional controller. 


Since integral action is absent, the offset is not entirely 
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Fig. 7.4 Closed Loop Control of Polymerization Reactor-1 
Increase in Initiator Level=0.025 Mol/1 
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eliminated by the proportional controller. However, integral 
action by introducing additional state variables can be 


easily incorporated to remove this offset. 


7.3.4 The choice of Q and R 

Q is a 3x3 matrix and was chosen as a diagonal matrix. 
The q33 term was set equal to zero to indicate that we did 
not wish to explicitly control deviations of Mw (weight 
average molecular weight) from nominal values. The remaining 
diagonal terms q;; and qz22 terms were chosen after numerical 
experimentation. They also show the relative importance 
given to the deviations of monomer concentration and zeroth 
moment from nominal values. Final values chosen were q,;=1.0 
G22=10°, and q33=0. In the proposed run R is a scalar. The 
value of R decides the magnitude of control action. A _ zero 
value of R_ results in large control action and drives the 
controller to the limits of saturation. In our case we used 


a value of R=107°5. 


Stochastic Estimation and Control: 

Even though we developed control strategies only for 
the case of deterministic model of the reactor, a realistic 
model would be stochastic in nature. The model prediction of 
states and the measurement of states could be used to design 
a Kalman filter. The state Estimation and regulator theory 
could be used realistically to design the feedback control 


of the reactors. 
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7.4 Conclusions 

In this chapter the nonlinear model of the reactor was 
linearized around the nominal trajectories obtained from 
open loop control policies. The linearized model was further 
discretized andthe optimal regulator theory was used to 
derive the feedback control policies. The simulation results 
show that the feedback control counteracts the effects of 
disturbances and brings the conversion and number average 


molecular weight back to the nominal values. 
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8. CONCLUSIONS 

The on-line measurement and computer control problems 
associated with the solution polymerization of MMA in a 
batch reactor have been studied in detail in this thesis. An 
experimental batch polymerization reactor with facilities 
for on-line measurements of density, viscosity, torque on 
the stirrer and molecular weight of the polymer was built. 
All these meaSurements were interfaced to a HP/1000 process 
control computer. An automatic sampling system was designed 
to collect samples of the polymerization mixture 
periodically from the reactor and to automatically inject 
them into the size exclusion chromatograph (SEC) after 
dilution with a solvent. Computer programs were developed 
for automatic acquisition of SEC data and processing of 
these data to obtain molecular weight averages of the 
polymer. 

Conversion measurements were made using three different 
techniques: gravimetry, gas chromatography and densimetry. 
The analysis of the results showed that gravimetry showed 
appreciable error in conversion measurement particularly for 
high conversion samples that are viscous. The densimetry 
results agreed very well with the GC results that were shown 
to be more accurate than the gravimetric results. It is 
recommended that on-line densimetry be used to monitor 
monomer conversion during the course of polymerization. The 
analysis of the SEC results showed that the number and 


weight average molecular weights could be measured with 
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reasonable accuracy during the course of polymerization. 

A mathematical model of the solution polymerization of 
MMA in a batch reactor has been derived. Since the 
parameters available in the literature did not predict the 
experimental results well, these model parameters were 
re-estimated using our own experimental data. Verification 
of the model (with the re-estimated parameters) was 
demonstrated to be excellent by comparing it with the 
experimental data. 

Optimal control theory was applied to the reactor model 
to derive optimal temperature and initiator concentration 
policies for four different problems. The optimum policies 
corresponding to minimum time and minimum polydispersity 
temperature profiles were implemented on the experimental 
reactor. The comparison of model and experimental results 
showed that conversion results followed the model results 
very well whereas discrepancies of the order of 25-30 
percent were observed in the number and weight average 
molecular results. These errors are well within the range of 
model prediction errors reported by other investigators in 
this research area. 

Linear regulator theory was also applied to the 
linearized model of the reactor for the design of closed 
loop optimal policies. Implementation of these policies on 
the simulated reactor showed that the feedback controller 
was able to track the nominal trajectories of conversion and 


number average molecular weight in the presence of 
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disturbances due to initiator concentration changes. 

As an overall conclusion this study demonstrated that 

with a proper choice of on-line measurement instruments it 
is possible to do advanced control of a _ polymerization 
reactor. 
Future Work: This project was started to provide the 
necessary groundwork for a series of application studies in 
the area of on-line measurements and computer control of 
polymerization reactors. Though closed loop control was 
studied only on a simulation basis, computer implementation 
and verification of these policies are recommended for 
future work. An energy balance of the reactor can be 
performed from the meaSurements of hot water and cold water 
flow rates and their temperatures. This energy balance 
information could be used to predict and control the 
reaction rate during the course of polymerization. The 
experimental system can be easily modified to. study 
continuous polymerization reactions. Also the reactor system 
can be used to study on-line measurements and control 
problems associated with emulsion and suspension 
polymerization reactions. 

In this study only optimal control theory has’ been 
applied to optimize the operation of the polymerization 
reactor. Since the polymerization reactor is highly 
nonlinear and contains time varying parameters, adaptive 
control techniques such as self-tuning regulators (STR) and 


self-tuning controllers (STC) could be applied to study the 
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APPENDIX A: EQUIPMENT AND INSTRUMENTS REQUIRED 


Reactor System 
1. Five litre cylindrical jacketted glass reactor. 


2. Hot water tank.(25 litres capacity, 1500 Watts) 
Blue M Electric Company, Blue Island, Illinois, USA 


3. Hot water circulating pump (centrifugal) 
MMiVolts S20 Watts, 3000 rpm. 
Littie=Giantpump Co; Oklohamar city, OK,US. 


4. Rotameter (Hot water flowrate) 
0-2 gallons per minute 
Brooks Instrument Division 


5. EMF Converter (millivolts to 10-50ma signal) 
MIS ivolts 960 cycles, 9 watts. 
Foxboro Company, Mass.,US. 


6. Foxboro M/62 Controller.(Model 621-5L) 
Input(10-50ma) -Output(10-50ma) 
Foxboro Company, Mass., US. 


7. Electro Pneumatic Transducer(10-50ma to 3-15psi) 
Foxboro Company, Mass.,US. 


8. Control Valve (for cooling water) 
Research Meter Inc, Tulsa, OK, US. 


9. D/P Cell(To measure cooling water flowrate) 
Foxboro Company, Mass., US. 


10 Densitometer(DMA 45) 
Anton Paar Ltd, Austria 
Distributed by Mettler Ltd, USA. 


11. Stirrer-Motor Assembly. 
130 volts. OvaAmpo/San.p- 
G.K. Heller Corporation, Floral Park, New york. 
12. Liquid Chromatograph (Size Exclusion Chromatograph) 
Model 244. 
Waters Associates, Mass., US. 


13. Gas Chromatograph(Model 5710A) 
Hewlet-Packard C 


14. Viscometer(Model M8B) 
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Norcross corporation, Newton, Mass., US. 


15. Reciprocating-Revolving Pump. 
Fiuvd Meteringeincs, N.Y .)@US? 


Automatic Sampling 
it 50 ml conical flask 


2. Solenoid Valves. 
Skinner Valves, CT, US. 


3. FMI Pumps(See Item 15 Reactor System) 
4, Magnetic stirrer 


5. Automatic Sample Injection Valve 
Valco Instruments Co., Houston, Texas, US. 


6. Programmable Sequence Timer. 
Potter and Brumfield, Guelph, Ontario, Canada. 


Computer System 


1. HP/1000 Minicomputer 

2. LSI 11/03 as 1/0 subsystem 
3. HP Plotter 

4. Terminals 
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APPENDIX B: SEC DISPERSION CORRECTION 


Let Mw(v) and M,(v) be the weight and the number 
average molecular weight respectively of the polymer within 
an infinitesimal fraction at elution volume,v. 

If there were no column dispersion i.e at infinite 
resolution the weight concentration and the molecular weight 
at any SEC elution volume, y, are Wy) and M,(y) 
respectively, where M,(y) is the true calibration 
relationship. If we let a Gaussian function G(v-y) describe 
the fraction of species at y which gets spread over to v, we 
can consider the molecular weight mixture at any v to. be 
made up by species coming from all different y elution 
volumes. The molecular weight of these species are M,(y) and 
their weight concentrations detected at v are W(y)G(v-y). 

By definition of the weight average molecular weight, 


we can write the Mw of the mixture at v as follows: 


+o 


JWly)G(v-y)M, (y) dy 


Mw(v) == (Bit) 
+o 


J W(ly)G(v-y)dy 


where 
G(v-y) = —— exp[ - ——— ] CB 2) 


and 
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M.(y) = D,exp(-D2y) (Bia) 


In Equation(B.2) o is the peak standard deviation of the 
assumed Gaussian distribution. In Equation (B.3) D, and Dz, 
are constants of the true molecular weight calibration. The 
denominator of Equation (B.1) gives the overall weight 
concentration detected at elution volume, v, which is F(v). 


Hence F(v) can be obtained by the following equation: 


+o 


F(v) = JWly)G(v-y)dy (B.4) 
—-& 
Let us now evaluate the numerator of Equation (B.1). 
Substituting for G(v-y) from Equation (B.2) and for M;,(y) 
from Equation (B.3), the numerator of Equation (B.1) can be 


written as: 


ee (eo4 (vey) 2 ) 
N = f W(y)(——— exp[- —————])Diexp(-Dzy) dy 
ae ( of2n on? ) 


After considerable manipulation the above equation can be 


evaluated as: 


N = M,(v)exp{(D20)7/2} F(v-D207) (B.5) 
substituting Equations (B.4) and (B.5) into Equation (B.1) 


we obtain: 
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F(v-D.07) (D507 
Mw(v) = —————————. exp { ————__},, (v) (B.6) 
F(v) 2 


By using the definition of the number average molecular 


weight, the number average molecular weight of the polymer 


sample eluting at the elution volume, v, can be given as: 


+o 


J Wly)G(v-y)dy 


—c 


My, (v) = ——— (Bay 
+o 


f{wly)/M, (y) }G(v-y) dy 


The denominator of Equation (B.7) can be evaluted by 
substituting for M,(y) and G(v-y) from Equations (B.3) and 
(B.2) respectively. Let denominator of Equation (B.7) be D. 


Then 


(Vey jc 
W(y)exp(- ————) 
+o oV2n 
D = J SY 
—o of2na D,exp(-D2y) 


The above equation can be evaluated as: 


1 (-p-0) + 
exp[———__] F(v+D207) (B.8) 
ee 


De 
M, (v) 


Substituting Equations (B.8) and (B.4) in Equation (B.7) the 
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instantaneous number average molecular weight at the elution 


volume, v, iS given as: 


F(v) (Dea) * 
M,(v) = —M————— exp{ ———— } M,(v) (B93) 
F(v+D20) 2 
Since F(v) is the weight fraction eluting at elution volume 
v, and M,(v) and Mw(v) are the number average and weight 
average molecular weights respectively at elution volume, v, 


the overall number and weight average molecular weights can 


be obtained by the following relations: 


+o 
JS F(v)Mw(v)dv 


—ocO 


Muwite ee (By 10) 
+o 
J F(v)dv 
+o 
f F(v)dv 
Mon at aoe Oras 0, 
+o 


S{F(v) /M, (v) tdv 


Equations (B.10) and (B.11) give the number and weight 
average molecular weights after the dispersion correction. 
These corrections are implemented in the computer program 
developed for the analysis of SEC data. 

For constant o7, Equations (B.10) and (B.11) can be 


integrated analytically and simplified further. 
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Let M,(uc) and Mw(uc) be the number average and weight 
average molecular weights of the uncorrected chromatograms, 
respectively. Then, Equations (B.10) and (B.11) can be 


rewritten as 


Mw = Mw(uc)exp{-(D20)?/2} (Bitz) 


K< 
= 
i] 


M, (ucyexp{ (D,0) 2/2)} (Bi 13) 


where 


oo 


SF(v)M, (v) av 


—c 


ot) 


i 
wo 


Mw (uc) 


jos) 


JF Cv dv 


—c 


oO 


J F(v)dv 


=o 


M, (uc) (B.15) 


SiF(v)/M, (v) fav 


Equations (B.12) to (B.15) could be used to correct for 


the dispersion effect. 
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APPENDIX C: SIMPLIFICATION OF REACTOR MODEL 


The mathematical model of the solution polymerization 
of MMA ina batch reactor has been derived in chapter five 


and is given as follows: 


Gifs = ke CO) 

aM/dt = -k,M/I (Cs 2) 

dug /dt = “Ck SRAM) Ao t(Rygtiar/ 2) io 7 (Ge3) 

ite /dtie= SOS Tt kyM) Xo Paki ANb As + kee AQ? (C.4) 
where 

Neate V2Ekgiy ks (oa) 


ZEK gi + CKOMtKAM Fakes) AG 
1, =O C6) 
kp Met—k 2S hiekoteg:) 
rA2 = Aat 2kpMA ,/(KkmM + kes ws Keto) (C7) 
Thomas(1981) carried out an order of magnitude study 
for the terms appearing in Equations (C.6) and (C.7). In 
Equation (C.6) it was shown that the term (2fk,I) is 
negligible compared to (kpMAo). Neglecting the rate 


constants for transfer reactions, the following expression 


for A; is obtained: 


rA4 =kpM/k , LCS) 


Similarly an expression for dz can be obtained as: 
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kido 
Substitution of Equations (C.5), (C.8) and (C.9) into 
Equations (C.2) through (C.4), the model is simplified to 


the following set of equations. 


aI 
eer eh! CE 1:0)) 
dt 
dM 
Nes okyMy.l (6011) 
at 
duo 
= Chel a k 2MyI +k3Sy/I (C.12) 

dt 
dus = k,SM? k 5M? 

—————- + —— + kM? Cees) 
dt VI VI 

where 

k; = kpY2fka/k; = A,exp(-E,/RT) (Gata) 
kgit= kpwV2ZEka/k, = A-exp(-B2/RT) CCRAt 5) 
k 3 = Kee V2Eka/k, = A3;exp(-E3/RT) COAG) 

2kp?k, 
ky = ———————— = A,, exp (-E,, /RT) (Cai) 


B2ZLKek, hse 


ae. Oe, “ie " a * 
a he : (a 


tt 


ai ) | See CR yy anat pe 30 


Su oe a 


hs my - high ; i | rhs (e.9) | F 


‘i iimie &e Demi tet? 


te ~ 


+ 


ro on 
“Maal * c<oe eebeeo * a . 


( | 


4} | \ (TR eee -aNieas vad = | 
pal.) (TH\ git lgepiga = ) RTL RS ae © gt 
BW. . OT aD\ et axe = (Ry hw ; | 


i ae 


(TH y Se Toeeg A | ened ant * | 
; Wie 


223 


2kp? Km 
ks = ——————_- = As exp(-E;/RT) LE. 8) 
(Otek, s 4 
(2+v) kp? 
k ppeodeet— tease tragexp (E67 RP) (C.19) 
Ke 
C= 24 t=v7 2) p= kao Vkh Ceq 208) 


The above simplified equations were derived by the order of 
Magnitude study on the terms occuring in the expressions for 
Ay anders. 

Further simplifications of the model equations can _ be 
achieved if the transfer reactions are neglected in the 
model in which case the terms kz, k3, ky and ks become zero. 


Then the model equations are: 


aI 
——— = —kd! COZ) 
dt 
aM 
Ry LM (G22) 
at 
dio 
— Ckyl (C. 23) 
at 
du2 
= k,M? (C.24) 
dat 


These models are used to derive the open loop and control 


policies in chapter six and also the Appendices D, E and F. 


as 


JR oF bs - 
. oe | 
sic vd Soot) wpe ee ee i 


joO¢esasges ory pig roe earigainde ag wboae | 


oo | 


bas yelpar even, anqiaces om 
ener a 5ns at od ee one me ones » sotite me , 


: . Send : 
: ; stings bor at rao ai etn yam me 


APPENDIX D: MINIMUM TIME TEMPERATURE POLICY 


The objectiveis to find the optimal temperature policy 
for a. given initial initiator concentration to produce a 
polymer product with desired conversion and number average 
molecular weight in minimum time. 

The simplified model equations have been derived in 


Appendix C. They are: 


aI 
— eee CD.st) 
dat 
aM 
—————s (= key ft (De2) 
dt 
duo 
= Ckgl i+ k3ZMyI o+k,Sy1 (De3) 
dt 
du2 k4SM2 k5M? 
= Sti a2 k eM? (D 4) 
dt VI v1 


The desired conversion X* and desired number average 
molecular weight M,“ can be expressed in terms of desired 


monomer concentration (M*) and desired zeroth moment of the 


polymer(uo*) as: 


M* <= Mo.( 1-35) (D.5) 


Mo™ = MoX* MW/M, * (D.6) 
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where MW is the molecular weight of the monomer and Mo is 
the initial monomer concentration. If tf is the final time 
and M(tf) and wuo(tf) are monomer and zeroth moment of 
polymer concentration respectively at the final time, the 


following objective function can be formulated. 


we 
min J= wiiM(t£)-M*}? + WolHo ttl) —koyt? t octet (P27) 
0 
where w, and wz are weights. 


The Hamiltonian for the above optimization problem can 


be written as follows: 


H=1-P,kyI-P2k,MYI+P3{CkyMVYI+k3SyI1 } (D8) 


In Equation (D.8) P,, Pz, P3 are co-state variables which 


must satisfy the following canonical equations. 


GP, 0H 1 P2k,M 1 keM + k3S 
sees EIS =kyP,+ -_ —P.(Ckaet- tic ) (D..9) 
dt oI 2 VI 2a VI 
aP2 OH 
wh we ky PeVIA eee (De 1.0) 
at OM 
aP; 0H 
a’ = 0 CDi. Th) 
at dio 


The co-state variables will also satisfy the following 


transversality conditions: 
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P,Gef) = 0 (D.12) 
oJ 

P,(tf£) = ——] = 2w,{M(tf)-M*} Gaya. 
OM |tf 
oJ 

P,(tf£) = =<2wis hits (tik-ug tt (D.14) 
Oho Cf 


The’ \«intrtial sand thinahcconditions sforcthe-state sand .co-state 


variables for the above problem can be written as: 


variable t=0 t=tf 
I le free 
M Mo free 
Lo free 
P, free 0 
P2 free 2w, (M(t£)-M*) 
P3 free 2W2 (ite (CE) no: 


If the temperature is constrained, the optimality 


conditions that must be satisfied by the optimal policy are: 


H = 0 
dH/dt < 0 Bae ae dM Mio) 
dH/oT > O es 
dH/oT = 0 T is in the interior (D.15) 


where TT" and T, are upper and lower limits of the 


temperature 
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A two point boundary value has to be solved with the 
Simultaneous satisfaction of the optimality conditions to 
obtain the optimal temperature profile. 

A discrete control method algorithm (Thomas(1981)), 
Thomas and Kiparissides(1984) is used to solve the above 
problem. It should be pointed out that the minimum time 
problem is not directly solved in this method. First a 
series of fixed terminal time problems are solved and _ then 
the smallest time for which convergence is possible, is 
selected as the minimum time. 

The essential steps of the discrete control method are 
given under: 
1.Select a fixed time. 

2. Guess a piecewise constant temperature profile T(t). 

3. Using T(t), integrate state equations forward. 

4. Compute P,(tf£), P2,(tf) and P,{tf) from Equations (D.12) 
to (D.14) and integrate the co-state equations backward; 
simultaneously compute g(t)=0H/0T 

5. If the final conditions of state variables are satisfied 
within a tolerance limit, stop. 

6. Otherwise compute new temperature history using the 


following equation 


Coy =T Ct) = eag Gti) 
new old 


where a is selected after some numerical experimentation. A 


small value of a would lead to long computing times’ before 
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reaching the solution and a large value can lead to 


oscillations or even divergence from the solution. 
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APPENDIX E: MINIMUM TIME TEMPERATURE POLICY SIMPLIFIED MODEL 


AS given in Appendix D, the objective in this problem 
is to find the optimal temperature profile at a given 
initial initiator concentration to produce a polymer with 
desired final conversion and number average molecular weight 
in minimum time. 

In this development the simplified mathematical 
obtained by neglecting the transfer reactions are used. (see 


Appendix C) The simplified model equations are: 


di 

= led (E 1) 
dat 
aM 
—— = -k,M/I (Eee2) 
dt 
dio 

SUC ky tk (B23) 
dat 


The objective function of this problem is: 


TE 
min tf = max -f -—dt (E.4) 
0 


The Hamiltonian for the above problem can be written as: 


H = -1-P,kgI-P2kiMVYI+P3Ckgl (E.5) 
InvEquation (B.5)P9, P2, and Ps are costate variables which 
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satisfy the following equations. 


GP, 0H Po kM 

See P,kgt — P3Ck, (E.6) 
at OM 2b oft 
dP2 0H 

=- —— = P2k,yI CEST) 
dt OM 
GaP; 0H 

=- = 0 (E.8) 
dat dio 


The optimality conditions assuming no constraint on the 


temperature are given by: 


H=0 (E.9) 
oH P,Eygkgl a, P.E,k,MyYI it: P3;CEaky 

5 ———— } = 0 (E.10) 
aT RT? 


From Equations(E.2) and (E.7) we obtain 


aM aP, 
+ = 0 (E.11) 
M P» 


From Equation(E.11) the following equation is easily 


obtained. 


P.M = constant = a CBrst2) 


From Equation(E.10), Equation(13) can be derived. 
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P;CBak, = P,Eakgl + P2E,k,MyI1 (Es) 


Substituting (6.13) in (E.9) we get: 


P,MVIk,(B,-B,y)-Ea = 0 (BE. 14) 


From Equation (E.12) and (E.14) we get: 


Eq 1 
k,VIl = ————— — (EB. 15) 
(E,-Ey) a 

Equation (E.15) can be used to get the optimal temperature. 

The optimal temperature can be obtained by the following 

procedure. 

4. Guess the value of Pz at time t=0 

5. Calculate a = P2M = constant at time t=0 

6. Calculate the temperature T from Equation(E.15) 

7. Integrate Equations (E.1)-(E.4) forward and at each step 
calculate temperature T. Proceed until the desired 
conversion is reached. 

8. When conversion has reached the desired value, compare 
the number average molecular weight with the desired 
value. 

9. If they are within the tolerance value, stop 

10. Otherwise go to step 1 and adjust the value of P2 at 
time t=0 

A bisection method is used to iterate upon the initial 


values of P, in this program. Computer programs have been 
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developed to solve the above problem. The results are 


presented in Section 6.6. 
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APPENDIX F: TEMPERATURE POLICY FOR MINIMUM POLYDISPERSITY 


This Appendix solves the problem 4 described in section 
6.3. The objective is to produce a polymer product with 
desired conversion (x*) and number average molecular weight 
and having minimum polydispersity. 

The model equations required after neglecting the 


transfer reactions (see Appendix C) are: 


dI 

= -kyl CE a) 
dt 
aM 

= -k,My/I1 CE 2) 
i 
duo 

= Ck (F383) 
dt 
du 2 

= kM? (F.4) 
dt 


Note Equations (F.1) and (F.3) are linearly dependent and 
hence uo can be calculated from I. 


The objective of the problem can be written as: 


Ming GbE) (F.5) 


The Hamiltonian for the above problem can be formed as: 
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H=-P,kyI - P2k,MVI + P3k,_M? (FG) 


In Equation (F.6) P;,, Pz and P3; are co-state variables which 


must satisfy the following equations: 


aP, oH P2k iM 
Soe Pk (Ee?) 
dat dl 2u i 
aP2 0H 
= - _—_ = P2k,yVI -2P3;k._—M (F.8) 
dt OM 
GaP 3 0H 
ee) (F.9) 
dt OMe 


The initial and boundary conditions of the problem can be 


written as follows: Boundary Conditions: 


Variable t=0 t=tf 
I Io I o-Mox*MW/CM, * 
M Mo My (1-x*) 
be 0 free 
P, free free 
P2 free free 
Ps =i =" 


(Faro 
Along the optimal path the following two equations will hold 


good. 
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= ( =P,Eakal -P2E,k,MVI+P3E,k,_M2 }=0 Ce 1) 
oT RRA 
H = -P,kgI--P2k,;MYI+P3k_M? = 0 (P 5 42:) 


Substituting Ps=-1 and Equation (F.12) in Equation (F.11) 


and rearranging we get: 


Pok,M/I{E,-E;} ae KeM24B eka = 0 CE is) 


But we know thats: 


kg = (2+v )kp?/k. = Asexp(-E,/RT) (F.14) 
k, = kpV2£kg/k; = A,exp(-E,/RT) (F.15) 


With a value of f=0.5, from Equations (F.14) and (F.15) we 
get: 
(B.-Ea) = -2(E,-E,) (F.16) 
Ke = (24p) ka 7/k. CER 7s) 


Substituing Equation (F.17) in Equation OMe P28) and 


rearranging we get: 
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Ky 2M 


Se i eee (F.18) 
ke P2VI 


From Equations (F.14), (F.15) and (F.18), the optimal 


temperature T can be obtained as: 


(E,-E, ) 
ee (F.19) 
-2A 5M 
R in( 
A,P2vVI 


The optimal temperature profile can be calculated by the 
following procedure. 
1. Assume P2(0) 
23 “Calculate T from Equation (F.19) 
3. Integrate Equations (F.1), (F.2) and (F.3) until desired 
conversion x* is reached. 
4, Check if M, has reached the desired value. 
5. If M, has not reached the desired value go to step 1 and 
adjust P2(0o). 
By using the simplified model, we avoid the solution of 
the two point boundary value problem. We have to iterate 


only on the initial value of the co-state variable P2. 
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APPENDIX G: EXPERIMENTAL RESULTS 


ac 6 Ty _ ; i : 


i : ri ee . 
a 


- i Na . 
ae. ey ye y : 


STUART JT eT 


7 


t t 
1 : ‘a 
on TT at 
ale 
i 
u 
f 


= 
: 
3 7 
i> 
4 j 
5 U 5 Jt We. , 
} au 
j 4 > 
“ 
‘ 
} 
a 
i 
j 
4 
1 
i 7 
' 
4 
wf 
i 
i 
ty 


238 


T/1OWSO°O =NOILVALNAONOD YOLVILINI IILINI 


9°O =OILVY LN3SA10S 


9,°S9 =deNLVaAdWw4al 


Gs) 


AYLAWIAVAD 
IN3A0¢4d 


ALISNSQG ALISOOSIA NOTSYSANGO 


b 


=ON LNAWIdSdx4 


Vos “OOE 
b 9S “GLE 
© 6S “OSC 
LE OSv. SCc 
v'SY melee 
6° OV SLI 
v'9E oct 
ELE Sch 
b 96 “OOt 
Lobe O08 
9° Lt 09 
6 bt “OP 
Gal “OC 
e}0) NIW 
LN30d4d NI 


NOTSdYSANGO SWIL 


T- eaessa Pease 


oS PETES Waele 


EVAN gE ROTATE LAT ees 


Zog 


TV/1OWOL*O =NOTLVALNAZONOD AOLVILINI IWILINT 


9°O =OILVY LNAA10S 


Oo, G9 =3YNLVaAdW3AL 


AYLAWIAVAS 
LN30d3d 


ALISNAQ ALISOOSIA NOTSASANOO 


c 


=ON LINAWIYddx4 


Ue 2S, “OLE 
Be) “Ove 
v°c9 “OCC 
t° 6S “00% 
6°0S Ost 
c°bkG “O9t 
c¢ Ot Ort 
by’ CcV Oct 
c' LE “OOt 
bE “O08 
G 9¢ O09 
So Lt OV 
O'bt O¢ 
rO)2) NIW 
LINd0d3d NI 


NOTSYSANO9 SAWIL 


se 


ami f 
NS See: 


- AL ee 


&a2 


eae 


Se ek ~ ew an ane fee ee 


os 


ae 
. © ¥? 


tet 
m. 
ats aa 
= 


<i TRU TEATS 


Ag PAS Sev ste 


val 


240 


TV/1OWGO*°O =NOILVALNAONOO YAOLVILINI QVILINI 


9°O =OILVYe LN3ZA10S 


Oo, GL =3aeNnlVaddwal 
8L56'0 Seo b S8 Ease Ost 
Orsce'O SwLES €° O08 8° SZ o9t 
G6r6°O SOONG G OL tock Ovt 
EEvE O Sec SZ E99 Oct 
CLEB O 8° 6t 6 £9 Smeos OO} 
E8c6°O vet iP (SiS € OS 08 
88l6°O Sg €°vv CCV 09 
GLO6°O E58 cle ple OV 
vS68°0O 9 O° Lt L’° 6k O¢ 
an 2 MLAWTAVID 5 NIW 

LIN30du3d IN39d4d NI 
ALISNAQG ALISOOSIA NOISYSANOO NOISYSANO9 AaWIL 


€ =ON LNAWITdSadx4 


pene teed Od aes See ae eee eRe eee ee ee eee=- 
i Sie : we 2. wa x 
- —— y, ee Soe 9 otek a gpa ep celine mepapted 


=a oe : ot eh ox Set - 


| 
Ss 


art 
; WI 
aa 


= 7 - tone 2 a s.hr ane te oe 
- : | : a ere % 8.3 : ; ean Ee eo 
/ 2 s ae eS c908.o pt. ee ead 
7 a oe : . ne ee Sn ci 
ma — a feue 6 ee 


oe . ee ¢. 93 tar 


; a a5 a bet ad 
Z f = a7 es =. —= Sg ee a 


: ) . / = eae 


=.= ot ov. BG =fiT%* ea) a ee 


sO) FT aos Nes Bet eI het 


241 


T/1OWOL’O =NOTLVYLNSONOD AOLVILINI IVILINI 
9°O =OILVY LNSA10S 


9,° SL =3saNLVasAdwal 


SLS6°O L Ee ¢° 98 LV 8e Oct 
esos6'O (6 SiG S° O08 9° Sz SOL 
v8sr6E oO O° st jh Ae, Sis O6 
91 v6°O O° 9k VEL SY GS) GL 
Gve6 O pot v°6S G& SS 09 
veces O 8G €° SV iS) Sy SV 
LEE O Gas | SE 6 Ve O€ 
9668 °0 Ob 6°9t + 6b St 
‘we ee AMLAWIAVED 2 NIW 
LN30d4d LN30d4d NI 
ALISN3Q ALISOOSIA NOISYSANOD NOTS&YSANO9 SWIt 


bp =ON LNSWI1dddx4 


S poneca we STS gig? martes * 1A TE 


242 


T/OWSO’O =NOILVAYLNSONGO AOLVILINI IILINI 
p°O =OILVH LNAATOS 


0,°S9 =3sYaNLVaAdWsl 


ETI6'O = CLS b’Ov O00 
€1G6°O c° Ot | 8Y vy OV OLt 
tcve'O SSIs 8° BE Vise Svt 
OOE6 'O Be b CE G°6C Oct 
99¢6°0 O's 8° GC Sve OO} 
0O0c6'°O Gro a a Eb? 08 
Et16°O Oe) eG Om Er OS) 
EL06°O (one) Ett 8° EL Sv 
+706 O ee) ie vb b°6 O€ 
9L68°0 0°O O°¢ Oe Gt 
ane a AMLAWTAVD ee NIW 
LN30da3d LN30u4d NI 
ALISN3G ALISOOSIA NOTSYSANOO NOISa4SANOO SWIL 


GS =ON INSWIYSdxS 


a | 


si 


: an 


S23 me & 


+ 


tip qv) (0s 


3 


Lies 


teh e 
Sur A 


») B23 


aed 


t) 


ag 


rac 


7 
~ 
v7 


en a. 


243 


T/1OWSO"O =NOILVYLNAONOO YOLVILINI IWILINI 
b°O =OILVA LN3A10S 


0, °S9 =JUNLVasdWsl 


S986°0 = 6 IES) = Ook 
VEG SO = 8°19 6 ES Ort 
tv96°O O° S6I iy OE hy Oct 
cvae Oo 0°¢9 bk vv v'OVv OOt 
8crE'O L° Ce G SE L es 08 
6LE6°O bck E°Sc v' Le 09 
00¢c6°O Sy i Sy Sark OV 
€806'°0O b°? G9 G*Ot Oc 
we oe. AMLAWIAVUD 2 NIW 
LN30dAd LN30d4d NI 
ALISN3AG ALISOOSIA NOTSYSANO9 NOTSASANOO SWIL 


9 =ON INS3WId4ddxX4 


“_” atewa: anarigder 
> on — = - ae 
> = SeRaTEG TARY: 


OS <PoreRnreyaie aoTAr ree? set Fiat: 


SS! 


244 


1/1OWGO°O =NOILVYLNAONOOD YOLVILINI IVILINI 


v°O =OILVa LN3A10S 


9,°GL =34yNLVaAdW4al 

O886'°0O = Gis = SO} 
Ovl6'O = 8°19 = O06 
9E96°O G Let c GS Cc 8Y GL 
L£4+S6°0 O°SS QV G° Cv 09 
SYLETS) vot CG €€ 8° CE Sv 
OETEB 'O G9 G°O0% O° ES O€ 
€806'°0O (are 8°6 Get St 

Sue ae 203 AMLAWIAVUD 2 NIW 
LN30d4d LiNd9ddd NI 

ALISNAQ ALISOOSIA NOISYSANOO NOISY3SANO9 SWI 


Z =ON INAWITaSdx4S 


; 


= 


245 


T/1OW GO°O =NOTLVYLNSONOO YOLVILINI WWILINI 


9°O =OI1LVe LN3A10S 


o,° G9 =JaNLVaadWw3al 


cst “OOOT | 
cL’ i meleleraa! 
OL’ tl “OOOLTEt 
69°14 “OOOCEt 
S7uit “OOOOEt 
plot “OOOVEL 
OL’ + “OOO6Et 
a “OOOGrt 
plot “OOO8T I 


L° 6k 


9° 


qui /6 


ALISYHAdSIGA1IOd MW SAV LM MW DAV WAN ALISNAG 


8 


=ON LNAWIYadxX4 


do 
ALISOOSIA 


LNAO add 
NOITSYSAND9 


me = 


>. te -Sornagae 


4g arte pele mf 
- = S 


1 ARTMS HAT GOSAET HL “faa Te 


‘~ PTA. 


246 


T/1OW OF °O =NOILVALNSONOO YOLVILINI IWILINI 
9°O =OI1Va LN3A10S 


9,°S9 =daNLVasdWal 


GOS OOEES OO90V LOS6°0O G° O09 Susy OLE 
00°? OO6ER OO6I EEvPE' O O° Cr Saiz Ove 
cet 00¢98 OO8Pr cLEB'O Ly OG Z ES Ol? 
O8't 00888 OOC6Y EOCE6 'O v'6t LE eS Ost 
vet 00688 OOLSY Lvc6'O O'El O° SY Ost 
6L°t OO6LE OO8rs G916°O OE 6° 9E Oc 
Le°t “OOCLE “OO6+1S c606 0 BV 9° 8C “O06 
6L° I “OOO90+t “OO68S tbO6'O 66°C 6° St AGE) 
cok “OOOSCt “OO69SL yO68°O 974 O°6 “OS 
een ae fe ie u/s dSSS*~*«iN NTN 
ALISHYSdSIGAIGd MW’ SAV LM MW DAV WON ALISNAQG ALISOOSIA NOISaYSANOO awit 


6 =ON LNAWIdadxX4 


eae 75 : : an ; z 
— ee ae as 


A" BRO 


_ ote wank me 


css a 


aon ire 


wrens ee eee J 


L.. “he S€ 


ie WH AAI = 


ewes aay 


WANT REIEMND Ut elTine saree 


— ee 


247 


T/1OW GSO°O =NOILVALNZONOO YOLVILINI IVILINI 
9°O =OILVa LN3SATOS 


0, OL =3aNLVasdW3al 


Ost “COOSE6 “OCO8cS LE€v6'O O°SS EAS “OSt 
pat “OOOTOt “OOS6S LLE6°O L°8€ S ZS “O9t 
CLS “OOOLO} “OO C9 ECEB O 88°C VES “Ot 
Bso°t “OOOO! | “OO0CS9 cSc6'O bv O%@ O° 8P “Ot 
99°F “OOOTL + “00699 LOcE'O pv Ot Savy “GOt 
QE “VBCCht “b98E9 9e+46°O O°Ok 6 SE “08 
ELI “O8S6ttt “CESS 1SO06°O Ses GS’ 8c nee) 
Bot OOOE St “OOVEL tL68°0 vE L’6t “OV 
pool “OOCO9TI “OOrSL 9888 °0 Bot 8° Ol “Od 
carey ee oe ON ae w/B = d—SSs~*~*«iNG SON 
ALISH3SdSTIGAIOd MW’ SAV LMA MW SAV WAN ALISNAG ALISOOSIA NOTSYSANOO AWIL 


OL=ON LNAWIYSdx4 


etn ob ee ee ee 


~ - - ee? Pa eres ere Penne ee = 


ed . VT RERTNE Nae 5A ae ee parm rt * ; * ataganeis a ie 
cae “S a = as a Ole nance edit se es e 7 Fe inure Eee = 
sts ia : <= ee . “ay i = 
a a . a = - 
i == - = = = 
sn. eS is 5 a . ; zs ae - ‘Bt = ive =f 
a =e) se 7 qi mM ye ie Ps... o 7 a co 
aoe =< or: _. eee Eeee B-8 = es 8 
a ae 1 a 135 Lae sae 5. ¢ roa oa 
* Co A bee 2 EO! 
/ ey nee 5 asus OS © Js ee: 
i es 


& ee 3 
- 
ae Soe _y Soe tea nts rom oe: 
y = as 
, Tt : 
. y = —_ ~ S P-- — 2 a =— - - 
2 3 8h 2 3S ASST 
3 ar-uitat TeRgg0RF 
D Gy OT TAR, Tea es 
K44-< ) LT AG DEsiHa, BAS Sat areral 
* 


248 


TV/1OW OF’ O =NOILVALNSONOD YOLVILINI WILINI 
9°O =OILVa LN3ATOS 


D9, ° OL =3aNLVYAdW3Al 


98°t OO8c9 OO8EE L956 '°0 GG Syk voEL Ost 
Ost O0O0S19 OO ve LSS6°0 0°09 SL9 Oot 
pst OO0OP9 “OO8TE pere oO O° SE Ve) Ort 
[Silo OO6L9 OOE8E plye oO (ana4 S’°6S Oct 
Le OO0C69 “OOO6E SVE6B'O Sy 7h} 9 eS “OO! 
OL} OOLILL OOETCY 8Sc6'0 G°Ok c° Gv 08 
pL OOLSL OOSEYV 9916°O 8° VLE 09 
OLS OOlSZL OO0O9P 8906 0 8c 6° 0c Ov 
Som OOSS6 OO6LS cS68 0 674 vet Oc 
tee eee een en) ee, cae IN30 Wad NIW 
ALISHYAdSIGA1IOd MW SAV IM MW DAV WAN ALISNAG ALISOOSIA NOTSaY3SANO9 SWI 


tL=ON LNAWIYSdx4 


a . d 4 


rp ae oes ee me eB owe: ena ks cal ed ee ee 66-4 Os VES 66 Aes 
‘Srfessterovses ‘ane ey “= chen sail aq. ‘ ¥Trbsoery seth = 
mg ~* Pia ; ; : 
So =, ° hoe — _ i soil $ 
oe eee a : eRe Ste ees ee oe = side a> eee bd = 9 ene e > 6 = Oe es ee 


| aa 
et”: 
: 
aa 
i 


a a a 


eg Gorge fe ae 


wa Ne 
ef 
‘ “ 
~ 
, 
i 


5 
e 
i 

2 
5 
3 
a 
& 


' %° pr S3Rippensaese 
> @ -prrent rash apy 


‘ty cr ®& sagLeiadetos Seria - sare ied. 


249 


T/1OW GO'O =NOILVYLNAONOO SOLVILINI IWILINI 


9°O =OI1Va LN3A10S 


9, OL =3YyNLVaAdWSL 


68° 1 O0OGl8 
68°} OOO! 8 
L8°b OOC68 
Qe |} “OO6TE 
8° t “OO686 
S/S AY “OOOEOt 
ost “OOO6O0t 
bO°4 “OOO9t + 
vot “OOOVEt 


OOLEV 9G76°0 
O06 CL 68E6 O 
OOGLY GTE6B'O 
OOLLS LvV~6O 
OOL’S EGlI6 O 
O00¢8S 8806 0 
00309 Oc06'O 
OOOTL yS68'O 
OOEO8 1888°0O 
qu /6 


ALISYAdSIGAIOd MW SAV IM MW SAY WON ALISNAQG 


clL=ON LNAWITYSdx4 


do 
ALISOOSIA 


LN39 d3ad 
NOISY3SANO9 


SWIL 


2%, -Wuranewes 
29 -orfar We 


, — : 
j* Pry ua | LARA Ings eae MARTI 


250 


T/1OW OF °O =NOILVALNAONOD AOLVILINI IWILINT 


9°O =OILVa LNSAIOS 


Oo, OL =3eaNLVaAdW3al 


g°O¢ 


Sot 


est 00819 OOLEE éLG6'0O 
es8°t OO8Z9 OOTVE O€G6'O 
pst OO8VT9 OOCSE SOG6'°O 
Ost 10S") 74°) “OO6SE eve oO 
6L I 00049 “OOE9E 09¢c6°0 
Ost OOS9EL OOCOV 09c6'O 
bLo4 OOEEL “OO8CY DMA (0) 
vL't OOEOS “OOL SY 1906'°0 
69+ OOLG8 OOLOS OG68 0 
ao i ree ee er ate ee ww/6SSstidSs*~<“‘é‘NO ONE 


ALISHAdSIGA1IOd MN SAV LM MW SAV WON ALISN3AQG 


do 
ALISOOSIA 


IN39 Yad 
NOISdYSANOO 


NIW 
AWIL 


E€+=ON LNAWIYSdxX4 


| — e- 


(aes BAG = eee ss = i 


See, See 29% ” ie 


yrs a = a 4 . = ’ 
des eee +a : a se 


_.” 3" ai =a. adanrer 
| ~OLFneG lee d 


1% oF Swarr SO Bares Peds? Waee gas 


a 


208 


T/1OW GO°O =NOILVALNZONOD AOLVILINI IILINT 


9°O =OTILVa LN3SA10S 


9, °GL =3ayNLVAsdWw3al 


9°0¢ 


9° Tt 


96°41 OOOES “OOL CE EtG6 O 
be 4 00689 OOOSE O00S6°O 
88°t 00879 OOrre LEVv6 O 
6L OOvsL OO8EV OcvE'O 
SLL OOrls OOLGY E0E6 ' O 
BSL >t OOLC8 OOG9PV vySc6 'O 
BL + OOL98 “OO88P LEGO 
pL t OO8E6 OO0OvS ycO6 O 
O09} OOOSOt “OOSS9 8l68°0 
Sag ge ne a det epee et ww/6—SdSs*~<“<«*“‘«‘i NSC 


ALISHYAdSIGA1O0d MW SAV’ LM MW SAV WON ALISNAG 


do 
ALISOOSIA 


LN30 ddd 
NOTSYSANOO 


vl=ON LNSWIYd4dxX4 


ee 


—- - Ex —: - 5 a nn Lee ant ie ee a4 te - 
a epee We 2a We et ak eTlomotty ee Sate 
=a Berar 2. as 
- ~ 7 : = ee oe, oes ee ae Pe 
: | : : és st a 
ie Se, : ae > ate na | 
ay Sea coe tes 7 y a f 
aa ae me. 3 e ce 
- a 4 © =u 4 
Fe = nr 
7 a reat oe shaky 2.46 % ae 
; #8: {oy et bec * 26 ‘ a4 
£ COE comeegs ee oe %) y : 


— 


oF, =tals7 Og 


ee 7 te eA By 


> 


rey. 2 realy | ee aS aetartnar warried 


Zon 


T/1OW OF °O =NOILVALNAZONOO YAOLVILINI IILINI 
9°O =OILVa LN3SATOS 


9, GL =3aNLVasadWal 


set OOGLG OOGOE S6S6°0 O° O€ O'SL Oct 

EB i OOVSss OO6ILE +cS6 0 O°e? 8° UL GOt 

Ost OCOL8S OO9TE 9EvVE'O L° St io) JAS) O6 

Ost OOCr9 OO8GE O9E6'O O°EL 6°09 GL 

6L I OOOES O01 GE 9SEV6E'O y's g GS 09 

OL’ tL OOE69 OOLOY O616°O SS) eV Sv 

vl t OOVEL OOOLV vit6°O G GCE O€ 

69° OOcSL OOGTY E668 0 Siar 0° Od St 
Pe ee ae ae qw/B = dSSS*S*CNZO. BSN 
ALISHYSdSIGAIOd MW SAV LM MW DAV WON ALISNAG ALISOOSIA NOITS&3ANO9 SWIL 


Gt=ON LNSAWIYddx4 


ges | 


a som Y aT 1) SOrersaiges 


a : eal wig 


ee 2 Sm ENS <P Sd = ee ee wn «9 Le ae eo we 


‘ei one. Ow bacilli oil vies: erraaeiv ene 


Sogn Rem ae ee eat ye 
- oe : - : ~~ e 

as 

Py 

a 

Le 

oot 


SS Prt Ae sae 


ithe reaps 


£ & 
x 


eee tease 


3 


RHINE: 
Nate en 
A 


Whey 


te 


P 


waive tse pe) 
2) 


